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CXXXVIIT. Concentration Contours in Grain Boundary Diffusion 


By R. T. P. Wurerie 
Atomic Energy Research Establishment, Harwell, Berks.* 


[Received June 10, 1954] 


ABSTRACT 


Formulae are obtained for the concentration in a poorly diffusing half 
space bisected by a thin well-diffusing slab, at different times after the 
edge of the half-space has been suddenly raised from zero to unit con- 
centration. 

This represents a simplified model for studying ‘ grain boundary ’ 
diffusion of one metal into another. The problem also has obvious 
application to the study of heat flow in a poorly conducting solid with a 
metallic fin. 


§ 1. INTRODUCTION 


Ir two different metals, A and B are in contact along a plane surface 
(the surface y=0) and are then heated, it is found that the metals diffuse 
into one another. Barnes (1950) has developed a technique by means of 
which concentration contours may be made visible: the penetration can 
also be studied by ‘tracer’ technique. It is found that the penetration 
takes place much more readily along the ‘ grain boundaries’ than else- 
where. It is the object of this paper to investigate this problem mathe- 
matically. 

In the idealized problem it is assumed that the half space y>0 is filled 
with a material of diffusivity D, except for a thin slab of width 2a parallel to 
the y axis, which has diffusivity D’ (D’>D). At times ?<0 the concen- 
tration C of the other metal is everywhere zero. At time t=0 the 
concentration on the plane y=0 is suddenly raised to unity and maintained 
at this value. The problem is to find the concentration C=C(z, y, t) 
elsewhere. 


In the slab C=C’ satisfies z 
Co 


DNC = Ae ae (1.1) 
Outside the slab 
pvc= ee ere te cs 411.2) 
at 

The boundary conditions at the edges of the slab (v= -:a) are 

CEU. es MEN Sie eo iret (L3) 
, 00" oC 

D i Sab) aa (1.4) 


* Communicated by the Author. 
4 L2 
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When y=0 
C=C(a, 0, t)=H(t) 0 ta) oor se 


where H(t) is the Heaviside unit function. 
In actual fact the D’s are functions of the concentration C, so that 
the differential eqns. (1.1) and (1.2) are non-linear. It is not proposed to 


consider this complication here. 
The problem as stated above is obviously equivalent to a heat flow 


problem, where a thin well-conducting fin projects into a block of poorly ~ 
conducting material.. In this case the equations would be exact. 

The same problem has been considered by Fisher (1951), but his solution 
appears to be incorrect. Roe (1951) has also done some calculations on this 
problem. He considers the case of a number of parallel grain boundaries. 
His calculations have not as yet been published. 


§ 2, APPROXIMATION WITHIN THE SLAB 


C’ is evidently an even function of x. Hence if y><a it is a fair approxi- 
mation to take 


2 : 
C'(w, y, )=Co'y, t) + 5 04! (ye >. 


or) ¢ , 


Then (1.1) gives 


oy? 
Neglecting terms of order a”, the boundary conditions (1.3) and (1.4) gives 
C=C,’ 
oC pe 
D Ox =D aC, . 


Hence eliminating C,' and C,' 
D' eC DC oe 
oy2 "adx oat ° 
This represents a single homogeneous boundary condition satisfied by C 
when x=a. Using (1.2) it can also be written 


2, t 
pve _ Da é je 


(2.2) 


7 (2.2.1) 


§ 3. Exact SoLuTION By MEANS OF A FoURIER-LAPLACE TRANSFORM 
Let 


W(x, uw, A)= I, exp (—At) ar{ sin pyO(x,y,t)dy. . . (8.1) 
0 
Then from the theory of Fourier and Laplace transforms 
CG > 12 i i} 4.00 
ey, t)= -{ ’ sin wy du 5 hey exp (At)b(x,,A)dA . . (3.2) 


when the A path of integration is to the right of all singularities of (A). 
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Multiplying eqn. (1.2) by exp (—Aé) sin py and integration with respect 
to ¢ and y 


a2 r co , 00 
= — (w+ 5) p+ I. exp (—At) dt sin BY oy COS MY o| 


y=0 
1 co 
ee) sin py dy [exp (—At) C],-, =0. 
But on using the boundary condition (1.5) 


: oC © 
sin LY = HE COs wy o| =U (ae, 0, 2)==uH (¢), 
Y y=0 


While [exp (—AZ) C}? 5 =—C(a, y, 0)=0. 
Hence Os r im 
Fa (+5) b=—§. PAC a DrPRy at at (3.3) 
The boundary condition (2.2.1) transforms into 
VO eee Os ey) 
D Ox aoe a Ox aaar (F —1)\y=0 ° . . ° (3.4) 
when x=a. 

The solution of (3.3), consistent with (3.4), and having % finite when 
ear ae p=dy+ibs, eee eae ies 87% (3-5) 
where 

b= 1GELAD)’ eae OG) 


oF (D'/D— 1) wexp| — V(r? at > “ie 
~ eB lees) N59 


At this stage it is convenient to introduce the following dimensionless 
parameters 
é= 


(3.7) 


ico a 
eee ev Di 
a 
= TB 3.) 
A=D'/D, 
g=(4—1) a. J 
Then applying (3.2) and replacing w by u/1/Di, A by A/t 
C=C,40C,, ; ee ee (38) 
(= expaA da 
where Ca al. sin (7) a (2A) = 


A—1/° A—V p2-+d €) dr 
Sere 21 )) (4042 Vie Pata) 


(3.11) 
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3.1. Evaluation of Cy 
C,, which is independent of &, evidently represents the concentration 
which would be found in the absence of the slab. 
In fact, performing the A integration 


1 ri ‘ dn 


1” 0 
=1—-| ary | exp (—p”) cos py’ du 
0 — co 


=—| exp (—n’2/4) dy’=erfe 7/2. . . . (3.12) 
n 


§ 4. TRANSFORMATION OF INTEGRAL FOR Cy 


Although not directly expressible in terms of elementary functions, the 
integral for C, can be transformed into a form which is more convenient 
for Smet SUS Te 

Let A=v?—p?. The path of integration for A can be taken from 
co exp (—iz) to oo exp (iz), passing to the right of all singularities of the 
integrand. Hence keeping p fixed, the path of integration for v is from 
co exp (—?m/2) to o exp (i7/2), keeping to the right of all singularities. 
Thus 7 


Oma [ sin yn exp (—p") ) udu f aE 
—ino pet = Vv 


+ 
Noting that gp. 4-1 


(4.1) 


. ld : : 
Begin a as (exp [tuq]+-exp [—i9]) 


—— [exp (—vé’) dé’ 
tits iJ, er [- (erat ae i) a 


one can write 


ra) ie) : co 1? 
O=—[ dé iy do.-| exp (—p? o-+i pn) du 


<font (4) Ju 


Noting that the last ee al vanishes if o> 4, and performing the third 


and fourth integrations 
o—l\2 
a+) J 


Cm F anf a | el gasew[— 7 {E+ 
= ste esp —n?*/40) )enfe| 5 (a= a= (#4 = - (4.2) 


which is in a form suitable for numerical integration. 
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Asa check, suppose 4 — oo, whilst c.=a/\/Di=8/(4—1) remains finite. 
Then for finite o the last factor is simply erfe €/2._ For large values of o 
the last factor is unity, and the integral converges owing to the factor 
go... Hence 


E of] “S 5 do 
asa ny aye ; exp (—72/40) 45 


Oo 
Meee) 
=erfe 5 erf5 
i.e, 
C=0,+0,=1—erf$ orf? 


Thus C=1 when £=0 or 7=0, i.e., x=a or y=0, and C0 when 
E> co and n> 0, ie., > 0 and y-> o. This is clearly correct, 
since if 4 -> oo whilst « remains finite, the slab can be regarded as a perfect 
diffuser, and hence attains unit concentration immediately. 

If 4 o, whilst B=« (4—1) remains finite, (4.2) simplifies to 


2ae Taye 
Cy= 7 j 5378 CXP (—7?/4<) erfe 5 (+ +6) 


= 1.) exp(—™)erte| 5(ag— +8) |on veo ( 43) 


Under certain conditions the integral (4.2) can be approximated by the 
method of steepest descents. It is found, however, that this method is 
only applicable in regions where the concentration is already very small. 
The formulae are given in the next section, and the analysis, which is rather 
tedious, is given in the Appendix. 


§ 5. SumMaRY AND RESULTS 


We have considered an idealized problem in which a slab of width 2a 
and diffusion constant D’ bisects the half space y>0, the rest of which 
consists of material of diffusion constant D. 

It is found that if the concentration on the plane x=0 is suddenly 
raised to unity at time t=0 and maintained at this value, the concentration 
in the quadrant y>0, x>a at time tis a function of the four dimensionless 
variables €, 7, 8, 4 


and is given explicitly by the formula 


4 do 7”? ] A—1 ee! ) 4 
cmertont sa Sexn(—Z)ertos Fa B +& eee (d)) 
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Fig. 1 
B=0-1 


Concentration contours for B=1-0. 
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Concentration contours for B=10-0. 


Fig. 4 


Owens 40 150 Oe 8 FOO 
h. tio of 
i i Graph shows ratio 0 
Penetration of C—0-6 contour along grain boundary. 
; penetration along grain boundary to penetration far from boundary as 


function of B=(4—1) te : 
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Unless 7 is comparable with 4 this simplifies to 


n ae. ae Ke 1 fo—1 
C=erfe n+ sel =373 CP (- r) erie 5 (+ +é 
o 2-2 ] l l 
71 5 (ee eee 
=erfc 4 Je I. exp( is ) entes (ea +6 )dr. . (5.2) 
and this integral is convenient for numerical calculations. 
If € is rather small, and 7>1 but < 43/8 the integral may be approxi- 
mated by the ‘ method of steepest descents ’, and we obtain 
Cy221-159 BU3 n-2/3 exp —[0-473 B-2/3 44/3 0-396 B-4/8 428 (1—Bé) 
+0(8-? (1—BEé)?)]- sop Se Me go vee ob on ee 
To the same degree of approximation the gradient of the concentration 
contours is 


(=) =tan 0 —0-63 79 B-2/3 (10-421 (n8)-2 (I1—Bé)). (5.4) 


These formulae are only valid, however, in regions where the concentration 
is already very small. 


Fig. 5 
65° 
2 
1.8 bw 
1.6 
1.4 55° 
| 1.2 50° | 
@D 
Son © 
Se 0.8 40° 
oO 20° 
0.2 10° 


=P Asse 
B=( HS 7iDe) cae 


Gradient of C=0-6 contour at point of contact with grain boundary 
as function of B 


The formula (5.2) does not involve A explicitly, but only through 
B=(4—1) a//Di. By plotting a series of contours C—constant in the €, 
7 plane for various values of f it is possible to follow the history of a given 
system at different times. This follows since at a given time 


p=(4— 1) = =constant, 


t—a=V Dt &, 
y=V Dt q; 
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so that at later times B diminishes, and the diagrams may be regarded as 
maps of the actual concentration contours on a diminishing scale. 

In figs. 1, 2, 3 the contours C=0-2, 0-4, 0-6, 0-8, 1:0 have been plotted in 
the €,7 plane for B=0-1, 1, 10.* 

Figure 4 is a graph of the ratio of the ordinates (7(0) and 7(0o)) for the 
C=0-6 contour where it meets the slab and far from the slab, plotted as 
a function of 8. It illustrates the extra penetration along the grain 
boundary. 

Figure 5 is a graph of 0, the angle at which the contour C=0-6 meets the 
» axis, plotted as a function of 8. (The contour C=0-6 has been chosen as 
this is the contour rendered visible by Barnes’ technique.) 


APPENDIX 


APPLICATION OF METHOD OF STEEPEST DESCENTS 


Under certain conditions the integral (4.2) can be approximated by the 
method of steepest descents. 
If X is large, 

1 exp (—X?) 
Te OE 
Hence in (4.2), if the argument of the error function is large for those 
values of o which make the most important contribution to the integral as 
a whole 


y fee) 
erfe X= Fe i exp (—n?) dn (A.1) 


A 
Cy | gle) exp [—fle)] do, ee eh ALS) 
where 
Ags! Clerc Blog 3 ee As) 
g(c) TIT 
lfy? d—lfo—l 2 
plo= 7h 4+ 5 3 +8) }- Pee st (As) 


Now 72/o is a monotonically decreasing function ore in’(1, 4). When 
A=" fi@==il , d 
o=1 its gradient is —7”. Also = (= +2) is a monotonically 
TO Sa eR AY 
increasing function of cin (1, 4). Wheno= 1 its gradient is § 3 ae tah 


Hence provided 


2 
vt>t(g+ gai): ep ae GAS) 


* The asymptotic formulae were not used in computing for ake igure 
and in fact it was found that they would only be useful for vasues 0 
considerably less than 0-1. 
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f(c) has a minimum in the range (1, 4), say at oop. Then if this 
minimum is ‘ sufficiently sharp ’ 
C= V 20 go) {f" (o9)} 2 exp[—floo)]. - - + (A.6) 

It is not easy to state in precise terms the conditions under which the 
steepest descents formula is valid: it depends on the variation of g(c). 
However, a fairly good criterion is that 

FOoy ee 

Ne a Bes 

Foe” 9 Sex 


In addition it is necessary that 


X(o0)= 5/4 (“> +£) > Sf eee 


before the approximation (A.1) is valid. 
To the same degree of accuracy as (A.6) 


oC 0 1 4—1 —1 
a re: eal S +8) Cp cae) 
reac oe ei 5) ees . . . . . . . . (A.10) 


so that the gradient of the contours of equal concentration is given by 
fen ee ale eee (oF +)", sal Aa 


d= > 00,jen ~ | Ace B n 


A complete discussion of all possible cases would be rather involved. For 
simplicity we confine ourselves to the practically important case when 4 
is large but B=(4—1) «=(4—1) a/1/Dt is of order unity. It is easy to 


see that if 
n<max ( As2B-1. 4 (3)) 


the minimum of f(c) occurs for values of o<A, and in this case the 
factor (4—1)/(4—c) occurring in (A.4) and (A.3) may be replaced by 
unity, without appreciable loss of accuracy, i.e., 


hares ery See age oe CANT, 
foy= 512+ (G+) }- ot AO AETS) 
Thus f'(o,)=0 


; 292 
if oo? (o9—1+Bé)= 1. 
In particular, if 


[Bol < Reso. > eee nr Atay 


2/3 
y( 7) a Ves j em es, ae \al 
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In substituting in g(o), f (oo) and f’’’(o)) we may neglect the second 
term in oy. Then 


7B 95/6 
g(o9)= Bao Pe pape es Oe) a ee te te) ste (Al 16) 
3 i 
floo)= seq P21? — sz 8-4 9%? (1—BE)-+0 B-* (1—BE).. (A.17) 
Applying formula (A.6) 
ee eet Bee a OXP | —f (04) |a me . =  (A.L8) 
To the same degree of accuracy 
oc 1 1 2 2/3 
oP 5 5 (001-8) Ox — sent #9 (1 5X2)" 1-88) Co 
area (A.19) 
oC 1 n i 1 9\2/3 
on ee a mee PoE peas (1- 4S) (188) Ce : (A.20) 


If C,>(C,, then the gradient of the concentration contours is 
Leos 0C',/0€ 
dC/dn~ ~—- AC, 0m 


=— saan? (.—53(2)" 0-20). per ae Acoh) 


This implies that in the region where the above approximations hold the 
steepness of the concentration contours gradually increases with 7, i.e., 
with y. 

Now consider the range of validity of the above approximations. From 
(A.5) the method of steepest descents is certainly not valid unless 


o> /#(3+ 75): Meee) fe: 1 (2) 


We have restricted ourselves to the case 


n<max (4° 2 5 4()") A he Se oe UD) 


and from (5.15) we have further restricted ourselves to the case 


pp Toe. aS eet em a. (C) 


tan G== 


The approximation to the error function is not valid unless 


la, 1 
Xo 53 Dg’? BOP>1, 


i.e., n> BU? < ata. wea: mantel) 
Now to the same order of approximation 

f(o9)=3X (oo)* 
so this condition also implies that C, itself is small. 
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Combining conditions (c) and (d), it follows that if £=0 we must have 


y>max (5 sr ) 


so that in any case 7>1. In this case C,>C,, since 
1_ exp (—7°/4) 
Va 3] 
for large 7 and this tends to zero more rapidly with 7 than 
exp (—f(o)) exp (—0-473 8-2/9 4%). 
The steepest descents ee is ee valid if 
9f'"oy)®_, 27 | 
2 f'"(o Sao Bl (nB)*°=3f (09) > 1 je a SUE aes 
and this is essentially equivalent to (d). 

Since, by condition (d), C,>C, in the region of validity of the approxi- 
mations, we are justified in replacing C by C, in the formula (A.21) for 
tan 6. 

Summarizing, the formula (A.18), (A.19), (4.20), (A.21) are only valid 
if € is rather small, 7>1 but « 4%/8, and under these conditions C,>(,, 
but <1. In order to investigate the region where C, is comparable with 
unity, numerical integration of (4.2) or (4.3) is required. 

Using the approximate boundary condition (2.2), Fisher derives the 
asymptotic formula 


C,=erfo3 = 


1 
~ exP{ an cm anf erfe €/2. 


As stated above, this formula appears to be incorrect. 
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CXXXIX. The Physical Nature of a Metal Surface in Conduction Theory 


By H. A. Mtszr* 
British Council Fellow at H. H. Wills Physical Laboratory, University of Bristol} 


[Received July 20, 1954] 


SUMMARY 


A metal electron striking the surface of the metal from inside will be 
teflected into the metal. In the present paper this process of reflection 
will be treated as a phenomenon of electron diffraction. The atoms in a 
metal surface form a regular array, and this array may be considered as a 
cross-grating. As all free electrons in a metal have practically the same 
(Fermi) energy, their wavelength is uniform, and the ratio between this 
wavelength and the cube-edge of a face- or body-centred cubic lattice 
can be calculated. This gives a relation between the electron wave- 
length and the dimensions of the cross-grating, and formulae for the 
distribution of intensity over the various orders of diffraction can be 
derived. The main result is that the specular reflection is generally 
suppressed over a wide range of angles of incidence. This is true for any 
angle between the crystal axes and the surface. Generally most of the 
intensity is deflected to the first order of diffraction, but in the case of 
glancing electrons (i.e. electrons striking the surface at very small angles) 
most of the reflected intensity is found at a direction almost perpendicular 
to the surface. 


§ 1. InTRODUCTION 


Recent experimental results on electrons in metals involving the effect of 
the metal surface on the mean free path have been accounted for theoreti- 
cally on the basis of a simple free electron model, in which the reflection of 
electrons at the surface is assumed to be non-specular (Sondheimer 1952). 
This means that the simple idea of electrons in a box made of smooth walls 
(i.e. potential walls) has to be supplemented by the assumption of some 
mechanism which deflects the electrons. Tt is recognized of course, 
that even the best polish technically possible may leave a sub- 
microscopic roughness. Ifthat isso, specular reflection by the plane of the 
macroscopic surface of the metal cannot be expected, but by improving the 
technique of polishing we might hope for a gradual approach to the ideal 
case. The object of. this paper is to suggest a mechanism of electron 


ae 


ee unicated by the Author. ; rag en 
i} ae and ed Physikalisches Institut der Universitat Frankfurt-Main, 


Germany. 
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diffraction, which would lead to the conclusion that specular reflection is 
rendered impossible even in the case of a ‘ perfect ’ polish, i.e., a polish 
leaving only the natural atomic roughness of the surface atoms. 


§ 2, ELECTRON WAVELENGTH AND INTERATOMIC DISTANCE 


For a metal with one electron per atom the free electron model yields for 
the wavelength of an electron having Fermi limiting energy in a cubic 


lattice Sar 
m\ “/ a . . . . . . . 1 . 
c 3n > = | ( ) 


where a is the length of the cube-edge and m the number of atoms per 
elementary cube, ie., n=4 for a face-centred lattice. We focus our 
attention on this case and find numerically 


A= 1-279 a. 


The distance between nearest neighbours being a/1/2, A; is too big to 
allow any diffraction to occur, but only by a factor of less than two. 


Fig. 1 


b. 


Arrangement of atoms in metal surface. 


Even ‘ perfect’ polish as defined above, however, will not produce a 
regular sequence of nearest neighbours in the metal surface. A slab of 
metal will consist of crystallites of random direction, i.e., in a two- 
dimensional picture a plane will look like fig. 1 (b) rather than 1 (a). This 
gives a regular array of atoms in the surface which looks roughly like what 
is known as an echelette grating in infra-red technique ; the grating 


constant is obviously bigger than the distance between nearest neighbours, 
and electron diffraction becomes possible. 
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The resistance of a metal is determined by the scattering of electrons 
with energies very nearly equal to the Fermi limit, and thus with wave- 
lengths very nearly equal to \,. If in the Schrédinger equation the 
energy is a constant, it becomes analogous to the differential equation of 
an optical wave, and the electron waves behave like optical waves in a 
medium of constant refractive index. Therefore we may use the well- 
known optical formulae for echelette gratings, especially for the 
distribution of intensities among the different orders of diffraction. We 
are interested in the question, how much of the intensity is to be expected 
outside the zero order. 


§ 3. QuaLiTaTivE ARGUMENT (Face-CENTRED CuBic LATTICE) 


A qualitative argument will show that there is a fair chance of finding a 
considerable part of the intensity in the first order of diffraction rather 
than in the zero order. 

Let s, d and ¢ define the dimensions of one element of an 
echelette grating as indicated in fig. 1 (6), and let us idealize it as indicated 
in fig. 2. Looking at the left part first, we find the condition that the 


Echelette grating as surface model. 


zero order will be suppressed by the usual procedure used in optics : we 
combine ‘ beam’ | with ‘ beam’ 2, 3 with 4 etc., and find the path length 
difference between any two of them is simply 4=s. The corresponding 
beams will cancel each other, if 4=/2, 3A/2, .. . ete. 

With an echelette grating we find a maximum of intensity in one order, 
if the direction of the diffraction order coincides with the direction of 
specular reflection from the planes of the lamellae, i.e., a direction at the 
angle ¢ from the plane of the grating. The right part of fig. 2 shows that 
the path length difference 4 in the case of this specular reflection is 
A=d sin 24 ; if this is to be a direction of diffraction, 4 must be equal 


to A, 2d,... etc. Now dsin 26=2dsin¢ cos 4, or, as .s=d tan $, 
A=2s cos?¢. Hence the condition for suppression of the zero order is 
EV) Fa ee ROR, RM. eh CRE) 

and for strong intensity in the first order 
A=28 COS? Os es. See, cathe «onl: 0) 


¢ is a small angle, so the two conditions give nearly the same answer. 
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Looking back at fig. 1 (b) for a moment we find that the depth s of the 
furrows of the grating is roughly the distance between two (111)-planes 
for a face-centred cubic lattice, ie., a/,\/3 or 0-577a. Comparing 
s0-577 a and A,=1-279 a, we find our conditions for annihilation of the 
zero order and high intensity in the first order nearly fulfilled ; this state- 
ment is true for a wide range of angles ¢. _¢ must not be too big, of course, 
but if ¢ is big, d becomes small, and if d<A,, there will be no diffraction 
anyhow. 

§ 4. QUANTITATIVE ARGUMENT 


Let us represent the plane surface of a face-centred metal by the 
intersection of an (hkl)-plane with the lattice. If the centre of a certain 
atom is a point of this plane, the same will be true for others, and all the 
atoms in the (hkl)-plane will form a regular array (fig. 3). There are 
various ways of arranging the points into parallelograms (some ways 
indicated in fig. 3) ; if there are only atoms in the corners and none inside 


Fig. 3 


° 


Array of atoms in (hkl)-plane ; proper combination to parallelograms (see text). 


the parallelogram, the areas of all these parallelograms are obviously 
equal. The area can be shown to be 


Fojo/(W+P+LP), . . ... . (Ba) 


if h, k, l are all odd numbers ; otherwise it is 


F=f (i+ k2--B) ee 
Any of our parallelograms may be considered as the basic cell of a cross- 
grating ; we choose the appropriate parallelogram in the following way. 
Starting from an arbitrary atom and calling it atom no. 0, we look for the 
nearest neighbours having a (111)-plane in common with point 0. There 
are two of them; we call them 1 and 1’. Then we choose among the 
atoms of the next row parallel to 1/01 the atom 2 in such a way that the 
angle 102 is as nearly as possible 90°; again there is a corresponding 
atom 2" in the parallel row on the other side. To make the choice 
unambiguous, we may add the condition «(102)>90°, leaving only the 
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identical parallelograms (0132) and (01'3’2’) as basic cells of the cross- 
grating. Figure 4 shows such a parallelogram in the case of a (431)- 
plane chosen arbitrarily as an example. Beside the atoms in the corners 
it shows the atoms inside the parallelogram as seen from the surface, but 
lying slightly below it. 


Basic cell of surface grating, for (431)-plane as an arbitrary example. 


In the (111)-plane the atoms are tightly packed, and their distance is 
smaller than A;,. Consequently, from now on we consider the (111)-planes 
as ‘ optically ’ dense reflecting planes. By so doing we may draw fig. 5 as 
a picture of the two-dimensional grating we have to consider. It is a 
series of flat* kinked steps. We may neglect the kinks, if the steps are 
broad enough, but we shall simplify the problem in a less serious manner : 
we neglect the small blackened area indicated in fig. 5, cutting off a tri- 
angular area at the bottom of the basic cell of our cross-grating and adding 
it on the top. For the resulting arrangement of plane lamellae the 
distribution of intensities among the various orders of diffraction can be 
calculated from well-known optical principles given in any advanced. 
textbook on optics (e.g. Sommerfeld). If the directions of 2- and y-axis 
are chosen as indicated in fig. 5, the z-axis perpendicular to the plane of 
the grating, maxima of intensity will be found, if 


r A r 
A= Ay tN and B=BotM 5 cot CARTE rcs . (4) 


Here a», By, Yo are the direction cosines of the incident radiation, «, 8, y the 
direction cosines for the maxima, A the wavelength of the radiation, a, b, 


le ee 
* Generally they are flat ; it can be shown that there are sub-steps, if (h—k) 
and (k—l) have a common denominator other than the factor 2 occurring if 


h, k, l are all odd numbers. 
4M2 
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as shown in fig. 5, n, and n, any positive or negative integers so that 
a2+fp2<1. Obviously we have 

© yg=— (lag? Bo?) and y= (1a? Bi). 
For the directions given by (4) the intensities become 
J=Cl(y—7o)—(B—Bo) tan ¢)° 


« fot] [Famer eto] 


[oetovoG] (ot Fat oeto v4) 
mee 


C is a constant, given by the condition that the sum of all diffracted 
intensities must be equal to the incident intensity, the other quantities 


Fig. 5 


Idealized surface grating. Blackened area neglected in calculations. 


being defined as before. We see that J will have its maximum value, 
if n»=—(y—ypo) s/A and n,=0, because in this case both factors in the 
denominator become zero. From our qualitative argument we expected 
strong intensity in the first order for A=2s cos? d. Now we find the same 
condition again: in the case of perpendicular incidence, which we 
assumed there, xgj=f)=0, yo=—1; with n,=0 and n,=—1 we have 
from (4): «=0, B=—A/b sins, or, as b sin ys corresponds to our grating 
constant d from fig. 2: B=—A/d; hence y=+ 4/(1—A2/d?). To make 
the denominators of (5) 0, we must have 
A=s (y—yo)=8 [1+ -V(1—A*/d?)]. 

Solving this equation, we find 

2s 

— 1+(s/d)?" 


Since s/d=tan 4, this is exactly condition (2 6). 
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Apart from this simple case eqn. (4) is too complicated to be visualized 
easily. An example may illustrate its results. Figure 6 shows the 
percentage of intensity in the zero order, i.e., in the direction of specular 
reflection, for two values of fy, viz. 0 and 4, as a function of cos”! a, ie., 
the angle towards the x-axis. The values are calculated again for the 
(431)-plane used as an example above. Figure 7 shows the intensity for 
n,=0, %)=—1, under the same conditions. For angles above about 
+20° and £,=0 the condition «?+-6?<1 is no longer fulfilled, the first 
order does not occur. The figures show that specular reflection is strongly 
suppressed over a considerable range of angles of incidence ; specular 
reflection is not the predominant behaviour of electrons striking the metal 
surface even if there is only the natural atomic roughness of the surface, 
and even if the most perfect polish imaginable has been applied. 


§ 5. RANGE OF VALIDITY OF THE APPROXIMATION 


This conclusion is more general than the special example of the (431)- 
plane might suggest. Calculations of other cases show qualitatively the 
same results in spite of quantitative differences. This is to be expected, 
as the quantity s, the depth of the grating furrows, mainly determines the 
conditions (2a) and (26), and this is the same for any (hkl). Generally 
it may be said that the range of angles of incidence, over which specular 
reflection disappears, increases with increasing area of the basic cell of the 
grating. 

Another objection against the calculations might be that the steps with 
their sharp edges are an inadmissible idealization of the actual form of the 
boundary of the metal. This is not too serious, however. In an early 
treatment of diffraction problems Voigt (1911) showed for the case of the 
one-dimensional grating that the exact form of the grating is meaningless 
to a first approximation. He expands the periodic function determining 
the cross section of the grating into a Fourier series, and shows that only 
the first N terms of the series influence the diffraction phenomena if N 
orders of diffraction occur ; the higher orders of the series are negligible 
to a first approximation. In other words, alterations of the form of the 
grating elements can be neglected as long as they remain below the order 
of magnitude of the wavelength. In this form we may use Voigt’s 
argument for the two-dimensional case : as the diameter of the atoms is 
smaller than 4,, it does not matter whether we consider sharp steps or 
some rounded-off structure, which would certainly be nearer to the truth. 

A more serious objection is this. The optical formulae used for the 
calculations are based on Kirchhoff’s treatment of diffraction. This 
implies the assumption that the incident radiation reaches the grating 
undisturbed. This is true only if A<d, in other words our calculations 
may lose their quantitative character in the case of small (Akl), where the 
dimensions of the grating elements are bigger than, but still comparable to 
the wavelength. But even then the calculations will be a fairly good 
approximation, since no major discrepancies between theory and experi- 
ment are known from infra-red technique, where the same formulae have 
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been used successfully for the design of echelette gratings, and the 
condition A<d is by no means fulfilled either. 


Fig. 6 


specular 
on reflection 


[%] 


-80° —*0 ae 20 20 8640 60 80° 


Percentage of intensity in the zero order (specular reflection) for different 
angles of incidence, (431)-plane as an example. 


-60 


~60° 4 40 -20 V 20 40 60° 


Percentage of intensity in the (0; —1) order for different angles 
of incidence ; (431)-plane as an example. 
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§ 6. GLANCING ELECTRONS 


Finally it is a shortcoming of the calculations that they do not include 
the case y)~0, i.e., the so-called * glancing electrons ’, which are the most 
important ones in certain surface phenomena (Pippard 1947). Obviously 
we must take into account the fact that a considerable part of each 
element of a grating is ‘ shadowed ’ by the preceding one (fig. 8). We can 
try to overcome this difficulty by leaving the shadowed parts out of the 
integration leading to formula (4). It must be realized, however, that 


‘ Shadowed ’ parts of grating under glancing incidence. 


this procedure makes the difficulties of the preceding paragraph all the 
more serious. Very narrow beams are cut out of the wave front by each 
grating element now, and the width of such a beam may now easily be 
smaller than \;. In this case diffraction by the edge of the (n—1)st 
element precedes the reflections by the nth element of the grating. 
Bearing this in mind we work out the following formula for the limiting 
case y 0, extending the integration over the parts not shadowed. As 
the result will be only qualitative anyway, we omit the kinks of fig. 5 and 
write down the formula only for the one-dimensional grating with furrows 


of width d: M f iar 
7 YO 
J=C | o—105 —n tan | a2 tan? ‘ 5 0 dj 5 (6) 


As the zero order is now to be found at a glancing angle, only negative n 
are possible in this case. For the case of specular reflection we have 


n=0, y=—ypo- hence 
“A 4 Any 4d? 


Lome? Na»? tan? d 


We obtain a simple expression for J/Jo : 


i 2 d ; 
oral om tan. 6 ° . . ri i 5 (7) 
0 


The predominant term is yd/X: the first order (n=—1) is no longer 
preferred ; there is a strong preference for those orders reflected almost 
perpendicular to the grating (y=1). That means that the glancing 
electrons will preferably be thrown down into the metal at large angles 


from the surface. 
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§ 7. THe Bopy-CENTRED LATTICE 


So far we have dealt with the face-centred cubic lattice. We found in 
§ 2 that the main condition for suppression of specular reflection and 
enhancement of the first order of diffraction was A-=2s. In a body- 
centred cubic metal the (110)-planes are the most densely packed ones ; 
the distance between two consecutive planes is a/\/2=0-707 a. A; now — 
becomes ~1-612 a (see (1)). The ratio A,/s becomes approximately 
1-612/0-7072-28 as compared to 1-279/0-577=2-22 for the face-centred 
case, i.e., again close to the required value of about 2. The area of the 
basic cell (cf. (3)) becomes 


F=he0/(V++P), 
if two of the (Akl) are odd and one even, otherwise 
FH y/(?+h+0). 
Apart from these minor differences and the fact that we have to consider 


the (110) instead of the (111)-planes, all the conclusions remain valid for 
the body-centred lattice. 
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ABSTRACT 


From the.experimental data at present available it is shown that the 
Reynolds number in the upper F region (300-400 km level) is of the order 
of 300. The region may therefore be turbulent. It is suggested that the 
high level diffracting screens which give rise to spread F echoes and to 
radio star fading are caused by non-laminar flow, and that their non- 
appearance during the daylight hours may be due to the inhibition of 
turbulence by large temperature gradients, by lower drift velocities, or 
by an increase in the kinematic viscosity. 


$1. INTRODUCTION 

REcENT observations by radio-echo, radio-astronomical and rocket 
methods have shown that in the F region of the ionosphere there 
are drift movements with exceedingly high velocities (Maxwell and Dagg 
1954 a, Briggs and Spencer 1954). As most large-scale fluid motions in 
nature are turbulent it is interesting to enquire whether these movements 
in the uppermost ionosphere may also be turbulent, and whether the 
diffracting screens which constitute sporadic F phenomena (the screens 
which cause radio star fading, for example, and the ‘ cloudy ’ structure 
associated with spread F echoes on h’-f records) are the result of non- 
laminar flow. 

Radio observations of movements in the ionosphere depend on the identi- 
fication of irregularities in its structure and then on measurements of the 
rate at which they move. It is not yet certain whether the movements 
of these irregularities in the ionization represent true movements of the 
atmosphere, whether they are electrodynamic ‘movements of the ionization 
alone, or even whether they should be interpreted simply as a redistribution 
of matter by progressive hydrodynamic waves. A complete investigation 
of the problem of turbulence in the ionosphere will not be possible until 
these questions are resolved. Rocket measurements of air movements in 
the E and F regions and photometric observations of the night sky glow 
suggest, however, that whatever relation may exist between the move- 
ments of irregularities in ionization and the neutral molecules of the 
atmosphere, the latter have systematic motions of comparable magnitude. 

This paper is essentially concerned with (i) an examination of the tur- 
bulence criteria for levels of the ionosphere above the maximum of ioni- 
zation in the F region, and with (ii) turbulence as a possible origin of the 
gas Sie Se ae 2 oe a 


* Communicated by Professor A. C. B. Lovell. 
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irregularities which diffract the radiation from the discrete sources of 
extra-terrestrial radio emission (radio stars); these irregularities are 
believed to be above the level of maximum ionization in the F region and 
are remarkable for their correlation with the appearance of spread F 
echoes on ionospheric sounding records, and for the fact that they appear 
mainly at night. 
§2. CRITERIA OF TURBULENCE 
(a) Reynolds Number at the 300-400 km Level 
The main difficulty in applying the Reynolds criterion for turbulence 
to the F region is the selection of appropriate numerical values for the 
parameters from the imprecise data about the F region which are at 
present available. The Reynolds number, Re, is defined by 
ee prly 
be 
where v=total velocity difference occurring over a depth /), ~=viscosity, 
p=density, u/p—l=kinematic viscosity.* It may first be questioned 
whether the term ‘viscosity’ has any real meaning when referred to the 
uppermost ionosphere, because the mean free paths there are so large. 
Estimates of the mean free paths at the 300-400 km level vary consider- 
ably, as shown in table 1, but a probable value is 200 metres. The 
thickness of the region under consideration is some 500 times this path 
length and it may therefore be assumed that viscous effects will be evident. 


Table 1. F Region Parameters (300-400 km Level) 


> 


Mean Kinematic 
Reference free path viscosity 
(cm) (cm? sec~+) 


Gerson, N.C.,1952, Advances in Geophysics I, 
p. 165; 1951, Rep. Phys. Soc. Progr. Phys., 


14, 316 2x 104 10° 
DE JAGER, C., 1952, Mém. Soc. Roy. Sci. 

Inége, 12, 223 6 x 108 4x 108 
Mirra, 8. K., 1952, The Upper Atmosphere, 

2nd Edition, p. 28 5x 10° 
Yera, D. G., 1952, J.Geophys. Res., 57, 217 6x 105 aI ra 
YerG, D. G., 1951, J. Met., 8, 244 8x 105 3x10” 


The present uncertainties in the temperature, pressure and density data 
for the F region lead to a similar diversity in the estimates of kinematic 


acEO avoid confusing ‘kinematic viscosity’ with ‘collision frequency’, for 
oth of which the symbol v is widely used, it is suggested that the symbol ¢ 
be used for kinematic viscosity in the ionosphere. 
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ao ¢, but the values listed in table 1 suggest that it is of the order 

The length dimension /, appropriate to the Reynolds number might be 
taken as the thickness of the region, but this would imply that the drift 
direction throughout its vertical extent is the same ; as this may not 
be the case, a value corresponding to half the thickness=2 scale 
height=270km will be assumed. Velocities in the region, estimated 
by visual and radio methods, have magnitudes of the order of 200 m/sec. 
These data lead to a value of Re=300. 

There is no universally accepted critical value of the Reynolds number 
for shearing flow in a free atmosphere, but present evidence suggests that 
it is much lower than the classical value of 1000, estimated for the flow 
of liquids through tubes; for example, experiments with jet streams 
(Andrade 1939) give critical values of 30. It would therefore seem highly 
probable that in the absence of strong damping factors and at the Reynolds 
numbers observed in the upper F region turbulent flow will result. 


(b) Microscale of Turbulence at the 300-400 km Level 

The fluctuating velocity field in a turbulent medium may be regarded 
as the result of the superposition of periodic variations with all possible 
wavelengths. These constitute a hierarchy of ‘turbulent eddies ’, the 
maximum wavelengths, /) (macroscale of turbulence) being limited by the 
dimensions of the medium under consideration, and the minimum 
wavelengths, 7 (microscale of turbulence) by viscous dissipation in the 
smallest eddies. In the case of the upper F region, where the Reynolds 
number is only of the order of 100-500, it is not permissible to apply the 
spectral laws (Kolmogoroff 1941,and others) which refer to fully developed 
turbulence, that is Re—> oo, but it is possible to calculate the cut-off in 
the hierarchy of eddies due to viscous damping. This is given by the 


equation c3\ 1/4 

Gis 
where C=kinematic viscosity and «=rate of dissipation of energy=v3/lp, 
that is, £3],\ 1/4 

= & 


For /,=140 km and Re=300 

n=2 km. 
It is interesting to notice that the effective wavelength of the F region 
irregularities which diffract the radiation from extra-terrestrial radio 
sources is of the same order (Maxwell and Little 1952). 


(c) Turbulence at Lower Levels in the Ionosphere 
With decreasing height the kinematic viscosity (the parameter 
essentially controlling turbulence) decreases very considerably, so that 
turbulence should become increasingly evident in the lower levels of the 
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ionosphere. In the lower half of the F region (200-300 km level) for 
example, (=107, J,=100km, v=100 m/sec, whence the Reynolds 
number—10! and the microscale of turbulence=100 metres. In this 
connection it may be remarked that the diffracting screens used for 
radio-echo measurements of drift movements in the lower F region 
(Briggs and Spencer 1954, and others) are always observable—as distinct 
from the high level diffracting screens which in general occur only at night 
—and have wavelengths of the order of 100 metres. 

In the E region (80-120 km level) €=10°, v=50—100 m/sec, /j=30 km 
and the Reynolds number=2 x 10°. It is to be expected therefore that 
the E region will invariably be turbulent, and observations of the decay 
of meteor trails confirm this.* The microscale of turbulence in the 
E region=10 metres. Scattering of electromagnetic radiation by turbu- 
lent irregularities in the E region has been investigated by Villars and 
Weisskopf (1954). 


§3. Rapio Srar FADING 


It is now generally believed that the ionospheric diffracting screen which 
causes fading, or scintillation, in the signals from radio stars occurs above 
the level of maximum ionization in the F region (Little and Maxwell 1951, 
Hewish 1952), but its origin is not understood. The maximum of radio 
star fading occurs near 014, and from fig. 1 it will be seen that this 
coincides in time with a maximum in the occurrence of irregular ionization 
in the F region. Ryle and Hewish (1950) have sought an explanation 
of these phenomena in terms of some ionizing process which has a 
maximum shortly after midnight. They suggested that the additional 
agency might be corpuscular matter of interstellar origin moving under the 
gravitational attraction of the sun. The radial velocity of such particles 
when they intersect the earth’s orbit is 42 km/sec, and as the orbital 
velocity of the earth is 30 km/sec their motion is normal to the surface of 
the earth at 025 20 local time. The mechanism should exhibit an onset at 
20 20, a maximum at 02) 20 and a cut off at 08" 20, and Ryle and Hewish 
pointed out that this is in rough agreement with the known variations in 
the radio star scintillation amplitude and in the occurrence of irregularities 


i It would appear that the Richardson criterion for turbulence, Ri, which was 
originally developed for tropospheric problems and which is defined by 


elf aT’ du\2 
mimo (Etr)/(G) 


is inapplicable to the ionosphere. In the case of the E region the velocity 
gradient du/dz=83 m sec~1 km (Greenhow, J. S., 1954, Phil. Mag., 45, 371), 
the adiabatic lapse rate [=1-5 10-4 deg. kK em=!, the temperature gradient 
dT'/dz=3-2 deg. Kk km~! (The Rocket Panel, 1952, Phys. Rev., 88, 1027), and the 
temperature 7’==250 deg. k, whence the Richardson number—10. As this 
value is much greater than unity it requires that the E region should not be 
turbulent. This is in contradiction to the experimental evidence. It may 
therefore be concluded that, unlike the tropospheric case, turbulence in the 
E region is not chiefly maintained by heat convection and buoyancy processes ; 
this further suggests that vertical movements in the E region may be restricted. 
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in the F region. Ginzburg (1952) has developed this theory in more detail 
and has considered not only the effect of incoming interstellar gas particles, 
diameter 10~8 cm, but also the possible effects of interstellar dust particles, 
diameter 10-?cm; he concludes that the interstellar dust will pass 
unhindered through the F region and that it can therefore be ruled out 
entirely as a generator of sporadic ionization in that region. 

The suggestion that interstellar gas is the ionizing mechanism must, 
however, withstand criticism on the following counts: (i) The inter- 
stellar gas attracted towards the sun comprises mainly protons, and their 
energy when moving at the maximum velocity relative to the earth, 
52 km/sec, is only 14 ev, which is insufficient to ionize the particles of the 
upper atmosphere. Ginzburg has considered the possibility of ionization 
by the helium atoms of the interstellar gas ; their energy is 57 ev, which 
is adequate for ionization, but it is questionable whether their density is 
sufficient for such a mechanism to be completely effective. (ii) It 1s 
difficult to understand why the diffracting screen should be formed for an 
average of only one night in three, since the interstellar matter impinges 
on the atmosphere every night. (iii) It is also difficult to understand why 
the interstellar matter should ionize in irregular patches. (iv) Radio star 
scintillations and sporadic F are often observed well before the predicted 
onset time, 20" 20, and rarely after 06, even though the cut-off time is 
08» 20. 

Observations of radio star fading during the recent eclipse, June 30, 1954 
provide further evidence against an accretion mechanism. It has been 
mentioned above that the diffracting screens causing radio star scintilla- 
tions and spread F echoes normally occur at night ; during these hours 
there is a decrease in the F region critical frequency. At the time of the 
eclipse, 12" 30 local time (=U.T.), the line of sight to the Cassiopeia radio 
star from Jodrell Bank intersected the F region at a point 100 km from 
the area of totality in the region, and the temporary simulation of night 
time conditions in this area caused the F region critical frequency to drop 
from 4:6 to 3-8 Me/s.* On the equipment used at J odrell Bank for con- 
tinuous observations of the Cassiopeia radio star the records showed a 
sudden onset of scintillations at 13210, that is 40 minutes after the 
maximum phase of the eclipse, and these continued for the following 
90 minutes. This appearance of scintillations during an eclipse occurring 
at 12 30 could not be attributed to an accretion theory. 

As an alternative mechanism, it may be suggested that the irregular 
distribution of ionization in the upper F region is caused by turbulence, 
and that the appearance of the diffracting screens represents a transition 


* Inverness values. The mean value of the F region critical frequency at 
noon during June 1954 was 4-6 Mc/s, and it did not fall below 4-2 Mc/s on any 
other day during the month. It will further be noted that the eclipse day was 
geomagnetically quiet (K indices <3) so that the decrease in F region critical 
frequency at the eclipse time could not be attributed to magnetic activity. 
The author is indebted to the Department of Scientific and Industrial Research 
for providing the ionospheric data in advance of publication. 
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from laminar flow during the day to non-laminar flow at night. These 
irregularities may be distributed throughout the whole region, or they may 
be localized near an interface between different circulatory systems in 
the F region. The arguments advanced in §2 show that the Reynolds 
number in the upper F region is, in all probability, above the critical value. 
It therefore becomes necessary to account for the disappearance of the 
smaller scale sporadic phenomena during daylight hours in terms of some 
process which decreases or inhibits turbulent flow during the daytime. 
This diurnal variation could be brought about by three mechanisms acting 
separately or in combination. (i) Repression of turbulence by large 
temperature gradients : if temperature gradients in the F region are much 
larger in the daytime than at night, there will be a greatly increased 
probability of turbulent flow occurring during the night hours. (ii) In- 
crease in velocity : for example, measurements of drift movements by radio 
astronomical techniques show that drift speeds increase significantly at 
night, from an average of 60 m/sec at 19 to 290 m/sec at 01" decreasing 
again to 170 m/sec at 05" (fig. 1). (iii) Decrease in kinematic viscosity : 
the night-time contraction of the ionosphere may lead to a decrease in 
kinematic viscosity but at present there is no experimental information 
about this. Finally it will be noted that turbulent flow will to some extent 
be continually damped by the terrestrial magnetic field. 

The association of sporadic F with geomagnetically disturbed conditions, 
and the occurrence of radio-star scintillations in the day time during severe 
magnetic storms may be accounted for in a similar way. There is a 
fourfold increase in the drift speed of the F region irregularities during the 
most intense magnetic storms (Little and Maxwell 1952, Maxwell and 
Dagg 1954 b) and at the same time the ionosphere is irradiated by excess 
ultra-violet and x-ray emissions from active areas on the sun, or is 
penetrated by fast travelling auroral particles, so that violent changes in 
the kinematic viscosity and temperature gradients may also occur. 


§4. FurtHeR EXPERIMENTAL Ports 

The rapidity of the fading in radio star signals is proportional to the drift 
velocity of the ionospheric irregularities which diffracts the incident 
radiation (Maxwell and Dagg 1954 a). As turbulent effects become more 
pronounced with increasing velocity it would be expected that the 
amplitude of the fading should show some correlation with the fading rate. 
There is no one to one correlation between these two quantities, but there 
is some experimental evidence which qualitatively supports this prediction. 
Bolton, Slee and Stanley (1954) in an investigation of the fluctuating 
signals from the Cygnus source, found that variations of greater amplitude 
seemed to have shorter periods; Ryle and Hewish (1950) occasionally 
observed sudden increases in both amplitude and fading rate which then 
gradually decreased during a subsequent period of 30 to 40 minutes. Mills 
and Thomas (1951) found a correlation of 0-6 between amplitude and 
the average fading rate for the case of fast fading ; at Jodrell Bank 
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sudden and simultaneous increases in both amplitude and fading rate 
have often been observed, and it has been especially noticeable that the 
fading rate falls markedly as the fading amplitude decreases near dawn. 
It is the nature of turbulent flow, that the eddies are constantly forming, 
disintegrating and reforming; this process goes on spontaneously and 
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continuously. The decay time 7 of an eddy of wavelength, /, is given by 
l 


T= 


> 


where u—the turbulent velocity. The contours of the radio star diffrac- 
tion pattern have wavelengths of 3-5 km * and are moving at speeds of 
200 m/sec. Assuming they are caused by turbulence, and that they have a 
turbulent velocity of 100-200m/sec, these contours should have a decay time 
of 20seconds and in consequence should be individually distinguishable only 
over paths of the order of their characteristic wavelength. Experimentally 
this is found to be so. Fading records taken with receiving equipments 
arranged at varying distances along the direction of motion show marked 
changes at spacings of 4 km (due allowance being made for the systematic 
time delays between the records) and no correlation at all at 7km. It 
will also be noted that the decay time of an eddy is directly proportional 
to its characteristic wavelength. If the large-scale F region irregularities 
with fronts of the order of 200 km moving with velocities of 200 m/sec 
(observed by Munro 1948) represent a macroscale of turbulence in the 
region, they would be expected to retain their identity over paths of 
several hundred km. Experimentally this is again found to be the case. 
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ABSTRACT 


It is shown that Hartree’s equations for the negative hydrogen ion 
have several solutions, the usual one of which does not minimize the 
energy. This non-uniqueness is also shown to exist for a Fermi gas 
where the usual plane waves of the Sommerfeld model do not even give 
a stationary value for the energy. The energy of the gas can be lowered, 
by describing some of the electrons by wave functions that decrease 
exponentially. The energy is lowered by about 1 ev per electron of this 
type. If the electron density is sufficiently low (band width H<0-5 ev) 
all the electrons can be so trapped and the metal becomes an insulator ; 
this provides an a priori reason for Mott’s rejection of the Bloch wave 
functions for nickel oxide. For higher electron densities it is possible 
to trap a fraction of the electrons, so saving considerable energy ; this 
is compared with Wigner’s correlation energy. It is suggested that 
Friedel, in his explanation of the fine structure of an x-ray absorption 
edge, has used several solutions of Hartree’s (or Fock’s) equations that 
correspond to only one exact excited wave function, so invalidating his 
explanation of this effect. 


$1. INTRODUCTION 


THE variation method,‘and ipso facto Hartree’s and Fock’s equations, 
have generally been used without regard to three important points : 
The variation method selects from a given set of wave functions that one 
which, on being varied, gives a zero first order variation in the energy ; 
however : 

(a) The minimum of the energy found need not be the only one. 
There may be several others. 

(b) There is no @ priori reason why the energy should be & minimum 
(as it must be for a ground state) for any given turning points ; this 
depends on its second derivatives. 

(c) There is no reason why the lowest energy should be attained at a 
turning point ; it may be on the boundary of the set of wave functions. 

These points are probably of no importance for the inner electrons of 
atoms, but the outer electrons, especially when these are loosely bound, 
er OE 
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as in negative ions, may very well have several possible wave functions 
and the usual ones need not necessarily give the lowest energy. These 
three points are illustrated in the case of a two electron atom in § 2, 
where we show that for Z near to unity, e.g. for H~, some or all of these 
points need attention. 

In $3 we consider the solution of Hartree’s equations for the Sommer- 
feld model of a metal. In this model the electrons are supposed to move 
in the field of a uniform positive charge of density |e|/2. The usual 
solutions of Hartree’s equations for this problem are plane waves, but 
we shall show that there are other solutions in which electrons are trapped 
in a bound state, e.g. exponentially decreasing wave function from some 
centre. The free electrons decrease their density, 2] |?, in the neigh- 
bourhood of the trapped electrons, so screening them. In consequence 
the Hartree field for a trapped electron is an attractive coulomb one, 
because an electron does not move in its own potential field. This 
immediately implies that it is possible to write down a solution with any 
number of electrons trapped self-consistently. That is to say Hartree’s 
equations for this problem have a large number of solutions, all corres- 
ponding to the ground state (point (a) above). It now becomes essential 
to examine the energies of these solutions (point (6)) and in general 
the energy as a function of the bound wave functions (point (c)). This 
we do in §3 taking only Is states of the form exp (—Ar) for the bound 
electrons and using the Thomas—Fermi method to give the electron 
density and sum of the eigenvalues of the free electrons. 

We find that, no matter what the electron density is, the energy is 
always lowered if one electron is trapped. We can now trap an electron 
at each one of a set of sites in the metal, but we must not put the trapped 
electrons too close together or they will overlap and the energy will rise. 
We suggest that this is a perfectly feasible explanation of correlation 
energy simply in terms of Hartree’s equations and we compare it in 
§ 4 with Wigner’s (1934) calculation of the correlation energy using much 
more complicated wave functions. In our opinion the situation is 
exactly analogous to that considered by Chandrasekhar (1944) with H7: 
that whereas it is possible to get a better value of the energy in the 
variation method by making the wave function much more complicated, 
as Wigner has done, it may also be possible by quite insignificant complica- 
tions in the wave function to do just as well. The advantage of the 
latter process is not merely the mathematical simplicity of handling 
the problem but the greater physical insight which is possible for a 
simple wave function as compared with a conplicated one. Mott (1949) 
has proposed that NiO is an insulator, despite the fact that not all its 
bands are full or empty. He claims that the Bloch wave functions 
(plane waves in the Sommerfeld model) are not invariably the correct 
ones to use, but that in some cases NiO being one, Heitler-London wave 
functions should be used. We suggest in § 5 as an alternative description 
leading to the same result, that trapped electrons (corresponding to 
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Heitler-London wave functions) always exist but that provided the 
band is sufficiently narrow then all the electrons may be trapped (as 
was suggested by Wigner as early as 1934 and further amplified in his 
1938 paper) so producing an insulator. 

Friedel (1952 a) has suggested that the screening of impurities in metals 
may frequently be accomplished by trapping one or more electrons over 
the impurity. He has used this (Friedel 1952 a,b) to explain X-ray 
spectra, the atom of the metal which has lost a K-electron, say, being 
regarded as an impurity. He furthermore suggests that the wave 
function for the trapped electron is not uniquely fixed. Thus in lithium 
an atom minus a K-electron may be screened by a 2s, 2p, 3p... electron 
and each of these excited states have different energies. Friedel suggests 
that the fine structure near the absorption edge of x-ray spectra is caused 
by transitions between the ground state and these different excited states. 
We believe that these excited states are not really different states but are 
merely different approximations to one excited state (point (a)) and so 
Friedel’s explanation of the fine structure breaks down. We discuss 
this further in § 6. 


§2. Two ELEecrron ATOMS 


The simplest possible choice of product wave function for a two 
electron atom is exp (—Ar,) exp (—Ar,). Choosing a value of A that 
gives a stationary value for the energy H, we find H=—(Z—5/16)? and 
A=Z—5/16, in atomic units. We may compare this with the correct 
value of # by subtracting the energy of the corresponding one electron 
atom, to give the first ionization potential. Condon and Shortley 
(1951, p. 346) give a table showing the percentage error of these ionization 
potentials for Z=2,3,and4. Theerrors are 6-4, 2-1 and 1-0% respectively. 
For Z=1, i.e. H’, the calculated value of the ionization potential is 
—0-027 while the observed, or rather accurately calculated value is 
+ 0-028 making an error of 200°. As Chandrasekhar (1944) has pointed 
out, while helium is a stable closed structure, the negative hydrogen ion 
is an open structure, which exists principally on account of incomplete 
screening and polarization. Chandrasekhar has therefore suggested that 
instead of starting with the above wave function we should try 

exp (—Ayry—Agrs) + exp (—Ayra— Agha) 5 
he then finds that A,;=1-039 and A,=0-283 and that the first ionization 
potential is +0-013 so reducing the error for H™ from 200% to 50%. 
Chandrasekhar’s wave function is not of the Hartree product form, but 
this does suggest that we try for discussion the unsymmetrical wave 
function exp (—A,7;—A2"2). In other words we are not going to insist 
that electrons with the same quantum numbers n,/, and m, but opposite 
spin, should have the same space-like wave function. The energy of 
this wave function is easily computed ; it is ; 
AP+BAAe+Ag 
ahaa eA Mat Aa) Ae ea? 


4N2 
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For the turning points of this function we have the equations 0H/0A,=0 


and 0#/0\,=0. We find that \,=A,=Z—5/16 is always a solution of 


these equations (for Z>5/16), but that for Z<17/16 (= 1-062) it is a 
saddle point and only for Z>17/16 is it a mmimum (point (b) of §1). 
Furthermore for Z<1 there are no other turning points, the lowest value 
of the energy being found on the boundary of admissible wave fanchiom 
(point (c)), namely for (Z, 0) or (0, Z), the energy of which is —Z?/2, 
which is less than —(Z—5/16)? for Z in the range (2—/2)5/16 (<5/16) 
to (2+4/2)5/16 (=1-067). We see that even when the solution Aa 
is a minimum it does not necessarily give the lowest energy, namely 
when 1-:062<Z<1-067 (actually there is another minimum for these 
values of Z which is less than —Z?/2). For 1<Z<13/12=1-083 the 
usual solution A,;=A, found above is not the only turning point (point (@)). 

We now give a description of the motion of the turning points about 
the \,A, plane as Z is increased; it will be followed more easily by 
reference to figs. 1 (a), (6), and (c). For 5/16<Z<1 there is only one 
turning point’ (a saddle point) at the position 4;=A,=Z—5/16. At 
Z=1 (fig. 1 (a) two minima appear at (1, 0) and (0, 1) and as Z increases 


they move in towards the line Ay=A,. At Z=17/16 two new saddle 


points appear at (17—5)/16 and as Z increases they move away from the 
line \,=A, leaving behind them a minimum at A,=A,=Z—5/16, where 
there was before a saddle point (fig. 1(6)). At Z=13/12 the saddle 
points and minima coalesce leaving only the minimum at A,=), (fig. 1 (c)). 

This simple example demonstrates very clearly the three points made 
in the introduction. It is true that we have not solved Hartree’s equations 
exactly since we have limited our one-electron wave functions to the 
form exp (—Ar); there seems to be no reason however to suppose that 
the three points would not be confirmed by solving Hartree’s equations 
exactly. 

It may be objected that the solutions of Hartree’s equations are not 
unique and do not give the minimum energy in only very narrow ranges 
of Z, and that no atoms exist with Z in these ranges (except for H-), 
but that does not rob our results of physical significance. On the contrary 
we should say that when two or more electrons are very loosely bound, 
or as Chandrasekhar puts it : when incomplete screening and polarization 
are important effects, then their wave functions are not given uniquely 
by Hartree’s equations. The same situation is probably implied when it 
is found difficult to obtain self-consistent solutions of Hartree’s equations 
due to what Hartree would call ‘ over-stability ’ (Hartree 1933) a quite 
common occurrence for the outer electrons of atoms. 

What we suggest in these cases is that the one-electron wave functions 
should not be assumed to be of the form (Hartree 1946) 


b(r)=R(n, l, r)P, imi (cos 6) exp (1m) xi(s) 
but rather of the form 


£(r)=R(n, l,m, 0, r)P,," (cos 8) exp (im,¢)x,(s). 


tl ll I ae, ag 
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Fig. 1 
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(b) Z=1-07. Energy at A is —0-5738, at (c) Z=1-09. Energy at A is —0-6045. 
B —0-5759 and at C —0-5735. 


i d A’s in atomic units. 
Energy contours for a two-electron atom. Energies an 
purine are shown by 4’s and saddle points always have a ee 
crossing over at them. The contours are symmetrical about A,= 
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Condon and Shortley (1951) go some way towards this in working out 
Hartree’s equations (eqn. 814) for the case where the radial function is 
assumed to depend on ,, J;, and m; (cf. eqn. 813 where m,; has been 
mistakenly omitted). 

All our remarks apply with equal force to Fock’s equations, because 
Fock’s and Hartree’s equations are the same for a two-electron system 
when the electrons have opposite spin, and it is on such a system that 
we base our remarks, (The reason for this equivalence is that the © 
exchange operator only operates on electrons of parallel spin.) 

There does not seem to be nearly as much opportunity for this non- 
uniqueness in the case of Hartree’s symmetrized equations (in which an 
electron moves in its own field): since with these equations the one- 
electron wave functions must be independent of the spin, and if there is 
spherical symmetry the radial parts of the wave functions can only 
depend on the quantum numbers 7 and /, not on m or o, 


§3. SmLF-TRAPPED ELECTRONS IN A METAL 


In §2 we came to the conclusion that incomplete screening and polar- — 
ization, if present in too great quantities, were liable to make the solutions 
of Hartree’s or Fock’s equations non-unique, and, of course, the usual 
wave function, which is always the simplest, would very likely not be 
the best to take account of these effects. Polarization must be very 
large in a solid, for in this case the electrons must be able to be found 
with equal probability throughout the solid, but at the same time they 
must avoid each other as well as possible. The usual wave functions for 
a metal are Bloch’s, in which the electrons cannot avoid one another at ~ 
all. This then is an ideal example to test whether Hartree’s equations 
have a unique solution. (We have not examined Fock’s equations.) 

To simplify matters we consider the Sommerfeld model of a metal in 
which the metal ions are replaced by a smoothed out positive charge of 
density |e¢|/Q2. The usual one-electron wave functions are just the 
Bloch waves exp (ik.r). We have now to decide in what way to 
generalize these wave functions. We shall suppose that only one or two 
electrons cannot travel through the length of the metal and we shall take 
the wave function of these electrons to be of the form exp (—Ar). That 
is to say we take 


exp (— A,r) exp (—Agr5) ITp,(r;) 
v 
for our complete wave functions, where #(r) extends throughout the — 


metal. Using the variation method we find that the 7, must satisfy 
the Hartree equations 


h? e2 dr’ 2 |, P+ ?+ 2 
aaah Im Vb + |- 3) | r—r’ te] Boks se eee eG 1 dr’ pb, 
) 


rh | 


=<y,  (i=34...), (1) 
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where ,=(A,3/7)1/? exp (—A,r) (‘=1 and 2). The energy is then given by 
; h2 ; 1 e2 
BE be eee x me z 
Qrdd=— a5 2 [PV dr—2[ [YP 5dr, dry 


lS, Maher 
een |; | |; | pean aa: 


where the %,%,... have been chosen to satisfy the eqn. (1). We have 
then to choose A, and A, so as to minimize H(\,, A.) ; it may, of course, 
turn out that A,, and A, are both zero ; this would mean that our suggested 
wave function is worse than the plane waves. If we write 2 |, Pp and 
call the potential for the free electrons occurring in (1) U , the energy 1s 
(using eqn. (1)) 


h2 
Bly, Agd=— 5,2 [HY de +Be,— | pU dr— [| Gat+Hete) 
1 eé? 1 e2 
Batts Ura 5 || 2%? +202 Het ppl = dry dr 
(2) 
In other words we do not need to know the ,(i=34...), but only 
2 Eis Ps and U. All these can be obtained approximately by the Thomas— 


x 


Fermi method 
87 
P— 3/3 [2m(#— CB's 


and. 
Ee= [[(H—U)p+Up]ar-+H2,1,0),  » . + (3) 


where # is chosen so that 


Jo dr= [pare 1, 0). 


The expression (2, 1,0) means either 2, 1, or 0, depending on whether 
there are 2, 1, or 0, bound states, i.e. A; and A,~0, A, or Az=0, or A, 
and A,=0. The reason for adding the term E(2, 1, 0) to the sum of the 
eigenvalues is that the Thomas—Fermi method sums over all the eigen- 
values in the potential U including the Is. (We are supposing that 
boundary conditions are applied on the surface of a sphere centred on 
the trapped electron, so that the quantum numbers of all the electrons 
bound or free are the usual n, J, and m.) However, if there are two 
bound ls electrons, we must not fill the Is free state with electrons in 
order to satisfy the Pauli principle. We must therefore add to the 
Thomas—Fermi expression for the sum of the eigenvalues the energy 
necessary to raise 2, 1, or 0, electrons to the top of the band. The 
corresponding change in p is, of course, insignificant. It turns out that Z, 
the energy of the top of the band before the 1s electrons are put up there, 
is independent of the number of bound electrons, so this possible variation 
does not occur in the equations for p and 2«;. 
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In the usual way we take U=¢/r (definition of £) as the dependent 
variable in the Thomas—Fermi equation for p and use the fact that 


1 
VU =" /r=—4re? | — D totuttee | ewer 
We then obtain the following differential equation for ¢ 
LS D 87 om)8/2, B sa 
TELE ES +A, exp (—2Azr)+A, exp (—2Agr)=— a7, (2m) [o—¢/r Pes Q 


(5) 
where A,=A,3/7 is the normalization factor. One boundary condition 
is that > 0 as r—> ©, this is only possible if H is given by 

3 1 

aa (2mi)*2— 5) er 
irrespectively of whether A, and A, are zero or not, as we stated above. 

If we did not have to take (2, 1, 0) electrons from the bottom of the 

band and put them at the top we could state the conclusion that the 
top of the band is unaltered by self-trapping, in agreement with Friedel’s 
(1952 a) more general conclusion. When, however, we do have to raise 
electrons to the top of the band, it is clear that the whole band will be 
shifted upwards and narrowed. We consider this effect in greater detail 
in § 4. 

The Thomas—Fermi eqn. (5) that we have obtained would have to be 
solved numerically ; to avoid this we use the usual approximation 
(Mott 1936) of linearizing, that is to say we assume ¢/r<H#. This will 
be true, as we shall show, for A, and A, sufficiently small. Equation (5) 
now becomes 


C*— «l= —4re*[A, exp (—2A,r)+A, exp (—2A,r)],  . . (7) 
where 
: 327? 
3H 
Using the boundary conditions that £0 as r->0 and roo we find 
the solution 


(2mys25 eS ee 


ty] 4,4, rane 
an ene een (—«r)—exp(—2A,r)] (eae +r exp (—2),r) EL nat 


cial, eee 


If 2A,<k; which is always the case in our problem, ¢/r is a decreasing 
function of r with a maximum at the origin, where 
C 4¢?),8 4e?),8 


ies Creasy Gat Vee pare FS (10) 


roo (k++ 2A,)? (k-+2A,)? oe eae 


so ris for A; small enough this can be made as small as we like compared 
with FH, 


Self Trapped Electrons in a Fermi Gas 1263 


Now using the known relation between p and ¢” (eqn. (4)), but not using 
the approximate eqn. (7) for ¢ we find that eqn. (2) can be written 


Oe 2H 
By, Xe) —B(0, 0)=— 5-2 | she de+ (2 1,0)—5 | f Waa bah) 


e2 1 127 1 
x Fade, deat 53 ft" dr [rt (vittye— a ) ee Lt) 


The straightforward procedure would be to calculate the integrals 
involving ¢ from the approximate ¢ (eqn. (9)). This gives 
ne 12e? 


5e? 
AA, A,)—H(0, 0)= Im CES) eo 16 (Ay+Ag)+ R RA, A) 
Se? 
Er heen so. > (12) 
‘where 
24 dic(A,+A,)t4A 2+ 12A,A,+4A2 
yy ay 3k Pela, i ny j 
FUN Na) 2AM Se ORO 
and 


Re ght 6c(Ay Ay) +4AP + 12A,A,+- 40/7 
S(A,, A eset d; (-F 2A,)8(e + 2A, )8(A, Ay)? 

Another possible procedure would be to say that the second derivative 
of an approximately known function ¢ is likely to be very poor, and that 
therefore we should get rid of the ¢” in (11) by using the exact differential 
eqn. (5). In fact, one can show by using an extension of Milne’s (1927) 


integrations by parts and eqn. (5) that 
ik ah eee BI ee 2 er 
cor dr | reat toed) dr— 5 | rh (st de— g) ar 


in 
42/05 [ "(trie a) |e 


using this in (11) we obtain from the approximate C (9) 
h? 5e? 
E(A,, A,)—E£(0, 0)— mM (Ay?+A?) — ile (Ay+Ag) 
24e? 4e? 
4S Ry Na) FT Ar Aa) + (13) 
where R(A,, Az) is the same as in eqn. (12) and 7'(A,, Ag) is 
ZA PAs 
B18 + 61c2(2A,-+ BA,) + 44e(4,2-+ 16A,A;-+ 92,2) + 8(2A2-+8A2Aj+ 12A,AP+3A/) 
(FB, Pe FEA POE AG) 
The most striking conclusion that can be immediately drawn from 
eqn. (12) or (13) is that the energy, E(A,, Ag) —E(0, 0) is not stationary 


at A,;=A,=0 due to the term —5e2(A,-++A,)/16, which is half the self- 
energy of the two bound electrons. That is to say the plane wave solution 
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of Hartree’s equations does not give even a stationary—let alone a 
minimum—value for the energy with respect to variations in the one- 
electron wave functions. At first sight this is very unexpected since 
Hartree’s equations are deduced on the basis that variations in the wave 
functions produce no first order energy change. The explanation is 
that we have varied the plane waves by more than the first order in going 
from exp (ik . r) (or at least the corresponding 1s wave function containing 
the Bessel function J,/,) to exp (—Ar) no matter how small A is. This 
shows that great caution must be exercized in using the variation principle 
when the wave functions have almost constant density over the whole 
system, cf. Coulson and Fischer’s (1949) remarks on the molecular orbital 
method. 


Fig. 2 
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The reciprocal of the radius of the orbit A, and the energy saved (ev) when 
one (fig. 2 (a)) or two (fig. 2 (b)) electrons are trapped at a definite site. 
These are plotted as functions of the band width H (ev). The upper 
of each pair of curves corresponds to eqn. (12) and the lower to eqn. (13). 


Consider now the variation of E(A,, A;)—E(0, 0) along the lines A,=0 
Le. only one electron trapped, and \,=),, ie. two electrons of opposite 
spin trapped in the same wave function. In either case the energy has a 
minimum (notice that the discontinuities of the form (2, 1, 0) have dropped 
out of the expressions (12) and (13) for the energy), which is most easily 
found numerically by regarding «/A as the independent variable, when 
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0H/0A=0 becomes a simple equation for A as a function of «/A. The 
minimum is shown as a function of «, and H the band width (eqn. (8)), 
in figs. 2 (a) and 2 (b) for A,=0 and A,=A, respectively. The pair of 
curves in each case refer to the alternate formulae (12) and (13); it is 
seen that which formula is used makes little difference and this leads us 
to believe that our linearization eqn. (7) of the Thomas—Fermi eqn. (5) 
is quite accurate despite the fact that for band widths less than about 
0:7 ev C/r at r=0 as given by eqn. (10) is bigger than H. In the same 
figures we show the energy saved, nearly 1 ev per electron trapped. 

Now it may be asked : is the minimum along the line 4,;=A, a minimum 
for variations off this line (cf. §2)? The answer is easily seen to be yes, 
for the right-hand side of 

OPH (Ay, Ag) |? CH PH 
[ OA0Az | < OA? OAs? 
contains the dominant term (f?/m)?._ We can therefore say that if we are 
only interested in trapping electrons at one site we save most energy by 
trapping two in the wave function exp [—A(r,+7,)], when only wave 
functions of the form exp (—),7;) exp (—A9r2) are allowed. 


$4, CORRELATION ENERGY 


So far we have only considered trapping one or two electrons at a 
definite position in the metal. In order to obtain an appreciable save 
of energy compared with the Fermi energy for plane waves we must 
trap a non-zero fraction of the electrons at different positions ; for 
example, we might trap one electron on each lattice point of a body- 
centred lattice. It should be remembered that we are only considering 
a Sommerfeld model so this lattice has nothing to do with the position 
of the ions. For the relations between these lattices see §5. However, 
if the trapped electrons are too close together compared with 1/A we will 
loose energy due to the overlap. Let us puta Wigner—Seitz sphere around 
each lattice point, of which the radius is say, a/A. Then the amplitude 
of the wave function of a trapped electron on the surface of the sphere is 
exp (—a) and with a equal to 2 or 3 this is 0-14 or 0:05 ; ‘we may hope 
small enough. The corresponding fraction of electrons trapped is 
f= 23A3/4na° (if we trap one per site) since Q is the valume occupied 
by a charge e. So we have the following contribution to the energy 
saved [H(A, 0)—(0, 0)]f:N, where NV is the total number of electrons, 
but this is not the only term. It will be remembered from $3 that for 
every electron we trapped we had to raise an electron from the bottom 
of the band and put it at the top because of the Pauli principle. The 
Pauli principle was there used in the simple form that there must not be 
two electrons with the same set of four quantum numbers. This is no 
longer of any use to us because the quantum numbers of the trapped 
electrons are ”, J, m, and o, while the quantum numbers of the free 
electrons, which are now Bloch waves, are hz, hy, &,, and ~o. We shall 


1266 J.S. Plaskett on 


therefore suppose that /,N electrons—the number that are trapped— 
have to be removed from wave functions at the bottom of the band and 
placed in wave functions at the top, remembering that as this is done 
the energies of the top and bottom of the band are shifting. We shall 
suppose that the energy of a wave function with quantum numbers k 
is h7k?/2m. The original band width is # and so this transfer of electrons 
requires an energy 


3h? (3n?\ 2/3 
sin(a) ¥ [athena 
where we have supposed that all the trapped electrons have the same spin, 
so that the /,N electrons of lowest energy have the opposite spin. We 
must subtract from this 
h? , _(30?\2/8 
[NE= x fiN (=) 

(from eqn. (6)) the energy required to raise f,N electrons to the top of 
a band that stays the same width # ; this term has already been added 
in § 3, it was, in fact, the (2H, E, 0) of eqn. (3). 

If, on the other hand, we trap two electrons at each lattice site so that 


the fraction of electrons trapped is / 2—= 2.23A3/47a8 then we find for the 
total energy saved 


62= 58a) | fyb fe |+| 20, 9B, 0) | fay. 


(14) 
This formula is not quite correct because the \ has been worked out for 
a certain band width EH, whereas, in fact, the band width is less; in 
addition the free electrons that are doing the screening have a longer 
average wave length and this will also effect \. We suppose these 
effects small as they certainly will be for big #, and plot the expression (14) 
divided by N in fig. 3, which therefore represents approximately the energy 
saved per electron when a certain fraction of electrons is trapped (it has 
ee plotted for a=2, and 3). The corresponding energy saved, &,, 
2 only one electron is trapped at each lattice site is also shown in 
Coes 
hice the same figure the usual expression for this saving of energy on 
: e ele energy is shown; it consists of the exchange energy—obtained 
, y ; tae Sane of parallel spin apart—and Wigner’s expression 
r the “ correlation energy ’"—due to keepin i 
SS ping apart electrons of opposite 
7 0-46¢e2 0-29e? 
ie r+ 5: lag 
where 2=4rr3/3 and ay=h*/me?. 
less than this though it is compar 


presumably our expression could be 
exchange. 


g! 


(15) 


Our saving of energy is considerably 
able with Wigner’s correlation term, 
increased considerably by introducing 
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Note that our curves in fig. 3 stop for sufficiently narrow band widths : 
at these points all the electrons have become trapped (/', or fo=1) and the 
Sommerfeld model (i.e. uniform positive charge) of a metal becomes 
insulating in this approximation. This situation has been envisaged by 
Wigner (1938) who calls such a system an ‘ inverted alkali metal ’ because 
“the positive charges are uniformly distributed while the negative electrons 
are localized because of their repulsive action”. Wigner has worked out 
the energy of such a set of localized electrons. He assumes that any 
given electron moves in the field (Wigner 1938, eqn. 10) 


— 3e?/27r,--r2e?/2r,3. ee we lO) 


Fig. 3 


OL OS © Cy WL Tex) K au. 
Cyl OFF Kes 2Oy stay ek: E ev 
22 14 ie) 7 54 4.3 fr, au. 


Energy saved in atomic units when a fraction of the electrons are trapped on 
a lattice, labelled &, or & depending on whether one or two electrons 
are trapped per lattice site. a—3 and a=4 show how much overlap 
has been allowed between trapped electrons. &’ is the ordinary exchange 
and Wigner correlation terms added together (eqn. (15)), ¢@” is the 
energy saved by localizing all the electrons (eqn. (17)). 


For r<r, this is the field due to a uniform positive charge distribution 
occupying the volume 47/3r2; for r>r, the potential spon be 
—e?/r, but since the wave function is very nearly zero for these r’s Wigner 
uses the potential (16) for all r. The wave function for a trapped electron 
is p,ocexp (—ar?/2) where w=r, 3?aq_71/? and ag=h?/me?, and the eigen- 
value of this wave function is (Wigner 1938, eqn. 11) 


€,= —3e?/2r,+3e%ag)/?/2r,3”. 


The total energy of this system of electrons is 


1 AE 
U PY per, 
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whereas the energy of the free electron gas is 


h2 ] 1 é 
Sin = 25 BB 5 | | Reidy 


Subtracting these to give the saving of energy when all the electrons are 
trapped we find 


Be2 SB eFayzl/2 «= #2 2-21 Se? 
OC" =C ings —C con |-s. +3 per ~ Im re re e (17) 


Wigner (1938) takes into account the exchange energy of the free electrons, 
namely —0:42e2/r,. We must not do this since we wish to compare 
this calculation with our previous one in which only a fraction, f;, of the 
electrons are localized and where we used Hartree’s equations, not Fock’s. 

The saving of energy @j,,—E€ on given by eqn. (17)| is plotted in fig. 3 
(we have not made similar calculations for the case of two electrons 
trapped at each lattice site) as a function of the band width # for values 
of r, greater than 6ag; Wigner does not believe the ‘approximation is 
good below this. Note that the wave function has the value 


exp [— V/(7,/44z)], 
on the Wigner-Seitz sphere, which is 0-29 for r,=6ay. [This is bigger 
than our 0-14 or 0:05, but Wigner’s wave function varies like exp (—7r?) 
whereas ours varies only like exp (—r). 

Incidentally, the numerical comparison that Wigner makes between 
the correlation energy as obtained in his earlier paper (1934) (our eqn. (15)) 
and in the later paper (1938) (our eqn. (17)) where it is assumed that the 
electrons form a lattice appears to be wrong. He gives for the correlation 
energy using a lattice (footnote 25) 0-15/r, and 0-21/r, Rydbergs for r,=6 
and 8 atomic units, whereas it should be 0-03/r, and 0-10/r,. This 
considerably reduces the agreement with the 1934 paper, which gives for 
the same values of 7,, 0-32/r, and 0-36/r,. 


§5. NIcKEL OxIDE 


We have seen in § 4 that if the electron density is below some critical 
value—corresponding to r, between 6 and 15 atomic units—the Sommer- 
feld model gives an insulator. Before this result can be applied to any 
interesting case, say nickel oxide, we must, however, introduce the effect 
of the periodic potential. It should already be clear that our results in 
§4 for the Sommerfeld model are not very reliable because of the ad hoc 
method we used to introduce the Pauli principle. Now that we introduce 
a periodic potential we cannot expect to get good quantitative results 
but only an idea of what a crystal lattice may do to self-trapped electrons. 


The Hartree equations for the valence electrons in a periodic potential 
V(r) are 


h? » |b; [? 
— aa bet] Vinee fe E ae mes oe VeN 
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A possible solution of these equations is Bloch’s in which | #, |? is periodic 
and so the total potential for the ith electron is periodic. Let us imagine, 
however, that a solution exists in which two wave functions ¥%, and so, 
say, are localized, then the total potential is no longer periodic. However, 
we may write the potential in the form 
e277 dr’ 2” |, (?—-1/2 _, 
(oleae tert 
the first part of which is periodic and gives a certain band structure— 
nearly the same as that obtained in a self-consistent Bloch solution of 
eqn. (18)—and the second term is (we hope) a slowly varying perturbation. 
We now use the Wannier—Slater approximation (Slater 1949) in which 
the eigenvalues of our problem are given by the eigenvalues of the equation 
d” |b; ?—1/2 
[ ai) -ee] A ae | pee 

where E(p) is the energy of a wave function in the periodic lattice with 
momentum p=/k and p is replaced by the operator —ihV. We shall 
furthermore assume that the potential due to 2’|%, |? is about the same 
as that produced by 2"|¢,|?. If now we suppose that the energy of a 
Bloch wave in our crystal is simply 1?k?/2m*=p?/2m*=—h?V?/2m* we 
have exactly reproduced eqn. (1) with m* for m. In other words, all 
our graphs, and in particular figs. 3 (a) and 3 (6), remain the same and 
only the units of energy and distance are changed. Our new ‘ atomic 
unit’ of length is f2/m*e2=m/m* times the usual unit (0-5 x 10-8 cm) 
and our new ‘ atomic unit’ of energy is m*e*/i2—=m*/m times the usual 
unit (27 ev). 

If in NiO we say that r,=2-5/?/me? (this is on the assumption that 
we need only trap two electrons per nickel atom) then we have to make 
m*>6m in order to have an insulator. Notice the peculiar result of the 
Wannier—Slater approximation : the trapped electrons have no preference 
in where they are localized ; over the nickel ion, the oxygen ion, or in 
between, all are equally good in this approximation. In actual fact 
there will be definite preferred positions for the trapped electrons, and 
it seems likely that this will immobilize them. 

A similar problem is the variation of the activation energy of semi- 
conductors with impurity concentration, which suggests that when the 
impurities are closer together than 100 x 10-8 cm they give a conduction 
band (see Mott 1951). If we replace the perfect crystal by a uniform 
dielectric, constant K=13 for silicon, then our previous results apply to 
the remaining ‘ crystal of impurities’. Wherever there is an ¢? in our 
previous expressions we must now put ¢?/K and this makes our system 
insulating for r, greater than some multiple of /2K/m*e?—=Km/m* 
0-5xX10-8 cm, ie. 7,>39 or 97m/m*10°* cm for the usual choices of 
multiple; 6 and 15. This would be in quite good agreement with 
experiment if m* were equal to m, but for such an extended crystal the 
band width is bound to be very small so making m*/m very big. Wemust 
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agree then with Castillan and Seitz (1951) and Mott eae that this 
explanation, as it stands, is not adequate. Of course, if m* is really 
about the size of m then, accepting this fact, there is no difficulty in the 
interaction of impurities up to 10~¢ cm apart, but we have merely trans- 
ferred the difficulty to understanding the low effective mass. 


§6. Frrepen’s THEORY OF X-Ray SPECTRA 


When an x-ray photon is absorbed by a metal an electron is knocked 
out of a low lying level in one of the atoms into a state of one of the 
conduction bands. The usual theory of the fine structure on the high 
frequency side of the edge is that it is due to the variation of the density 
of states in the bands and the variation of the transition probability. 
Friedel (1952) has pointed out, however, that when an electron is 
excited it leaves behind a positive ion that must be screened. For 
example, in lithium absorption Friedel suggests that the screening is 
due to a bound 2s electron. In addition, he supposes that it is possible 
to screen by a 2p, 3p,... bound state. Each one of these screenings 
gives a different absorption edge, and Friedel’s explanation of the fine 
structure is that each edge produces a peak in the absorption. In a later 
paper (Friedel 1952 b) he goes further and supposes that the two possible 
spins of the bound electron, say, the 2s, will give two edges. 

From our examples of the non-uniqueness of Hartree’s equations, and 
especially from § 3 we believe that many, if not all of Friedel’s solutions 
of Hartree’s equations (or Fock’s) in which he traps either a 2s, or 2p, 
or... electron, are merely different approximations to one exact excited 
state, just as our solutions of Hartree’s equations in which we have all 
plane waves or one or two trapped 1s electrons are approximations to an 
exact ground state. We believe this will be the case if wave functions 
with all quantum numbers n, /, m, and o, are occupied up to some energy 
E, the top of the band (except, of course, for the electron that has been 
excited). This may not be possible ; for example, if we bind a 2p electron 
there may be no 2s wave function in the band to fill; in this event we 
clearly have a new excited state different from the one in which the 2s 
wave function is occupied. In other words, we say that approximate 
many electron wave functions are the same or different according as they 
have the same configuration or not. 

When we screen the positive ion by a 2p bound electron and the 2s 
free wave functions exist, we may still produce a new state—different 
from the one in which we have a 2s bound electron—by merely removing 
one of the 2s free electrons and putting it at the top of the band. How- 
ever, this does not resuscitate Friedel’s theory because the new state has 
an energy greater than that expected by an amount equal to the band 
width. 

Furthermore, if the system with a 2p bound electron really approximates 
to a different excited state then it must have been reached by a two (or 
more) electron transition, because we have had to remove one (or two) 
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electrons from the 2s wave function in addition to the K electron. We 
are not convinced by Friedel’s (1952 b) calculations of the high probability 
of such multiple electron transitions. True the transition probability 
may be high if the wave functions are taken to be plane waves, but it 
seems likely that if correlation is taken into account so that the electrons 
avoid each other, the probability might be much reduced. 
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SUMMARY 


The analysis of Goldstein and Rosenhead for the boundary layer growth 
on an impulsively started body has been extended to three dimensional 
bodies. The results are used to test various critera of three-dimensional 
separation. : 


$1. INTRODUCTION 


GotpstEIN and RosENHEADt (1936) have studied the growth of the 
boundary layer on an impulsively started two-dimensional body. The 
initial motion may be assumed irrotational and without circulation. 
Thus initially the potential flow about the body will give a very close 
approximation to the flow outside the boundary layer. These_conditions 
and assumptions cease to be valid as soon as the boundary layer separates 
from the body to form a wake. The analysis therefore gives an account 
of the initial motion of the boundary layer and the time and place of 
initial separation. In this note the method is extended to three-dimen- 
sional bodies ; again the initial external flow may be assumed to be the 
potential flow about the body. The iterative form of solution employed 
in paper I may be extended without difficulty to include all three-dimen- 
sional effects ; the successive approximations however increase rapidly in 
complexity. The complicated form for the skin-friction makes simple 
deductions of time and place of initial separation difficult ; these may 
however be found easily, in any given case, by the data provided in this 
note. 

This extension was carried out as part of work to investigate three- 
dimensional effects in boundary layers. Although valid for small time it 
seems probable that it will be the only solution in which it is possible to 
find curvature and cross-flow effects explicitly. Results of the solution 


are used to consider the conditions of three-dimensional boundary layer 
separation. 


ee 


* Communicated by Professor L. Howarth, F.R.S. 
{ Referred to as paper I. 
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§2. Bounpary LayvEeR EQuaTIons 


The three-dimensional boundary layer equations in their most general 
form may be written 


Tt Te set Te apt ge pKa Kyun 2 pS, (2) 
nett yt pe Kin K y= erase: (3) 
where eee ee ok = and K,=— a ot Sree (4) 


are the geodesic curvatures of the lines x=constant and y=constant, 
and U,V are the mainstream values. It is convenient to use the following 
abbreviated notation :— 


OG) wa MEO a ors 
ee Ps E Gop eb - + (5) 
and similarly 
cree OL) ares rales) =(5) 
male ge) a0" alan [6 Be ean 
ee ee eek ye N= KU KV —U’—V,.. :- (7) 


since these latter combinations occur frequently. 


First Approximation 
Initially the boundary layer is of zero thickness and by a simple exten- 
sion of Goldstein and Rosenhead’s arguments the equations of the first 
approximation are found to be 
Ou, Ob, CO 070, 


eid Fs ae oS oan abs (8) 


The boundary conditions on w and v are u=v=0 at z=0, u/U, v/V tend 
to 1 as z—> 0, u=U, v=V at t=0. 
The solutions of eqn. (8) satisfying these conditions are 


Q(t 2, 
i=) Orla) a= ort 7 ; erfn= ap aha . (9) 


z 


where an 
Q= 2(vt) 1/2 : 


The first approximation to w, may be obtained from eqn. (3) and is 
1 2 
w= 2HAN fy erfn— a8 (l—e" } ee on (LO) 


As noted by Goldstein and Rosenhead this value of w is not completely 
satisfactory, since it behaves like 2(vt)”7.Nn as 7— ©, however they 
point out that u/U tends to 1 for moderate values of 7 so that within the 


boundary layer w,=NO(ut). 
402 
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Second Approximation 
The equations for the second approximation are found in the usual 
manner by substituting w,, v1, w, in the ignored terms of eqns. (1) and (2). 
In this way the equation for vw, becomes 


id ae py aon 
bata Ot Sama a ayuae siya Ox The oy 1 1B, 
— Ku, + Kyw,?. . 2 Vela Ga 


When the values of w,, v,, w, are substituted into the r.h.s. of (11) and | 
use made of the irrotational condition of the external flow the equation 
takes the form 


Oru, Ou 2n Pe tr ie 
y Se — Ft = UN fen?) —1— pert ne apallr =e: =} 
+-(VV'+UM—UV)(erf*n—1). 
Putting 


tg=t{— UN for(n) +(VV' + UM—UD) foa(n)} 
the following equations for /.,, fz. are obtained : 
2 ee : 2 
for” +20 fer’ —4foi1=4 {ert fie Nee mile erfye—" + “a je on.) i} ,» (13) 


foo’ +2n foo’ —4foo=4 {erf? n—1}, Mera Se 
where dashes now represent differentiation with respect to 7. Equation 
(13) has been solved in paper I. This solution together with that of 
eqn. (14) may be written in the form 
Soi=Seip t%i(foier +Pilfeider 


where /,;, is a particular solution, (f2;).1, (f2;)¢3 are complementary 
functions and «;, 8; constants which must be chosen to satisfy the boundary 
conditions f,;=0,7=0; f,;>0,7— 0. It may be shown that 


3 ; 4 eo ee 
(farlo=2(2y?—1) ef? n+ aye” erfy+1—se-*+—e-*", (18) 


4 A 2 A 
(fo2)p=2n" ert? n+ ane” erty ea, oe te eee ee 


(fei)er=2n*+ 1, re ee 
2 : a 
(fei)ea= =a (29?-+ 1) erf n+-ne™, . Saas ne 
and. that 

7=— jigs etl 1° 1 4 

1 37) 2122, Bi= =n ( 1+ =) =0-8036, 

Hee et: ¢ 4 i 
Ca hae) eee O S00, B= =0-7183. 


It is easily seen that the derivatives needed to find the skin-friction, are 


3) 2 4 of, 8 
mene ar —_— I+ —- = . J 22 — — 
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In a similar manner it may be shown that 


Ve=t{—VNfy,+(UU+VM—VU')foo}. =... (21) 
Equation (3) then gives 
W = 2(vi*)' 4 | far(s 6) ds+-BY fox) yas}, ee) 
where A=—N?—UU"—VV+UM'+VM—U(V)—V(U') 
and B=M?—M(U'+V)—VV"—UU—VM—UM' +07 
PVC kV EK UG V2—U220'7 
RHEE IS NAN 


The integrals in eqn. (22) are given by 
Prix {faa for(s) ds= §(24°— 3) erf?y— aa (4—(4 (4y?—11) e~"} erf 
, 4(2)12 
a 10 + gore ert V/(2)n + a {39° +77} 
priya 
+P (2y3+8n) erf n+4(y?2+ ner—sh, ss Be (ee 
and . 


n 93/2 
Fy = [ fast ds = 3° erf? y+-y+ s—apert (2)¥?y 
2 2 —1° 2 21? 23 
+ 5-773 (29°—1) erf ye" + g— ye" 1 a9 (59°49) 


qri/2 2 
+8a{7 Gr? +30) ert sor Ne *—¥} eee (35) 


Third Approximation 


The equation for w, is found in a similar way to that foru,. The right 
hand side may be most conveniently collected together as follows : 


e? 0 
ea | (Che fe) erfn+ U(AP,, + BP 22) ae 


“O22 at 
+ {nerf salle") fa, publ, (2s) 
where 


C=—2N{K,V*—K,UV+0V+ UU}-+U4U"—M'+(VY} 
+UV{V—-M+(U)}, 
D={(UU+VM—VU"}{2K,V—K,U+U}+UV{(V"— } 
{VV +UM—UV}{U'—KV}+UV 2+ VVV'+M{UU'+ VU} 
+U2M’+UVM—UVU'—VOV—UXV) + VV’), 
=—UN?, 
F=N{VV'+UM-UV}. 
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This grouping of similar terms in 7 rather than of terms similar in the 
external flow appears to be the best method, since this latter grouping 
produces 12-14 equations for solution. 

A solution of (25) is obtained by writing 


Us=t? {Cfe1t+Dfsot+ Efsst Fist UAfss+UBfsg}, - - (26) 
where f3, is a solution of the equation 
fi’ + 2nfsi —8fai=4rv - +. See 


with boundary conditions 
fsi=9, ii fai > 9, ae, 


and 
r= fo, erty, ro= fos ert n, af 
1 mC l _nty \ G22 
LS (» erfy— =a (l—e™ I) on (» orl) = 1a ae ) On’ 
2 3 _n? 
Cio Foi; Me Se Fy 7 Boo. 


ie) ” 2:0 


In paper I two equations of type (27) have been solved, these two solutions 
denoted in that paper by f; and F',—f, are equal to fsi—f35 and f5,—fo9g 
respectively in the present notation. Thus once f'3;, f39, fs a f3¢ have been 
obtained, f,5, fz, can be found from paper I. 

As in paper I solutions of (27) are sought in the form 


f3i=X (7) erf? n+ Y (7) erf? 7+Z,(7) erfn+T'(n)+ 4 m+ Viti (28) 


where m and n are the complementary functions 


1 2 
m= 4n4-+- 122+ 3, n=m erf y+ ar (473+ 10) e-". (29) 
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The solutions are shown graphically in fig. 2 and tabulated in table 1 ; 
for reference their analytic forms are expressed in the Appendix. 

The derivatives needed to find the skin friction can be obtained by 
differentiating the expression (28) and evaluating at7=0. The derivative 
is then 


9 16 
Fa) = sat (MotaaWe © + + + GO) 
n=0 


N.B.—f2, is taken from Goldstein and Rosenhead’s paper and is given here 
for easy reference. 
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Substitution of the analytic forms now gives 
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256 
+ 5-75 =0-0286, (34) 
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——0-0679. (36) 
In a recent paper on a related problem, boundary layer growth on a 
spinning projectile, Illingworth (1954) has solved some other equations of 
type (27) by the method of variation of parameters. The resulting 
integral form has been found quite difficult to evaluate explicitly but 
gives a very simple form for the wall derivatives. This method applied to 
our equations leads to the result 
) 
On Jo 
These integrals have been evaluated numerically and used to check the 
values given by (31)-(36). 
The third approximation for v, may be found in exactly the same way, 
the solution is 


Ug=E{O* fg, +D*fg.+L "fog tPF *fgat VAf35+ V Big}: eed OK. 


= 3\. (n—m)e" r,(s) ds. 
0 


where 

O*=—2N{K,U2?—K,VU+VV+UV3+UV{U"—M'—(V)'} 
+V+{V—M-+(U)'}, 

D*={VV'+UM—UV} (2K U—K,V+V3+UV{0—U"} 
+{UU+VM—VU'}{V—K,U}+V02400U' 
+M{VV+UV}+4+ 02M + UVM —VVU'—UV'U" 
—V{U')+U(U)’, 

E*=—VN?, 

F*=N{UU+VM—VU’}. 
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The total results are now summarized in the skin-friction components 
which are given as 
Ou nes 
Fe), = HOUT 12840441 60730 +1-4370(7 V+ UM—UYV)} 
+#{—0-1501C—0-1127D+0-0306#-+0-0286F 
—0-0827U A —0-0687UB}}], ey (38) 


ov =F 
) =4(vt)-¥/2[1-1284V +4{—1-6073VN+1-4370(UU + VM—VU’')} 


az], 
+#2{—0-1501C*—0-1127D*-+ 0:0306H*+0-0286F* 
—0-0827V A —0-0687V B}]. Pee er Pe Ge 39) 


§3. APPLICATION TO AN ELLIPSOID WITH A DISCUSSION OF THREE- 
DIMENSIONAL SEPARATION 


The results of the last section have been used to study the boundary 
layer on an ellipsoid of axes in the ratio 30:6:1. The ellipsoid was 
considered moving parallel to its middle axis so that the external potential 
flow was known and the various functions needed to find the coefficients of 
the skin-friction in eqns. (38) and (39) had been tabulated in an earlier 
problem. 


Fig 3 


Before giving the results for this body, some possible definitions of 
separation in three-dimensional flows will be mentioned. These include 
separation as defined by the disappearance of the skin-friction in the 
external streamline direction on a general surface and by the disappearance 
of the chord-wise skin-friction component on a wing. Hayes (1951) 
points out that these criteria are not always consistent and suggests 
alternatively that one must expect regions of separated and unseparated 
flow to be divided by a skin-friction line.t 

+ A skin-friction line is defined as a line whose tangent is in the direction 


of the resultant skin-frictional force and as such it may be regarded as the 
limiting streamline direction as the surface is approached. 
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The definitions and their effects are now considered for the infinite 
wing of fig. 3. In this figure let x be measured normal, and y parallel to 
the leading edge and let s be in the external stream direction. Also let 
Ty: Ty: 7, be the skin-friction components in these directions. Then the 
line of separation might be taken to be AB or A’B’ depending on whether 
+,=0 or 7,=0 is chosen as the separation condition. Now consider 
Hayes’s criterion applied to the same problem. Since the wing is infinite 
this criterion implies that separation is defined by a skin-friction line 
parallel to the leading edge, that is by a line on which there is no chord- 
wise component of skin-friction. By definition AB is the only line with 
this property and so in this case Hayes’s criterion gives two-dimensional 
separation. 

Returning to the ellipsoid it will be seen that there are two lines where 
the external and surface streamlines coincide ; these are the lines of sym- 
metry joining the front and rear stagnation points, ACB and ADB in 


Fig. 4 


D 


S 


A 


Incident Stream 


fig. 4. On these lines the first two definitions of separation lead to the 
same result ; the point of separation being the point where the skin- 
friction component along the lines disappears. Also on this body Hayes’s 
surface streamline would be expected to pass through these two points. 
This has been tested by the numerical method described below but the 
method is not accurate enough to give a completely conclusive result. 
The most that can be said is that a surface streamline starting from the 
separation point on ACB reaches ADB in the region of separation. 

In the actual determination of the time and place of separation by 
eqns. (38) and (39) the second approximation was used first and the results 
checked by including the third approximation. This procedure showed 
that accurate results could be obtained for the inner half of the ellipsoid. 
On the outer region the influence of the third approximation was so large 
often altering the initial time of separation by over 30%, that the results 
in these regions can only give a guide to the actual boundary-layer 
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characteristics. Also potential flow is used throughout the investigation 
although it ceases to be accurate as soon as any part of the boundary layer 
separates. However, since the times used in this investigation are small 
the wake will be thin and so potential flow will be a reasonable approxi- 
mation to the actual flow. 

Subject to the restrictions of the last paragraph the growth of the 
boundary layer on the ellipsoid will now be described. The boundary 
layer first separated on the line ACB (fig. 4) at the rear stagnation point 
at time UT/C=0-0196} by the second approximation and UT/C=0-0181 
by the third approximation. After this time the separation point on 
ACB moved forwards and separation began on neighbouring streamlines 
at points not coinciding with the rear stagnation point. Streamlines 
further along the wing remained unseparated right up to the rear stag- 
nation point. At 7’U/C=0-04 separation takes place on ACB at an 
eccentric angle of 6° from the rear stagnation point and streamlines over 
the centre half of the body were separated. Separation was complete at 
TU/C=0-05 by which time separation took place on ADB on the trailing 
edge close to the wing tip. At this time the separation point on ACB 
has moved to 8° eccentric angle from the rear stagnation point. 

In addition to the test described earlier the surface streamline given by 
Hayes’s criterion was compared with the separation line defined by 
7,=0 when UT/C=0-03. The lines coincide on ACB by definition and a 
maximum divergence of 2° in eccentric angle was obtained. Hayes’s 
line always gave earlier separation. 

The results of this last comparison together with that of the yawed 
wing will now be used to state the following tentative separation criterion: 


The line of separation on a general surface is the last surface streamline 
on which the skin-friction in the external streamline direction does not 
become negative at any point. 


Surface streamlines are chosen since they give no component of skin- 
friction normal to the separation line. This implies that the velocity 
profiles normal to separation are of similar form to the two-dimensional 
separation profile. The actual surface streamline chosen is based on the 
two examples above which show that Hayes’s separation line 7, is never 
negative. 

ACKNOWLEDGMENTS 

T would like to thank Professor L. Howarth for helpful discussion on 
this problem, also the Department of Scientific and Industrial Research 
for a maintenance grant. 

REFERENCES 
GoipstTEry, S., and RosenuEaD, L., 1936, Proc. Camb. Phil. Soc., 32. 


Haves, W. D., 1951, NAVORD Report 1313 (Nots 384). 
Intrincworts, C. R., 1954, Phil. Mag., 45, 1. 


UT/C is a non-dimensional number based on time 7’ and a standard 


velocity (U) and length (C). 


1282 L. C. Squire on 


APPENDIX. 
fs1 
Xy=— 3 (49+ 249243), 
Yym matter" (n+ 0) + dia? Ont Bn) 


al 8 
Z=—7 10-2" (y?-+5)+ o-” ae a —1(4oyy?+20,+1), 
9 
yarn (Ey op) 4 stern (2 49h 


Il 4 
— 71 (352+ ie) = ale ew 


45 
“a ee 
+ (44+ 127?+-3) ext 4/3n— ert / 2? 
oy 4 
4,=— ie (1+ ea), ; 
8 32/2 9/3 1 
Visio aval 3646 Shes ee SE ee On 
S32 
1 
X y= ¥G (20y4— 1292-3), 


¥ y= go" 1? (1099 Ty) + BPM? (2984-31), 
Zg= +h ea} (5y2—4) + Jo" By(4n2-+5n) — F(4era9?-+ 241), 


1 2 2 4 2 
Pa 3g" 79? (549 +1739) + $8q0- 2" 7 V2 (2+ 1)— saan Very 


aie 
V/(3)7-1 erf +/3n[4n4+ 1272+ 3], 


TR 
8 2 


oa) 


1 3 
a. — 5 (ton! — 6r? -5), 
va oC 


BY ee 1 
=te m* 7-1/2 (1073— 7) + a7 ae ace qBont? (44+ 472+ 1), 


Z,= je" 7-1 (107?—8)+e"" (Fert hel (28+) )) + 5Pe—(1+ 429 


Boundary Layer Growth in Three Dimensions 1283 
: 3n* 3/2 3 ‘ Saeed 2 64 -1 3 
Ls Tee ln (5473+ 173n)+-e-2" a Bon?— 457 (11y+ 47?) 
4 —* 7-1/2 : —1/2 
= Gy Se ath (B2—329)+ is Cate Oh 


64 
aegral E 1/3 erf 4/3n— Te 4/2 erf vn | [47n4-++ 127?+ 3], 
16 


16 3 64 5 
= Se ae asa ae we 
a (+5v3 v2)" aa 


1 
X= — zy l4y+127?+3}, 


il . ] 

¥ p= — gre (6n8+-11y) + 53 Bat tt? (4qt-+129?+3), 
4 eur 2 —17? (9,8 

Zy=— 37 Me" (y+ 1) + 75 Boe (29° +5y), 


5 4 4+ 5a 4 
T = Sy tito (B78n a 9iln) ie" 7 UF) Ba —*n— 27° | 


180 : 6 3 
8 2 
= 4/(3)7—1(4y4-+ 1277+ 3) erf V3n— TB (2)7-1 e-” erf a/2n; 
16 
A, car a7 


[ 1284 ] 


CXLIII. Delayed Disintegration of a Nuclear Fragment 


By A. Sotuem and 8, O. SORENSEN 
Fysisk Institutt, University of Oslo, Norway * 


[Received August 20, 1954] 


A NuMBER of reports on experiments with photographic plates exposed to 
cosmic radiation have recently described the observation of disintegrations 
produced by heavy fragments at the end of their range (Danysz and 
Pniewski 1953, Tidman et al. 1953, Freier et al. 1954, Cheston and 
Primakoff 1953, Hill et al. 1954, Bonetti et al. 1954 and Powell 1954). Itis 
believed that the phenomenon is due to the existence within the fragment 
of a neutral hyperon A°, and that the disintegration is caused by a trans- 
formation of the A° particle. This transformation may be either a 
‘mesonic’ decay A° > P+7-+@ (37 Mev) in which a negative 7-meson | 
escapes from the disintegration, or it may be a ‘non-mesonic’ decay 
giving a neutron and ~178 mev energy release leading to an excitation and 
subsequent break up of the fragment. An event bearing on this new 
type of nuclear excitation has recently been observed in our laboratory. 
It was recorded in a stack of stripped Ilford G5 emulsions, each 600 
thick, exposed at high altitude by means of free balloons during the 
Sardinian Expedition 1953. A photo-micrograph of the event is shown in 
Plate 35. 

From a nuclear disintegration A of type 12+-4p in the nomenclature of 
Brown et al. (1949) is ejected a particle of range 243. An exact deter- 
mination of the charge of the particle is not possible on account of the 
shortness of the track. But a comparison of the shape of the track 
with those of other charged particles stopped in the emulsion indicates 
that the track is very similar to that of particles with charge 2e. In B, 
where the particle is brought to rest, or is moving with very low velocity, 
it disintegrates into two singly charged particles of approximately equal 
range ~73004. We regard the equality of the two ranges as accidental. 
Both particles end in the stack of stripped emulsions, and show the char- 
acteristic appearance of tracks due to protons. In the following 
discussion the particles C and D will be assumed to be protons, each of 
kinetic energy ~44 Mey, although we cannot yet exclude the possibility 
of one or both of the tracks being due to deuterons or tritons. More 
exact mass determinations of the tracks by means of scattering and gap 
Sy are in progress. The angle between the two particles C and 

a827 1°57, 
ee ee eee 
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Discussion 


From reasons thoroughly discussed in previous reports we exclude the 
following three possible interpretations : 

1. Chance coincidence of the end of the fragment’s path with a two- 
pronged star. 

2. Collision of the fragment with a nucleus in the emulsion. 

3. Decay of an excited nucleus. 

We are left with various alternative possibilities, and will first discuss 
the event considering it to be an example of a delayed disintegration of a 
helium fragment due to the presence within the nucleus of a A° particle. 
As neither of the tracks C and D is a 7-meson we consider a ‘ non-mesonic ’ 
decay of the A° particle with an energy release of ~178 Mev. The 
energy and momentum balance in the disintegration process may then 
be used to determine which isotope of helium is emitted from star A. 
A certain ambiguity exists in calculating the energy release due to the 
unknown binding energies of the A° particle in various nuclear species. In 
the following, however, we take for the A° particle the same binding energy 
in the parent nucleus as for the neutron in the corresponding normal 
nucleus. 

Alternative 1: Assuming a single neutron to compensate for the 
momentum of the tracks C and D in the disintegration, i.e., the fragment 
being a *He,* nucleus, we have the following energy release : 


Sum of the kinetic energies of the two protons C and D : 88 Mev 
Kinetic energy of the neutron : 116 Mev 
Binding energy of *He,* : 8 Mev 

Sum : 212 Mev 


This energy is too high by an amount of ~34 Mev to be provided for by a 
Q-value of ~178 Mev. Even a strong reduction in the assumed value of 
the binding energy of the *He,* nucleus cannot remove this discrepancy. 
We may therefore exclude the possibility of the fragment being a ?He,* 
nucleus. 

Alternative 2: Assuming two neutrons to compensate for the 
momentum in the disintegration, i.e., the fragment being a 4He,* nucleus, 
the minimum energy release in the process is obtained if the two neutrons 
recoil in the same direction with equal energy. We obtain the following 
energy release : 


Sum of the kinetic energies of the two protons C and D: 88 Mev 
Sum of the kinetic energies of the two neutrons : 58 Mev 
Binding energy of *He,* : 28 Mev 

Sum : 174 Mev 


This energy is sufficiently low to correspond to a @-value of -~178 Mev. 
The minimum energy release may be even lower since the estimate of the 
binding energy is almost certainly too high. Hill e¢ al. (1954) estimate the 
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binding energy of a A° particle bound in a 4He,* nucleus to be ~d4 Mev, a 
value which has to be compared with the energy ~20 Mev required to 
remove a neutron from an alpha particle. The energy release in the 
disintegration is therefore consistent with the fragment being an alpha 
particle in which one of the neutrons is substituted by a A° particle. 

The energy and momentum conservation laws do not exclude the frag- 
ment being a heavier isotope of helium. However, the instability of 
these nuclei in conjunction with the very low binding energy of the AS 
particles makes this very improbable. 

Due to the short length of the track, and consequently the uncertainty 
in the charge determination, we can not definitely exclude the fragment _ 
being a nucleus of charge 3e. This, however, would be very difficult to 
reconcile with the conservation of charge, as only two singly charged 
particles are ejected in the disintegration. 

In the previously mentioned reports the ‘ mesonic ion *-hypothesis has 
been discussed, considering the fragment as ejected from the parent 
disintegration with a 7-meson around it in a bound quantum orbit. The 
disintegration is then attributed to capture by the fragment of the 
7-meson with an energy release of Q~142 Mev. In the present case the 
‘mesonic ion’ would consist of a lithium or beryllium nucleus with a 
negative 7-meson rotating around it in a Bohr orbit. The evaluation of 
the energy release under this assumption, however, gives the energy of the 
disintegration to be too high for a z-meson star. 

Assuming the emission of neutral particles of sufficiently high energy 
and momentum, the disintegration could be attributed to the nuclear 
capture of a negatively charged hyperon. We can, however, definitely 
exclude the fragment being of unit charge, and regard this possibility as 
very remote. 
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A nuclear fragment, ejected from the disintegration at A, produces a track of 
length 243 » and comes to rest at B where it decays in two charged 
particles C and D. 
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CXLIV. A Cloud Chamber Study of Cosmic Ray Interactions Underground 


By H. J. J. Brappicx and B. Leontic* 
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ABSTRACT 


A multiplate cloud chamber has been used to study the interactions of 
cosmic rays underground. The nature of events previously reported as 
pairs of associated penetrating particles (A.P.P.) is discussed and an upper 
limiting cross section for the formation of such particles is estimated as 
2x 10-°° cm? per nucleon. An explanation of the high cross section found 
previously is put forward. A cross section of 3-2+1-8 x 10-%° em? per 
nucleon for the nuclear interaction of cosmic rays underground has been 
obtained. 


§ 1. INTRODUCTION 


THe production of penetrating secondary particles by cosmic rays 
underground was first reported by Braddick and Hensby (1939), who used 
a cloud chamber at a depth equivalent to 60 m of water. The phenomenon 
was subsequently investigated by Braddick et al. (1951), George and Trent 
(1951) and more recently by Lovati et al. (1953), Deutschmann (1953), 
George et al. (1953), Appapillai et al. and others, all using cloud chambers. 

In some of these experiments the cross section found for the production 
of single penetrating secondaries was an order of magnitude greater than 
that associated with the production of penetrating particles in stars 
produced by «-mesons (George 1952); it was therefore suggested that 
a separate process was involved. The aim of the present work was to 
re-examine this phenomenon as well as other interactions of cosmic rays 
underground. 


§ 2, EXPERIMENTAL ARRANGEMENT 


A cloud chamber containing six lead plates, each 2cm thick, was 
operated at a depth equivalent to 26m of water. It was situated in a 
tunnel 80 cm below the rock roof. The lead plates were 30 cm wide, while 
the illuminated depth was about 14cm. The chamber was of the con- 
ventional cylindrical type, 45 cm in diameter, illuminated by flash tubes. 
Two cameras were used for the photography ; one was placed with its axis 
on the axis of the chamber, while the other made an angle of 15° with it. 


* On leave from the University of Zagreb, Yugoslavia. 
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The wide angle between the cameras enabled the reprojection of photo- 
graphs to be done with great accuracy. The chamber was placed within a 
counter telescope which selected penetrating particles traversing , the 
chamber plates and, in addition, a 10 cm lead block placed above the 
chamber. Later, the whole counter telescope was placed under the 
chamber to allow the observations of possible secondaries produced by 


neutral primaries. 
Fig. 1 
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§ 3. EXPERIMENTAL RESULTS 


With the first arrangement (fig. 1) 6841 photographs were taken showing 
at least one ionizing particle traversing the lead plates inside and the 
absorber outside the chamber. In the search for single penetrating 
secondaries only events occurring inside the chamber and in the first four 
plates were considered. From 550m of lead traversed by penetrating 
particles no certain cases of single penetrating secondaries were observed 
The results of several other authors are shown in the table for Bonanno 

Since the production of penetrating particles found in the present 
experiment, and indeed in most of the recent work, is much less than that 
reported earlier, it is important to re-examine the former work for 
possible spurious effects simulating A.P.P. The present observations giv 
data which assist this reinterpretation. ae 


aot main contributors to the high apparent cross sections are likely to 


(1) Accidental coincidences between penetrating particles, and 
(2) knock-on electrons not recognized as such. 
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We found that single secondary particles occasionally penetrated one 
lead plate without multiplication or conspicuous scattering. Fourteen 
cases were observed. The projected angle of scattering in the 2 cm lead 
plate was less than 5° in six of these, while there were five angles from 10° 
to 15°. 


Metres of Pb | No. of single 
Authors Experimental method | of u-meson | penetrating 
track secondaries 


Braddick et al. Single-plate 78 4 
cloud chamber 


Braddick et al. Multiplate 
cloud chamber 


Lovati et al. Multiplate 
cloud chamber 


Deutschmann Multiplate 
cloud chamber 


George et al. Two-plate 
cloud chamber 


Appapillai et al. Single-plate 
cloud chamber in 
magnetic field 


Braddick and Leontic Multiplate 
cloud chamber 


From the range and estimated ionization of these particles none of them 
could have been heavier than 7-mesons. If they were 7-mesons with the 
observed range, five of them should on the average have been recognizable 
by ionization greater than twice the minimum. The energy of these 
particles considered, as 7- or u-mesons, could not have been greater than 
50-70 Mev and it seems unlikely that production would be confined to 
m-mesons of this energy range. The expected multiple scattering of 
mesons of this energy does not agree with that found in the experiment 
(fig. 2). It is therefore probable that all the particles were electrons. 

Pairs of particles have formerly been considered to be associated if 
they were reasonably coplanar, contemporary, and appearing to come from 
a point in an absorber lying inside or outside the chamber. In the present 
experiment direct observation inside the chamber was used to estimate 
the number of cases in which chance coincidences could simulate A.P.P.’s 
coming from a point outside the chamber. This was done by selecting 
tracks of penetrating particles that appeared contemporary and crossed 


4P2 
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m in a volume in the chamber corresponding in shape and 


within 10 m This criterion wa made to 


size to 10 cm absorber above the chamber. 
correspond as closely as possible to that use 
AP.P.’s. Nine cases satisfying it were found. 


d previously in selecting 


Fig. 2 
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Distribution of scattering of secondary particles compared with multiple 
scattering calculated for mesons of mass 276. 


The ‘spurious pairs’ obtained in this way, and the ‘ electrons pene- 
trating one lead plate ’ were subjected to an analysis similar to that used 
by Braddick, Nash and Wolfendale (1951). They plotted the multiple 
scattering distribution of the more and less scattered members of their 
A.P.P. pairs and found that the two curves differed from those obtained 
from ‘ artificial pairs’ composed of single particles varied at random. 
This was used in support of the existence of A.P.P.’s. In the present 
experiment, the scattering distribution of the pairs of tracks considered as 
chance coincidences simulating A.P.P.’s did not differ from that of 
‘ artificial pairs ’, but when the electrons simulating A.P.P.’s were included, 


the distributions were altered to agree with those obtained by B.N.W. for 
A:P-PYs, 
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If the results of the early work were in fact affected by the spurious 
events described above, then the cross section calculated by taking all 
these events as ‘ genuine cases ’ should be comparable to that obtained in 
those early experiments. If we therefore take the penetrating electrons 
and the chance coincidences together, the cross section so obtained is 
4-5 < 10-29 cm2/nucleon. Taking the chance coincidences only, this 
figure falls to 1:9 10-29 cm2/nucleon. Both of these cross sections are 
fairly near the values obtained by Braddick and Hensby (1939), Braddick 
et al. (1951) and George and Trent (1951) and some others, and it seems 
that the results of the above authors can indeed be accounted for in terms 
of the spurious events described. If these are to be avoided, only those 
single penetrating secondaries that are formed inside the visible volume j 
of the chamber and traverse at least two lead plates about 2 cm thick can 
be considered genuine. 

According to the data of the present experiment the upper limit of the 
cross section for the production of single penetrating secondaries is about 
210-39 cem?/nucleon. Since the cross section for the production of 
energetic stars by «-mesons is slightly higher than this, it is reasonable to 
assume that the production of single penetrating secondaries is analogous 
to this interaction. 

Three cases of nuclear interactions initiated by ionizing primaries were 
found. In these interactions several secondaries were produced by an 
energetic primary, presumably a p-meson. Some of these secondaries 
were heavily ionizing while others were minimum ionizing and penetrated 
several lead plates. Two of these events (fig. 3, Plate 36) occurred in the lead 
absorber above the chamber, while the other one occurred in the chamber. 
These interactions are considered to be analogous to energetic stars 
obtained in photographic emulsions underground by George and Evans 
(1950-51). The cross section found by us is also consistent with that 
‘obtained by these authors, being 3-2--1-8 x 10-39 em2/nucleon. In calcu- 
lating the cross section no geometrical corrections were applied, since it 
was considered unlikely that any event of this kind had escaped detection. 

In order to investigate the possible interactions of neutral primaries 
inside the chamber, giving rise to ionizing secondaries, the counter 
telescope was placed under the chamber and about 2,000 photographs 
taken with this arrangement. No definite cases of interactions involving 


neutral primaries were observed. 


§ 4, CONCLUSIONS AND Discussion 


Contrary to the results obtained by workers who used small cloud 
chambers equipped with only one or two plates, the data obtained in the 
course of the present investigation indicate that there is no justification for 
considering the formation of single penetrating secondaries by «-mesons to 
be of a different mechanism from that responsible for production of ener- 
getic nuclear disintegration stars produced by these particles. 


1292 On a Cloud Chamber Study of Cosmic Ray Interactions 


While this work was in progress, results of similar experiments were 
reported by Lovati et al. (1953) and Deutschmann ( 1953). The former 
worked at a depth equivalent to 55 m of water. Using a simple counter 
telescope as well as a split tray arrangement, they obtained several 
examples of nuclear interactions initiated presumably by pe-mesons. The 
latter obtained no events in which single penetrating secondaries were 
produced from 422m of lead traversed by p-mesons. Similar results 
were obtained by Walker (1953) who used a cloud chamber in a magnetic 
field. 

From the arguments exposed above it is apparent that the results 
_ obtained lately with multiplate chambers and based on events occurring 
inside them are more reliable. Although they show that single penetrating 
secondaries may occasionally be formed by p-mesons, this reaction is essenti- 
ally the same as those in which more energetic nuclear events are initiated 
by -mesons and the adoption of a separate term for it seems therefore 
unjustified. 
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SUMMARY 


- The disturbance produced behind a shock wave of arbitrary strength, 
travelling along a tube of varying cross section, is investigated. The 
shape of the cross section is arbitrary and the problem is linearized on the 
basis of small variations in area along the length of the tube. 

It is found that if the initial cross sectional area is S, and there is a net 
change in area of amount [4S], then the consequent change in pressure 
behind the shock is 

—K(py—Po)[AS]/S, 
where (p,—,) is the initial pressure discontinuity across the shock. 
The parameter K decreases monotonically with the shock strength and 
0-5>K>0-394. 

More generally it is shown that the pressure is affected by two distinct 
perturbations superimposed on the initial flow, namely a permanent 
modification arising directly from the variations in cross section, and a 
transient disturbance reflected from the shock and propagated with sonic 
velocity relative to the main flow. Expressions for both these contri- 
butions are obtained. 


$1. INTRODUCTION 


In a recent paper (Chester 1953) the author investigated the flow behind 
a shock wave travelling along a two-dimensional channel of non-uniform 
width. In particular, when the variations in width took place within a 
finite length of the channel, the ultimate change in strength of the shock 
wave was calculated together with the form of the disturbance which is 
reflected back along the channel. A different approach is given here, and 
is applied to a tube whose cross section is quite arbitrary in shape, the 
perturbations on the flow behind the shock wave being produced by small 
variations in the cross sectional area. It is found that, ultimately, the 
disturbance depends only on the variations in the area of the tube and not 
on the actual shape of the cross section. 


§2. TRANSITION RELATIONS ACROSS A Unirorm SHock WAVE 


For future reference it is convenient to quote here the transition 
relations for a plane shock wave of uniform strength moving into a fluid at 


rest. 


* Communicated by the Author. 
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Let the shock velocity be U, and let the fluid velocity, pressure, density 
and sonic velocity in the region of uniform flow behind the shock be 
di» Pris Pr» 4% respectively. The corresponding quantities ahead of the 
shock will be denoted by 0, %, po: %- It is also convenient to introduce 
the symbol U,=(U )—q,) for the velocity of the shock relative to the flow 
behind it. 

If y is the adiabatic index and My=Uo/dp is the Mach number of the 
shock, we then have the following relations (see, for example, Courant 
and Friedrichs 1948) : 


2U. 
Pebeaen hooey th « + oS ec 
U,;=U ao 2M)? 
ay yal Feerrargii Plat fi Te} esl 
2 —l 
te P 
P= oo pol Ut as} = PBS (2 Mey +} eee 
Pee as y—1+2M,-? 1/2 
WMS a se | a 3 . . . . . . (4) 
Gla Mane 
=a ane eel ee © ee Sei a a 


§3. THe BounpaRy CONDITIONS aT A Non-UNIFORM SHOCK 
We next consider the effect of a small perturbation superimposed on the 
flow behind the otherwise uniform shock wave. 
Let (x’, r,@) be cylindrical polar co-ordinates relative to a fixed origin 
with the x’-axis in the direction of propagation of the shock. Relative to 
these axes, let the fluid velocity and pressure behind the shock be res- 


pectively V,=(q,+¥, v, w) and (p,+p), and let the normal velocity of the 
shock be Vy. Then the transition relations are 


2V 7 2 = 
vie obi _ a _ 2Po0 y—! 
1 yl 1 ra ? Pix = Plt So 5 (6) 


X 


and the equation of the shock may be written 

a’ —Ugt—fir, 6, t)=0, “so Wee: eae sete 
where f is assumed to be uniformly small. 
Correct to the first order of small quantities, V,— i 

, g , Vo=(Uo+f;,) and the unit 
normal at a point (r,@) of the shock is (1, —f,, aay It follows that 
ee. Wo= (oth ores —Uofeir)- xe Gta nvm (A) 
Substitution of relation (8) in (6) yields the following conditions to be 


satisfied at the shock front : 
2 
= y+1 (eis v=—4y),, w=—q, far, 


_ 4p0Uo 4p,U, ait) 
deere Ga lie Srreesy 
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Relations equivalent to these were first obtained by Ting and Ludloff 
(1952) and applied to a problem involving the diffraction of shock waves. 


$4. Tae Equations oF MOTION IN THE PERTURBED FLOW 


It is more convenient to derive the equations of motion relative to axes 
which move with the mean flow behind the shock wave. The co-ordinates 
relative to these axes will be denoted by (x, 7, @). The velocity vector is 
now simply (wu, v, w)=V (say) and the equations of conservation of mass, 
momentum and entropy, when linearized on the basis of small pertur- 
bations, are respectively 


0 
+ Lapesa) EUS, me nae trees KI) 
Ne ul 
ie ae 5 eae shee dwe) fel) ta, 1s ( ) 
os 
a acts SE. et eee) 


Here and elsewhere dependent variables without a suffix refer to the 

perturbed flow superimposed on the uniform flow denoted by suffix 1. 
Equation (10) with the aid of eqn. (12) may be written 

ae 

a,? ot 

where a,2—=(0p/dp)g is the square of the velocity of sound in the uniform 


flow. 
Between eqns. (11) and (13) the velocity V may be eliminated to yield 


the following equation for the pressure : 


Le GEED 
ay) 
Nee a Cia 


+piv.V=0, Rebs 22 (18) 


(14) 


Equation (14) implies that the propagation of small variations in 
pressure is still governed by the wave equation, notwithstanding the fact 
that there are entropy variations (and hence vorticity) present in the flow. 
It is not, however, true of the velocity components, as it would be in a 
genuine acoustical problem. This important fact was first realized and 


exploited by Lighthill (1949). 


§5. APPLICATION TO FLOW IN A SHock TUBE 


We now consider the application of the preceding results to the parti- 
cular problem of flow in a quasi-cylinder. The boundary of the quasi- 


cylinder is defined by the equation 
r= R(x’, 0)=R(x+qit, 9), ne. ae LO) 


where R is arbitrary except for the restriction that R, should be small. 
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Since the unit normal at the boundary has components 
fis, —R,/R, 1}R/(R2+ R,?), 
it follows that the condition to be satisfied there is 
RR? 


Vos qi (R2--R 2)e? (16) 
where V,, is the normal velocity. Since, by eqn. (11) 
oV 1 dp 
= SS es Oe td 
ot P1 on ; ( ) 


the corresponding condition for the pressure is that, at 7=R, 


ap Pda 0? 2 
et — GEER age (OP) . woe ga Ca 
correct to the first order. 

Now ultimately we may expect the pressure to become uniform over 
the cross section of the tube, and restriction of the investigation to the 
average pressure leads to a significant simplification in the analysis. We 
make use of this simplification in the following way. The two-dimensional 
form of the divergence theorem is applied to yp over the cross section to 


get 

Op Ot Ache ee Oy Oc 2) 

|e | [(Sa+3 r ar r2 002 Est phe circ 

Since, on a linear theory, it is permissible to apply the boundary 

condition (18) in the plane of the cross section, the left-hand side of 

eqn. (19) may be written 
0? " 

— pid" Ox (ZR?) d6=—pqPAS"(x+qyt), . . . (20) 

where AS has been written for the variation in cross section relative to 
some reference cross section.* 


Furthermore, since p satisfies the wave equation, 
LCP veel Ghe ae 1 Op ap 
I] (Gat ror’ pn) S= {J (= on =) as. (21) 


=| | pds, 
where the integral is taken over the cross section, it follows that 
OPRaalGorr 
dat a2 oe Pdr AS" (etait). an C2) 
The solution may be written in the form 


P=Ap,q,"[F(x+ayt)+ AS(a+q,t)], Seer)! 
eee 


* The entropy, for which the symbol S was al i 
the subsequent analysis. ; aiapeagasnncat es rege 


Hence, if 
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where a," 1 

A= a i ? 

ay7—qy 1—m? 

and F' is some function to be determined. The contribution to P from 
the term involving AS is the variation in pressure due directly to the 
change in cross section. This disturbance is reflected at the shock and 
gives rise to the term involving F(x+a,t). No contribution of the form 
G(«—a,t) is included since there is no mechanism for reflection upstream 
of the shock which could give rise to such a term. 

To obtain F we require the boundary condition to be satisfied at the 
shock wave, where eqns. (9) must be satisfied and from which a condition 
can be derived to be satisfied by P 

We note first that, as a consequence of (9), 


ap 4p,U, aV, 


(24) 


ee 25 
an ay) a’ . 
at x=U,t and r=R. 
With the help of eqn. (16), this gives 


an EL) CRE R GE ds Os) 
when «=U jt. 
If now we apply eqn. (19) in the plane of the shock, the left hand side will be 


me U. 
y+l 


where (U,+q,) has been replaced by Uy. To evaluate the integrand on 
the right of (1 ae we have, from a ) and (13), 


1 
rs ot naples re a {ti fo}= {u Aas piGy2 v4 ? : ; : (28) 


at x=U,t. Hence, by differentiation, 


eS Uat es oe tee. 2-6 @-. (27) 


rax(y+}) 1 not = aA J | 99 
oe 4p,U, Pr rica 4p,U, Di AF OS pee! : ( ) 
But, by (9) and (11), 
du 1 BA Ae (1+M,*) dp 
dt =Uy+ Uw, =a pee ae 5. Lane dt 
1+M,? 

= wet Uae) Og se eer 621) 
Hence SUS). Le . en 
Ug 29,0," {o, + Up, }+ PU ee ae ( ) 


and (29) becomes 
Tomes etal 4 F) a } 
Pe dommes eis LOB +5 Bee 


=9 a 2 
{EE (et Uap.) baht (n+ Fr.)}, 2) 
at «=U ,t. 
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This relation may now be used as the integrand on the right of (19). 
When combined with (27) the following boundary condition is obtained, 


to be satisfied at x=U,t: 


d) (1+, a? 
15 +U, ae (Pi+ U,P,)+uU/ (P+ iP.) 


2 
=— mn P1917 48" (Ut). . ° . . . (33) 
1 


The function F'(x--a,t) is now found from the fact that the relation (23) 
for P should satisfy (33). It follows at once that /’ is of the form 


—BAS{r(x+a,1)}, 
where 
es Us oa M,+m 34 
7] Up-at a ie a) 
so that 
P=Ap,q,[—BAS{r(a+a,t) }4 48S(a-+qy8)], Ps es Oh 


=— Ap q,[BAS {r(x' + (a,—9,)t)}7— AS(2')]. . + (36) 

Substitution of (35) in (33) now gives 

1+M,?+2M,m-1 
OF TEM aI 

Equation (36) illustrates clearly the two distinct perturbations super- 
imposed on the initial flow, namely a permanent modification from the 
term Ap,q,?4S(x’) arising from the variations in cross section, and a 
transient disturbance, —ABp,q,24S{7(x’+(a,—¢q,)t)}. The latter is 
propagated with sonic velocity relative to the fluid and is due to reflection 
from the shock of the perturbations in the flow caused directly by the 
variations in cross section. 

When the tube consists of two cylinders connected by a transition 
section, and the uniform flow denoted by suffix 1 refers to the upstream 
section, the ultimate pressure change behind the shock will be given by 


Po—Pi=—A[B— 1]p1917[ 48]/S 
=— K(py—p)[AS]/S 
(say), where S is the cross sectional area in the upstream section, [4S] is 
the net change in area and 
Kix ae ts La (39) 
(I--m) (1M, ?4+2M,)? = 


use being made of the relation 


(37) 


(38) 


Pi—Po=P11U 4. 

The parameters m, M, and M, are all functions of P,(Po through 
eqns. (3), (4) and (5), and the variation of K with this parameter is shown 
in fig. 1. It is a monotonic decreasing function with upper and lower 
limits respectively 0-5 and 0:394.... It is interesting that the pressure 
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increment has the same sign for both subsonic and supersonic flow 
behind the shock (m21) even though the contribution arising directly 
from the variations in cross section is of different sign for the two regimes 
(for there is a change of sign in the parameter A). This implies that to 
increase the Mach number behind the shock, even in supersonic flow, 
requires a contraction, as opposed to the expansion necessary in steady 
supersonic flow in a duct. 


Fig. 1 


eee ee 160 ATS. 49 10 
Pi/Po 


Variation of K with p,/po. (y=1"4) 


For subsonic flow behind the shock a disturbance will be reflected into 
the upstream section of the tube, and the pressure variation in the pulse 
will ultimately be given by 


feat Oe —ABp gqPASs (rz + (a,—4,)t) }/S 
= —N(py—9) AS {r(2' + (41-91) ) 18 » « (40) 


(say), where 


=~, ~~ MQ m+ Mo? + 2M) © ; 


But, when the flow behind the shock is supersonic, this pulse will be 
convected into the downstream section of the tube and the whole of 
eqn. (36) is required to obtain the average pressure increment, since 
AS(z') is non-zero in the downstream section. Furthermore the assump- 
tion that the pressure will ultimately be uniform over a cross section is no 
longer true in the supersonic case (except in so far as neglected effects, 
such as the action of viscosity, will tend to produce such a state), for the 
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full solution contains a fluctuating contribution which does not die out 
even at distances remote from the origin of the disturbance. However, 
the average pressure is a reasonable guide to an understanding of the state 
of the fluid and, in any case, the fluctuations do not penetrate downstream 
as far as the shock, so that this criticism applies only to the rearward pulse 
(for further details see Chester 1953). 

The values of t and N for various shock strengths are shown in figs. 2 
and 3 respectively. The parameter 7 is monotonic decreasing and takes 


Variation of 7 with p,/py.  (y=1'4) 


the values 2, 1 and 0-6076... respectively when the incident shock has zero 
strength, when the flow behind the shock is sonic and when the shock 
strength becomes infinite. 

For subsonic flow behind the shock N increases rapidly from a value of 
0-5 for a shock of zero strength, and becomes singular when the flow 
behind the shock is sonic. For supersonic flow N becomes negative (and 
singular at m=1) and increases monotonically to an asymptotic value 
of —0-1550.... 

The singularity in N at m=1 (which arises from the singularity in A) 
does not imply a similar behaviour in the solution for the pressure. For, 
by (37), B=1 when. m=1 and hence, in eqn. (36), the term in square 
brackets is zero. More precisely, if m=1-+8, then A=0(8"1) and the 
term in square brackets is 0(5), at least for moderate values of the time. 
But the tendency for A to become large for transonic flow can have an 

appreciable effect on the pressure when f is so large that 


AS {7(x’ + (a,—4,)t)} 
differs appreciably from 4S(a’). For subsonic flow (q,<4,) this can only 
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occur in and upstream of the transition section, for (36) is equivalent to 
(38) in the downstream section. On the other hand, if ¢,>a,, these large 
disturbances are excluded from the upstream section. This tendency 
of the disturbances to accumulate will undoubtedly give rise to non-linear 


Fig. 3 


me 
oe] a a 
Pine g 


Variation of N with p,/po. (y=1°4) 


phenomena such as shock waves and in practice becomes manifest as the 
familiar process of ‘ choking’. These remarks, of course, apply only to 
the rearward pulse and not to the region immediately behind the main 
shock wave where eqn. (38) is always valid. 
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ABSTRACT 


A multiple beam interferometer used for detecting liquid helium in 
thin gaps is described. Experiments performed with this apparatus 
have led to the conclusion that helium II will enter a narrow gap 
maintained above the liquid surface, providing the gap is bridged at a 
point below the surface tension height. A plausible mechanism is 
proposed and the bearing of the conclusions on experiments with liquid 
helium IT discussed. 


$1. INTRODUCTION 


EXPERIMENTS have recently been reported by Lane and Dyba (1953) 
concerning the rising of liquid helium II into narrow gaps by an 
extraordinary capillary action. Two parallel plates held in close 
juxtaposition separated by spacers were maintained above the liquid 
surface and in contact with it by means of the helium II film. These 
authors observed liquid to fill the gap between the plates up to the 
surface tension equilibrium height, and found further evidence to support 
the hypothesis that the film was the agent transmitting hydrostatic 
pressure. This is reasonable, yet. the mechanism by which liquid 
transported through the film is first established in an empty gap is very 
difficult to understand. 

Perhaps the simplest interpretation would be to assume that the 
appearance of liquid was in some way connected with the spacers which 
directly bridged the gap. However, experiments performed some time 
ago in this laboratory (Hatton and Rollin, unpublishedt) seemed to 
indicate that the phenomenon occurred even when the gap was not 
directly bridged by spacers. The experiments described in this paper 
using an optical method especially adapted to conditions in liquid 
helium, were therefore made in an effort to clarify the situation. 


§2. MeTHOD 


Liquid helium has a very low refractive index and a small absorption 
coefficient in the optical region. This means that observation of the 
meniscus between liquid and vapour in small apertures is particularly 
ee EE See eee 
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The interferometer illuminated with green mercury light. The letters refer to 
the parts indicated in fig. 1. The gap | has been filled through the film ; 
the meniscus is at G. The other two gaps m and n are empty. The 
wire bridging the gap is indicated by R. The transmission in this 
region has been increased by the removal of several small flakes of 
silver from the glass owing to the wire moving as the interferometer 
was assembled. The small variations of the gap width indicated by the 
fringes were caused by the rubbing down ofthe edges with jewellers’ 


rouge. 


N. G. McCRUM Phil. Mag. Ser. 7, Vol. 45, Pl. 38. 
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The interferometer illuminated with green mercury light. All three gaps 1, 
m and n are filled with helium, the menisci being at G. To fill the 


gaps m and n it was necessary to lower the interfer ometer until the 
edge AA touched the surface of the bath liquid, 
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difficult. As it was expected that minute transient volumes of helium 
were to be detected, an interferometric method was devised capable of 
doing this. 

Fig. 1 


PLATE Y PIRATE? eX 


The interferometer plates. 


Fig. 2 


The interferometer assembly. 
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The apparatus consisted of two optically polished glass plates X and Y, 
held together with springs S and separated by wire spacers W, as shown 
in figs. 1 and 2. In order to eliminate edge effects, horizontal grooves 
were cut into X so as to separate the area under observation from the 
spacers and to prevent liquid reaching the gap by ordinary capillary 
action. The edges of the grooves after the initial rough cutting were 
rubbed down with jewellers’ rouge to remove small chip marks and 
abrasive scratches. The lower spring bore a brass strut T as shown in 
fig. 2 which supported the apparatus vertically in the cryostat. The 
strut was bolted to a cylindrical wire cage which could be moved up and 
down inside the cryostat through the liquid helium, by means of a waxed 
hemp line and a windlass operated by turning a glass cone in a ground 
joint. Thus, the mechanical support for the glass plates passed through 
the helium bath. This method of support ensured that the heat leak 
into the glass plates was extremely small, thereby simplifying the 
interpretation of the results. 

By depositing a thin film of silver on each of the plates, the assembly 
was turned into a multiple beam interferometer. With silvering of the 
order of 90°% reflection and with the springing carefully adjusted three 
or four well resolved interference fringes per cm were observed. Very 
small droplets of helium were usually detected with unfiltered mercury 
light, changes of optical path then producing very clearly defined 
variations of tint. When using monochromatic light, liquid was dis- 
tinguished from vapour by means of a break in the interference fringes. 
The green mercury line or the sodium doublet was used depending on 
which wavelength gave the better fringe separation. The actual fringe 
separation was visually estimated and the interferometer spacing 
calculated, a rough knowledge of the spacer width resolving the 
ambiguity in the order of interference. Defects in the surface, plate 
distortion caused by the springs or thermal contraction and small 
foreign particles were all easily observed. 

An optical bench consisting of a monochromatic source, a collimating 
lens, several layers of Chances ‘ON 20’ glass and, when necessary, 
suitable Wratten filters, was mounted on one side of the cryostat so as 
to introduce light normally on to one of the faces of the interferometer. 
General features of the formation of liquid between the plates were 
conveniently observed with the naked eye held close to the opposite 
side of the apparatus, but small detail was detected by means of a 
low-powered microscope with x8 magnification. 


§3. EXPERIMENTS 


The first experiments were made with an interferometer formed from 
plates similar to X and Y (fig. 1). The horizontal groove AB in X, was 
in this case 1:3mm deep and 3-5 mm high, having a roughly triangular 
cross section. The other plate was identical to Y except that it did not 
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contain the vertical grooves C and D. It was then found that the empty 

interferometer could be filled with helium II but not with helium I, 

when the liquid surface was below the edge AA. The filling was done 

either by varying the height of the interferometer isothermally or by 

lowering the temperature to below the A-point with the edge AA above 

the surface yet below the equilibrium surface tension height appropriate 

to the gap width. The interferometer was emptied by raising it clear 

of the liquid. The maximum height of the edge AA above the free 

surface when on lowering the interferometer, the helium entered the 

empty gap, the entry height was less than the height of AA above the | 
surface, the exit height, when on raising, the last liquid left the gap. 

In this experiment the heights were 0-65 cm and 1-6 cm and the gap 36 p. 
The exit height was equal to the surface tension equilibrium height 

(Allen and Misener 1938) appropriate to the gap width. 

The shape and proportions of the cross section of the triangular 
groove AB were not considered satisfactory and it was therefore cut 
with a rectangular cross section, being 5mm high and 2mm deep. 
The plates were then resilvered, assembled and placed in the cryostat. 
In the first experiment, on cooling to helium temperatures, the plate X 
cracked horizontally at F forming two independently sprung segments. 
The upper segment formed a small angled wedge with the plate Y. The 
fringes were horizontal and closely spaced and on counting it appeared 
that the gap varied from approximately 30, at the top to 50m at the 
edge AA. The exit height was approximately 1cm and the entry 
height zero; the empty interferometer could now only be filled by 
lowering it until the surface touched the edge AA. 

For the next experiment a new plate was cut similar to X. It was 
mounted with plate Y in the interferometer, the gap being throughout 
30 wide. The entry height was approximately half the exit height 
which was 14mm. The liquid on entering the slit appeared first as 
a very small droplet just above the edge AA. Its initial growth was 
slow but increased as the perimeter increased, as would be expected if 
helium were being transported to it through the film. The time taken 
from the first appearance of the initial droplet to the final equilibrium 
position was of the order of seconds, being shorter the lower the 
temperature. 

The next experiment was conducted with the same assembly as the 
last one, but the gap reduced to 12 by using a finer tungsten wire 
spacer. A short length of this wire, U, was also placed horizontally 
across the gap as indicated in fig. 3. No precautions were taken to 
exclude dust particles and a careful examination of the surfaces revealed 
that the silvering, which was that of the previous experiment, had 
deteriorated. Small areas of silver were detached from the glass and 
lay as minute flakes between the two surfaces. In this experiment 
detailed observations were made of the formation of helium IJ in the 


gap. 
4Q2 
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On lowering the empty interferometer until the helium II surface was 
in the region of KK (fig. 1) the first sign of liquid was seen between the 
plates as shown in fig. 3 (A). This droplet then grew slowly out towards 
the edges. Subsequently a small drop appeared in the region L. The 
silvering appeared badly flaked in this region. This drop although of 
a circular section initially, grew into a pear shape and then moved down 
as indicated. The two volumes of helium then joined as shown in 
(3C+D). At the same time small droplets were appearing higher up 
in the gap, growing and then moving down either to the rising meniscus 
or to the wire (3H, F+G). The final equilibrium position is shown 


Fig. 3 
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Successive stages in the establishment of liqui i i 
, iquid helium II in an empty ga 
mace contains small quantities of bridging material. The nie aes 
rom the initial stage A, to the final stage H was of the order of seconds. 
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in fig. 3(H). The distance from the meniscus at this position to the bath 
surface was now equal to the appropriate surface tension height. It was 
noted that droplets apparently grew where the silvering was badly flaked 
or at the wire and that no droplets ever occurred above the equilibrium 
height. The inference was that liquid would not enter an empty gap 
unless there was a nucleation centre below the appropriate surface tension 
height. 

The next experiment was designed to test this hypothesis. Plates X 
and Y were assembled freshly silvered with a 30 spacing. A close 
inspection of the surface after liquefaction revealed a very small particle 
of dust at the point P (fig. 1). Isothermal variations of the height of 
the empty interferometer showed that after filling at the maximum 
entry height the meniscus in the gap, when finally in equilibrium, passed 
through the enclosed particle at P. If the empty interferometer was 
raised a small distance of the order of a millimetre above this entry 
height the empty gap could be filled by tapping the apparatus, thus 
causing the bath surface to oscillate. It would not fill at this height 
unless tapped and at greater heights would not fill under any circumstance. 
Further evidence was obtained by carefully stripping the interferometer 
so as not to injure the silvering and placing a small piece of wire about 
2 millimetres long and of the same cross section as the spacer at the 
point Q (fig. 1). The plates were then replaced and the assembly placed 
in the liquefier. Great care was taken to prevent dust particles fouling 
the surfaces in order that the experimental conditions should be identical 
with those of the preceding experiment with the exception of the small 
piece of wire at Q. As was to be expected the maximum entry height 
was now larger than in the former experiments. The final equilibrium 
position of the meniscus in the interferometer when filled at this height 
coincided with the wire at Q. 

For the concluding experiment a plate similar to Y as shown in fig. 1 
was constructed. Two vertical grooves C and D, 2mm wide and 2mm 
deep were cut and the edges polished with jewellers’ rouge. This plate 
was then assembled in the interferometer with plate X and separated by 
30 wire spacers. This method of mounting provided three identical, 
independent gaps 1, m and n. Gaps m and n were kept free of any 
bridging material but a small piece of wire was placed at R in slit 1 just 
above the edge AA. It was found that liquid helium II would not enter 
the empty gaps’m and n unless the bath surface was brought into direct 
contact with the edge AA. Gap | on the other hand, could be filled 
through the film. If gap | was filled so that the entry height was a 
maximum, the meniscus of the liquid in the gap when finally in 
equilibrium coincided with the wire. Photographs of the appearance 
of the interferometer under these conditions illuminated with green 
mercury light are shown in Plates 37 and 38. 

Since in these more detailed experiments a bridge at or below the 
surface tension height was found to have caused the appearance of 
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liquid between the plates it must be assumed that the formation of the 
droplets in the first and third experiments was due to the presence of 
unobserved bridging material. The positions of these bridges were 
evidently at 1-6 cm and 1:4 cm above the bath surface when the empty 
interferometer was filled at the maximum entry height. 


§ 4. Discussion 


Kelvin has discussed a hypothetical experiment bearing a marked 
resemblance to that just described (Thomson 1871). He considers a 
closed system containing two detached volumes of the same liquid, 
one with a large meniscus and the other in a small capillary, in communi- 
cation with each other only through the vapour phase. The capillary 
contains insufficient liquid to raise the meniscus to the height that 
would be attained if it were in direct hydrostatic communication with 
the large volume of liquid. It is argued that given sufficient time and 
some thermal contact between the two volumes of liquid, equilibrium 
would finally be achieved by the condensation of vapour into the 
capillary. The experiments reported in the previous section indicate 
that if the liquid were helium below the A-point, equilibrium would 
certainly be established quickly since the mass transfer would take 
place through the film acting as a hydrostatic link. The permanence of 
equilibrium between detached volumes of liquid implies the maintenance 
of equilibrium between the vapour and all the various liquid menisci. ° 

In common with other liquids, a concave liquid helium surface is in 
equilibrium with a vapour at a pressure lower than the equilibrium 
‘pressure for a plane meniscus maintained at the same temperature. 
The relation between p, the equilibrium pressure over the curved surface 
of radius r, py the equilibrium pressure for a flat surface, o the surface 
tension and V the molar volume of the liquid at temperature 7’ is given 
by the Kelvin equation, 


In Sj 9 a aes 


At a height H above the meniscus of the bulk liquid the pressure will 
be given 


P=Dp exp (—MgH/RT), PS ies sy OL ha BF) 


where M is the molecular weight, assuming the vapour obeys the ideal 
gas laws. Let us assume that there are a set of small irregularities of 
varying radii on an object projecting through the surface of a liquid in a 
closed vessel. If the object is lowered into the bulk liquid the height of 
any one of these irregularities above the surface is decreased. Consequently 
the pressure p and therefore the potential for capillary condensation are 
increased (eqns. (1) and (2)). Shown in fig. 4 is a schematic diagram 
illustrating the proposed mechanism by which a drop of liquid helium IT 
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is first established in an empty gap between two parallel plates. Bridging 
material is enclosed in the gap at B. As the plates are lowered into the 
liquid helium the height of B above the surface S is varied from H, to 
H, to H;. The corresponding liquid radii of curvature R, kh, R, 
respectively are in equilibrium at these heights. Once the radius &, 
is established then liquid from the film bridges the whole gap, the drop 
falls to the lower part of the gap at C and then proceeds to expand. 
As this case corresponds to filling at the maximum entry height, the 
meniscus, when finally in equilibrium, will coincide with B and H, will 
be the surface tension height appropriate to the gap width. It must 
be presumed that similar equilibrium conditions are never attained 
above the A-point owing to the immobility of the helium I film and the 


A schematic diagram illustrating the formation of liquid helium IT in a gap 
containing bridging material, B. 


relative inefficiency of mass transport through the gas phase. For 
under the experimental conditions previously described, the gap, when 
filled with helium II, could be emptied by raising the temperature above 
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the A-point ; the helium I would slowly evaporate away owing to the 
absorption of radiant energy. Presumably the isothermal conditions 
existing in the helium II region due to the mobile film are destroyed at 
-point. 

rca der Waals interaction of helium atoms with an ideally flat 
substrate has been estimated to permit a film of the order of 400A to 
exist above the liquid meniscus (Frenkel 1940, Schiff 1941). This liquid 
will be known as the film liquid. However, since it is very difficult to 
get a truly flat surface (Thomson 1953) it must be assumed that on any 
surface other than perhaps a cleavage plane, the total volume of adsorbed 
liquid will exceed the film liquid. This additional liquid will be known 
as the capillary liquid. The volume of capillary liquid will be a function 
of the microstate of the surface and will be temperature dependent. 
The experiments previously described indicate that the volume of capillary 
liquid will probably change discontinuously at the A-point, owing to 
non-equilibrium conditions in the helium I range. Schiff (loc. czt.) 
suggests that the volume of film liquid will change discontinuously for 
similar reasons. These experiments emphasize the significance of the 
microstate of the surface of vessels used in film flow experiments 
(Chandrasekhar and Mendelssohn 1952). It is possible that the wide 
divergence in measurements of the He II film thickness (Kikoin and 
Lazarew 1938, Daunt and Mendelssohn 1939, Burge and Jackson 1950, 
Bowers 1953, McCrum and Mendelssohn 1954) may be, in part, due to 
the use of different substrate. 

Most probably in the experiments of Hatton and Rollin small particles 
were trapped in the gap thus forming suitable nucleation centres for the 
formation of the initial droplet. The foil spacers used by Lane and Dyba 
(loc. cit.) would presumably also initiate the entry of liquid. Bowers (1953) 
has described difficulties encountered in using a loosely wound foil to 
measure the thickness of the helium II film. It appeared that at satura- 
tion the formation and retention of liquid between the folds rendered a 
precise determination of the thickness impossible. This was probably 
due to the surface tension rise of liquid into gaps lying between the 
foils. The nucleation centres would be provided either by direct touches 
of two foil surfaces or by bridging material. 

Lane and Dyba (loc. cit.) have already pointed out that the curious 
tendency of helium to form in narrow gaps above the liquid level, will 
have to be considered in experiments below the )-point. The present 
work must certainly serve to emphasize this caution. In addition 
however, these observations have now elucidated the mechanism by 
which this undesired formation of liquid takes place. The safest way to 
avoid such formation will be to see that the experimental arrangements 
used in future work do not contain constrictions in which capillary 
condensation may occur. If gaps of constrictions below the surface 
tension height have to be employed they will only stay free of liquid if 
bridging material such as dust or condensed gases are rigorously excluded. 


of Liquid Helium in Narrow Gaps 1311 


ACKNOWLEDGMENTS 


I would like to thank Dr. Mendelssohn for the interest he has shown 
in this work. I am grateful to Dr. D. F. Brewer for several valuable 
discussions and Mr. D. O. Edwards for assistance with the experiments. 


REFERENCES 


Atuen, J. F., and Misener, A. D., 1938, Proc. Camb. Phil. Soc., 34, 299. 

Bowers, R., 1953, Phil. Mag., 44, 1309. 

Burce, E. ie and Jackson, L. C., 1949, Proc. Roy. Soc. A, 205, 270. 

Daun, J. G., and MENDELSSOHN, ee 1939, Proc. Roy. Soc. A, 170, 423. 

GHANDRASEKHAR, B.S., and MrNDELSSOBN, K., 1952, Proc. Phys. Soc. A, 65, 
226. 

FRENKEL, J., 1940, J. Phys. U.S.S.R., 2, 365. 

KIKoIN, A. ie and LazAREw, B. G., 1938, Nature, Lond., 141, 912; 142, 289. 

LANE, G. 405. and Dysa, R. V., 1953, Phys. Rev., 92, 829. 

McCrum, N. G., and MrNDELSSOEN, K., 1954, Phil. Mag., 45, 102. 

ScuiFr, L. I., 1941, Phys. Rev., 59, 839. 

THOMSON, W., 1871, Phil. Mag., [4] 42, 448. 

THomson, G. P., 1953, Structure and Properties of Solid Surfaces (Chicago 
Gomer and Smith). 


(mela heer 


CXLVIL. The Alpha-Particle Component of the Cosmic Radiation 


By C. J. WADDINGTON 
H. H. Wills Physical Laboratory, University of Bristol * 


[Received July 14, 1954] 


SUMMARY 


In an investigation of the primary «-particle component of the cosmic 
radiation a nuclear emulsion technique has been employed. A value of 
60-17 g/em? was found for the collision mean free path of o-particles in 
emulsion. The justification for using this value to obtain values for the 
mean free paths in other media is discussed. 

A value of 320-+ 36 «-particles per m? per steradian per sec was found 
for the primary flux at a geomagnetic latitude of 55°N. This was in 
agreement with previous values. The energy distribution found was in 
agreement with that for the heavy primaries (Z>3) given by Dainton et 
al. (1952) and Gottstein (1954). 


§ 1. INTRODUCTION 


THE presence of multiply charged particles in the primary cosmic radiation 
was first reported by Freier et al. (1948), who showed that there were highly 
energetic nuclei of elements with atomic number as high as 20 incident on 
the earth’s atmosphere. Since then, mainly by the use of nuclear emulsion 
techniques, a number of investigations have been made of the various 
properties of these ‘ heavy primaries ’. In particular, attention has been 
directed towards the problems of determining their fluxes, energy distri- 
bution, charge distributions, and mean free paths. These properties have 
been shown to be of fundamental importance in any consideration of the 
possible origins of the cosmic radiation (Bradt and Peters 1950 a, b, 
Dainton et al. 1952). Some of the results obtained have, however, been 
conflicting, in particular those concerning the energy and charge distri- 
butions. 

In an attempt to clarify some of these problems an investigation has 
been undertaken of the primary «-particle component of the cosmic 
radiation, Values have been found for the mean free path in emulsion, for 
the primary flux at a geomagnetic latitude of 55°N, and for the energy 
distribution up to 3 Bev per nucleon. These results have been compared 
with those found previously (for the heavy primaries), and they have been 
shown to be consistent with those obtained by other workers using the same 


emulsions but different techniques (Dainton etal. 1952, and Gottstein 1954), 
ee aL en 


* Communicated by the Author. 


On the Alpha-Particle Component of the Cosmic Radiation 1313 


§ 2. EXPERIMENTAL PROCEDURE 


A stack of glass-backed 400u Ilford G5 emulsions was flown for four 
hours at an altitude of 98000 feet at a geomagnetic latitude of 55°N 
(corresponding for «-particles to a geomagnetic cut-off of 0-33 Bev per 
nucleon). This stack was the same as that used by Dainton e¢ al. and 
Gottstein in their investigations of the heavy primaries. It was arranged 
in such a manner that half the emulsions faced each other, with air gaps of 
approximately 1mm between them. 


(a) Track Selection 


A number of these ‘ facing ’ plates, which were particularly selected for 
their lack of distortion, were scanned for tracks satisfying the following 
criteria : 

(i) That the track should have a length in the emulsion greater than 
8mm. This eliminated from consideration all the slow singly charged 
particles, and ensured that there was sufficient length of track to allow an 
attempt to be made to determine the scattering parameter a. 

(ii) That the track should have a grain density greater than three times 
the plateau value. This criterion eliminated all the fast singly charged 
particles, while including all those with multiple charge. 

The plates were scanned along horizontal lines parallel to the top edges. 
These scanning lines were separated from each other by a distance of 
10 mm, and from the edges of the plates by a distance of at least 10 mm. 
When allowance was made for tracks crossing two or more scans, each 
scan could be considered as being independent. Thus the scanning was 
along a line rather than over an area, and the expected length distribution 
was of the form dl//? rather than dl/l?. 

Using these selection criteria in the scanning, there were a number of 
possible ways in which a track could have been missed due to inefficient 
scanning. Because of this possibility various checks were made on the 
experimental data to ensure that no systematic errors had been introduced. 

(i) Tracks with ranges only just in excess of the minimum length 
accepted (8mm) could have been missed. Whether this had in fact 
occurred was checked by plotting the length distribution of all the tracks 
found. This showed that, to within the statistical limits, few tracks had 
been missed due to this effect. 

(ii) Tracks with grain densities only just in excess of the minimum 
accepted could have been missed due to an underestimate of their grain 
densities. This possibility was checked by comparing the grain density 
distribution of the singly charged particles with that found by Camerini et 
al. (1950) for the singly charged particles ejected from stars. There were 
no apparent differences between these distributions which could be attri- 
buted to the missing of tracks. The statistical fluctuations were large 
however, because of the limited number of events, and the possibility 
could not be eliminated that a proportion of those particles with the 
minimum acceptable grain density had in fact been missed. However, 
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since the «-particle minimum grain density was 1-3 times the minimum 
accepted grain density, the possibility that any «-particles had been missed 
due to this effect was certainly small. 

(iii) There could have been a systematic missing of tracks, of all grain 
densities and all ranges, but of the fifty-four cases where a track was 
observed in one scan and was expected to be observed in another, only 
one was missed in a manner which could not be explained, indicating that 
any effect of systematic missing was negligible. 

(iv) Tracks which appeared in the scan to be steeply dipping and which 
were therefore neglected by the scanner, could have been scattered into 
the plane of the emulsion elsewhere along their paths, and have had a 
range greater than the minimum accepted. Calculation shows that this 
effect would only be of importance for those particles having a large 
scattering parameter, and hence, although it is probable that a certain 
number of the slow singly charged particles were missed due to this effect, 
the number of «-particles missed was extremely small. In addition, 
tracks with high grain densities and thus probably high scattering 
parameters were examined, even if they appeared to be quite steeply 
dipping. In fact one track, not that of an «-particle, was missed on a second 
scan due to this effect. 

From the above considerations it was assumed that the effects of bad 
scanning were small, and could be neglected in any analysis of the results. 


(6) Track Measurement 

A determination was made of the scattering parameter, %, and of the 
grain density, g, on each of the 345 tracks found in a scanning length of 
50cm. The scattering parameter was determined by the coordinate 
method of Fowler (1950) using a 4D block cut-off, while the grain density 
was determined using an arbitrary counting convention. In those cases 
where the grain density was greater than about eight times plateau and the 
particles were identified as being of charge two or greater, the delta ray 
density per 100 u was also determined, using the same counting conventions 
as those employed by Dainton et al. 

Due to the difficulty of making accurate and reproducible grain counts 
on tracks having densities greater than about 50 grains per 50 y, this para- 
meter could not be considered to be as reliable as the scattering parameter. 
While it was possible to resolve the doubly charged from the singly charged 
particles (fig. 1), it was not possible to draw any conclusions concerning 
the existence, or otherwise, of a ‘trough’ in the grain density of the 
relativistic «-particles, analogous to that occurring for the fast singly 
charged particles (Morrish 1952). This is in agreement with the results 
of Zorn (1953) who found that even when more refined grain counting 
techniques are used, there is no evidence for such a trough. Moreover, 
the question of whether there are any He?® particles among the primary 
a-particles must await better measurements. 

Since it was possible in nearly every case to distinguish by inspection 
between particles with charge one, and those with charge two or greater , 
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the measurements on the singly charged particles were made with less 
accuracy than those on the multiply charged particles. In those few 
cases where this immediate identification was not possible the parameters 
a« and g were measured as accurately as were those for an «-particle. As 
a result of this simplification, little reliance should be placed on the 
individual measurements of any of the singly charged particles. 


Fig. 1 


GRAIN DENSITY - GRAINS PER 50m 


SCATTERING PARAMETER o IN DEGREES PER [OOm 


Relation between the grain density—in grains per 50 »—and the scattering 
parameter a—in degrees per 100 py. Star producing a-particles are 
represented by the symbol @, all: other particles by O. The proton (P) 
curve was calculated from the range-energy relation, and from 
comparison grain counts made on -mesons ending in the emulsions. The 
deuteron (D) and triton (T) curves are constructed from the proton curve. 


§ 3. THe Mean Free Pata 


In calculating the flux at the top of the atmosphere it is necessary to 
know the value of the collision mean free path of «-particles in air (Asir) and 
glass (Aq). Unfortunately a direct experimental determination of Aair 
does not appear to have been made, and as a result it is necessary to obtain 
this quantity either by extrapolation from the values found in other media 
(Bradt and Peters 1950), or from the variation in star density with depth 
in the atmosphere (Segre 1952). The results obtained by these methods do 
not agree, and both are open to criticism. The extrapolation of Bradt 
and Peters is empirical, while the results of Segre could be strongly 
influenced by the difficulties of detecting a-particle produced stars. 
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In an attempt to distinguish between these various conflicting results, a 
direct determination has been made of the collision mean free path of 
a-particles in emulsion (A,,). In making this determination it has been 
necessary to neglect any variation of A,,, with energy, due to the limited 
experimental material available. The average energy may be taken as 
being of the order of 1 Bev per nucleon. 

In the various emulsions scanned a total length of 165-3 cm of «-particle 
track was observed, but due to the scanning criteria used, this length 
had to be corrected before the mean free path could be calculated. From 
the track length of each «-particle it was necessary to subtract a quantity 
¢@ which was either (i) the track length traversed before reaching the 
scan (3) or (ii) 8mm if 6 was less than 8mm, the minimum length of 
a-particle accepted. Because of this correction the track length of 
a-particles which were capable, under the scanning conditions, of pro- 
ducing an observable nuclear disintegration was reduced to 80-1 cm. 
Since all the plates scanned were exposed with another plate facing them, 
it was in certain cases possible to follow the «-particle tracks found in 
one plate into the facing one. This resulted in a further 21-4cm of 
track. Altogether six nuclear disintegrations were observed, four in the 
scanned plates, and two in the facing ones. 

Due to the poor statistical weight of this result, a supplementary 
investigation was undertaken to increase the significance of the value — 
found experimentally for A,,,. 

Using scanning criteria similar to those employed previously, the tracks 
of «-particles were found near the top edges of various plates of a 
stack of 600 uw Ilford G5 stripped emulsions, which had been flown at an 
altitude of about 80000 feet, at a geomagnetic latitude of 55°N. These 
tracks were identified by scattering and grain density measurements as 
being those of «-particles, and were followed through the stack until they 
either caused a disintegration or left it. A similar correction was applied 
to the lengths of the tracks as before. In addition a number of high 
energy heavy primary disintegrations were found from which relativistic 
a-particles were emitted. These were also followed through the stack. 

Altogether, in a total effective path length of 199-2 cm, thirteen dis- 
integrations were observed, which gave a final value for Aaa 


Aem*=60+17 g/cm? 


=15-3-+4-3 em of emulsion, 


- 


which does not agree with the experimental value of 120 g/cm? found 
by Segre. 

In order to convert this experimental value for the mean free path 
in emulsion to a value for the mean free path in air, it is necessary to 
find a relation between the mean free path and the type of absorber 


* It is interesting to note here that Zorn (1953), using a monoenergetic beam 


of a-particles of total energy 380 MeV, finds a collision mean free path of 
80+10 g/cem?, 
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through which the particle is passing. An empirical relation of this kind 
has been found by Bradt and Peters (1950 b) from a consideration of the 
mean free paths of particles with Z>2 in brass and glass. In this the 
cross section o, is represented by the relation 
op=n(rytre—24r), 
where the subscripts refer to the two colliding nuclei, 
r= 1-45 x 10-2 x A¥3 cm 
and 
Ar=0-85 x 10-38 cm. 

Using this relation, A,,, is found to have a value of 72 g/cm?, which 
is in good agreement with the experimental value found in this investigation. 
As a result of this agreement the values for A, (44 g/cm?) and Ay 


(53 g/em2) used in the calculation of the flux have been those deter- 
mined by Bradt and Peters. 


§ 4. THe Primary «-ParTIcLE FLux 


In order to obtain the flux at the top of the atmosphere it was necessary 
to correct the directly observed flux for the effects produced by the matter 
above the plates. To reduce the uncertainties introduced by these 
corrections as far as possible without unduly increasing the statistical. 
errors, particles entering the plates at zenith angles greater than 60° 
were excluded from the analysis. Calculation showed that for the re- 
maining «-particles the number brought to rest due to the loss of energy 
by ionization was small and could be neglected. 

The flux at the top of the atmosphere (I) in particles per m? per 
steradian per sec was then given by the expression 


r{azser(2)G)()-s}o 15+ 8) 


where N=the total number of «-particle tracks crossing a scan (104. 
ef. 82 individual «-particle tracks), 6,=the zenith angle of incidence of 
the ith track, 30/27—a geometrical factor, h=the time of exposure 
(4 hours) in seconds, a=the area of emulsion scanned (2 cm?) in m?, 
2 ,—=the average grammage traversed by an «-particle in air passing from 
the top of the atmosphere to the stack (18-5 g/cm”), «,=the average 
grammage traversed by an a-particle in glass in the stack before being 
detected in a scan (2-23 g/cm?), d=the average value of ¢ as defined in 
§3 (4:10 g/cm?) and z=the number of secondary a-particles per m? per 
steradian per sec entering the plates due to nuclear disintegrations 
produced by heavy nuclei. 

In order to calculate z it was necessary to assume values for the flux 
and mean free paths of the heavy primaries, and for the frequency of 
production of «-particles in nuclear disintegrations initiated by heavy 
primaries. The values used were those of Dainton et al. for the flux, 
those of Gottstein for the frequency of production of «-particles, and those 
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of Bradt and Peters for the mean free paths. The effect of re-absorption 
of these secondary «-particles was calculated making the simplifying 
assumption that on the average they were created mid-way between the 
top of the atmosphere and the plates. This factor was so small that a 
more refined calculation appeared unjustified. 

The influence of any secondary «-particles due to disintegrations initiated 
by singly charged particles was neglected. That this was justifiable 
could be seen from the data of Sdérensen (1951) on the frequency of 
ejection and distribution in energy of the helium nuclei emitted in 
disintegrations of this type. Even if the extreme assumptions were 
made that the probability of ejection of helium nuclei in air was the same 
as in emulsion, and the ionization loss above the plates could be neglected, 
it appeared that secondary «-particles of this type would not make up 
more than 5% of the total number of secondaries entering the plates. 
Since z was only about 8°% of J, the influence of the charge one initiated 
a-particles could be neglected. 

The value found for the flux was 

I[=320-+36 «-particles per m? per steradian per sec. 
The error quoted in this value is purely statistical and does not include 
the uncertainties introduced by the mean free path corrections, which 
will however be relatively small as the value of J is not very sensitive 
to small changes in the mean free paths. 

This value should be compared with the value of 280+8 «-particles 
per m? per steradian per sec, found by Ney and Thon (1951) at the same 
geomagnetic latitude and an altitude of 105000 feet using a scintillation 
counter inserted into a counter telescope. 

Four tracks were found of particles having charges greater than two, 
all of the Li, Be, B group, but due to the very limited statistics it did 
not seem profitable to give a value for the flux. 


§ 5. Tur ENERGY SPECTRUM 
(a) Introduction 

Kaplon et al. (1952) have obtained integral energy spectra for the 
various groups of nuclei composing the heavy (Z>2) component of the 
primary cosmic radiation. They were able to relate these spectra to an 
empirical expression of the form 
eG 
~ (1+ 2)" 
where N(H) was the number of «-particles having an energy greater than 
H, E was the kinetic energy in Bev per nucleon, n was a constant having 
the value of 1-35, and C was another constant whose value depended 
on the group of nuclei being considered. For the a-particles C was 
given as 400 particles per m? per steradian per sec. 

Dainton et al. and Gottstein in their work on the heavy primaries 
(Z>3) have reported energy spectra which fall off with increasing energy 
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more rapidly than that of Kaplon eé al. but which like that of Kaplon 
et al. appear to have a constant value of n when related to a similar 
empirical expression. 


Fig. 2 


INTEGRAL ENERGY SPECTRUM 


2 
&-PARTICLES PER m PER STERADIAN PER SEC 


O14 02 ©3 05 07 10 ee) 
KINETIC ENERGY PER NUCLEON (BEV) 
(a) The integral energy spectrum for the primary «-particles together with the 
statistical fluctuations for each point. For comparison the spectrum of 
Be, B, C, N, O, F, nuclei (x 10) found by Dainton ef al. is shown marked 
by X. Curves with n=1-9 and n=1:35 are also shown. 
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In view of the disagreement concerning the form of the energy spectrum, 
and the generally accepted constancy of the form for the various com- 
ponents of the heavy primaries, the energy spectrum of the «-particles 
is of considerable interest, since a determination of it should help to 
decide between these two conflicting reports. 


(b) Determination 


The values of the scattering parameters obtained gave the energies 
in the plates up to about 3 Bev per nucleon. At higher energies than this 
the measurements became unreliable as the determinations of the scattering 
parameters could only represent upper limits. In order to find the 
primary energy distribution it was necessary to extrapolate these individual 
energies to the top of the atmosphere. To reduce the uncertainties 
involved in this extrapolation, particles with zenith angles greater than 
60° were excluded from the analysis. 

The energy loss due to ionization as the particle moved from the top 
of the atmosphere to the detector was evaluated for each individual 
particle, the assumption being made that each particle was a primary 
one. Thus it was possible to obtain the energy distribution at the top 
of the atmosphere, and the resulting integral and differential energy 
spectra are shown in fig. 2. 

On the integral energy spectrum are shown curves given by the 
empirical expression of Kaplon et al. with a value for C of 505 particles 
per m? per steradian per sec, and with n=1:35 and 1:9. Although 
the curves are not very sensitive, particularly at low energies, to the 
precise value of 7, it appears that the experimentally obtained spectrum 
is best fitted when n=1-9. While Dainton e¢ al. and Gottstein have not 
expressed their energy spectra in this form it would appear that they 
are in good agreement with this result, as is shown in fig. 2, where the 
results of Dainton et al. for the heavy nuclei (4<Z<9) are compared 
with those obtained here. 

It would thus appear that, provided it is justifiable to assume that 7 is 
a constant from Z=2 to higher values of Z, the energy distributions 
obtained in the present experiment are in agreement with some other 
workers. However it does not appear possible to reconcile these with 
those of Kaplon et al. Further investigations into the form of the energy 
spectrum at low energies are required in order to resolve this discrepancy. 
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CXLVILL. Neutron-Widths and the Density of Nuclear Levels 


By J. M. C. Scorr 
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[Received July 7, 1954] 


SUMMARY 


The theoretical relation between the level width and the mean spacing 
D of the levels is investigated, using a model with strong configuration- 
mixing, not strong enough however to destroy the size resonance effect. 
Rough numerical values of J/D are obtained using a previous analysis 
of the size resonance effect in the thermal-neutron cross sections. The 
lack of sharpness of the surface of the nucleus has been taken into account, 
and is responsible for a factor of two or three in the calculated level widths. 


$1. INTRODUCTION 


In the present state of nuclear theory, the quantitative interpretation 
of resonance reactions can be reduced to the problem of explaining or 
predicting the partial widths of the energy levels of the compound nucleus. 
Although in favourable cases, with light nuclei, it may be possible to make 
a theoretical estimate of the width of an individual level for nucleon 
emission, it seems that usually we can only hope to predict the average 
width in a statistical sense. 

In view of the extensive data on level widths which it is hoped will 
soon emerge from experiments with fast choppers of the Brookhaven 
type, it seems desirable to consider what should be expected theoretically, 
taking into account the size-resonance in nuclei, in order that we may make 
the best use of the new observations in testing theories of nuclear structure. 
The model developed here, which may perhaps be called the medium 
interaction model, involves no new adjustable parameters, so that a 


comparison with observed level-widths ought to furnish quite a severe 
test. 


§2. WipTHS AND Mzran Spacine or LEVELS 


There is certainly a rough proportionality between the mean level- 
spacing, D, and the mean reduced level-width (when barrier-penetra- 
bility has been allowed for). The method of deriving this relation from 
thermodynamic reasoning (the evaporation theory of Weisskopf (1937), 
N. Bohr and Frenkel) does not seem capable of much further refinement. 
A more promising line of attack is suggested by the refraction picture 
(Feshbach, Peaslee and Weisskopf 1947, Blatt and Weisskopf 1952, 


* Communicated by the Author, 


a ae a 
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Kapur and Peierls 1938). In Weisskopf’s later treatment. it is supposed 
that one can write a wave function for an incident neutron, just inside 
the surface of the nucleus, in the form 


UneCir) cos pMra-C\ie Fay... 6. (1) 
and its logarithmic derivative at the surface is then expressed in the form 
fA/\=Rl(ré) /rb),-p=— KR tan [KR+¢(F)]=—KR tan z(£). (2) 


The angle z(#) which is thus introduced is monotonic and is equal to a 
multiple of 7 at each resonance. 

If we are prepared to assume that the average value of dz/dH at the 
resonances is equal to its general average, 7/D, then the usual formula 
for the average neutron-width of the levels follows, namely 


2D KAM arti gs a. (3) 


This assumption is sometimes described by saying that the behaviour 
of z(#) for an actual nucleus resembles the smooth behaviour that would 
be expected for a simple refracting sphere, but with a much longer 
optical path-length inside. 

In order to get a more precise relation one needs a more detailed 
hypothesis about nuclear structure. It is not difficult to show that with 
some models (3) can be quite misleading. For example, consider for a 
moment the most extreme weak-interaction models, which are used in 
the theory of atomic spectra. In such a model the wave function of a 
state of the compound nucleus is supposed to be quite well described by 
a zero-order approximation constructed out of single-particle wave 
functions, or else out of a one-particle wave function and a wave function 
for the target or core. Most of the levels will only be able to emit a 
neutron by virtue of configuration mixing, and will have widths much 
smaller than (3). These are the ‘ compound levels ’ of Bohr and Mottel- 
son (1953) or the ‘ Auger levels ’ of Scott (1954 a). If such a model were 
correct, the behaviour of Weisskopf’s phase angle 2(#) would deviate 
greatly from the smooth behaviour originally suggested, in a way depend- 
ing on the optical path length to the centre and back. 

In such a weak-interaction model, the perturbation term in the 
Hamiltonian must be so small that it shifts the levels of a configuration 
by an amount small compared with the distance of the next configuration 
with which it can mix. This criterion is practically never satisfied in 
nuclei. The level spacing revealed in neutron transmission experiments 
is often of the order of 1 kev or less, whereas the perturbing term is 
certainly far larger than would be needed to displace the levels from their 
zero-order positions by 1 kev. A statistical theory of perturbations is 
developed below, for the case in which the level density is great but the 
perturbing term not large enough to destroy the influence of the ‘ optical 
path’ 2KR on the behaviour of 2(E£), producing deviations from (3) 
corresponding to systematic fluctuations in dz/dE. 
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§ 3. WavE-FuNCTIONS OF COMPOUND STATES IN A 
MepiumM-INTERACTION MODEL 


In order to avoid complication, the following assumptions and simpli- 
fications will be made. A neutron with energy Z is incident on the target 
nucleus. The compound system can only break up in one way, by 
re-emission of a neutron with /=0; thus radiative capture is neglected, 
and inelastic scattering is only possible at higher energies. All angular 
momenta are taken to be zero. 

All co-ordinates describing the internal motions of the target nucleus 
are denoted by the single symbol r;, and the relative position of the 
incident neutron by (r,6,¢). Antisymmetry requirements are ignored 
for simplicity. Following current usage, that part of configuration 
space where the system can be regarded as completely separated, and 
where the wave function can be decomposed 


P(r, r= Xolt:) - (47) *77-*w(r) eek Moke BONS 

is called a channel. 
The Hamiltonian 

H=H®OLAYMLAH'=HO+1H . . . . (5) 


is split up into that of the target nucleus (H™), that of a neutron in a 
suitably chosen central field V(r), and a perturbation term. 
The equation 


HOJ\r,, t)=£;9,," 2 oe 
taken with the boundary condition 
(djor)rD.=-0 at r=; ee 


possesses a complete ortho-normal set of eigenfunctions ®, and resonance 
energies #, for the internal region r<R (the compound nucleus). There 
is another complete ortho-normal set of eigenfunctions ©,, associated 
with the unperturbed zero-order Hamiltonian, and satisfying 


HO® = (HOH) DS, OB OD (0), me (8) 


These two sets are connected with each other by a unitary matrix 
{Cs}, the elements of which determine the level widths’ of the compound 
nucleus (as indicated in § 5 (ii)). We therefore wish to estimate the order 
of magnitude of ¢,,.. 

Now one can easily deduce from the definitions (6) (8) that 


[ O*H'®,,0 dr=(E,—B,) | BP, dr, 
that is to say, 


1 
Cm Bago | OH ©, dre he cet ene 


This resembles the formula for first-order wave functions in perturbation 
theory, but is exact. If the perturbation is very weak and we can assign 


a small upper bound U to the factor on the right, (9) may give the order 
of magnitude of c,,,. But the plausible conclusion that 


leant a U/ |H,—£,, | . . ° e * (10) 
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cannot be right when this upper bound U is large compared with the 
level spacing D. For the completeness relation 

ones corel et ee Says e asp. (11) 
would then give oe 


1= [| ens [P(/D) dB ~ U/D2> 1. 


We should instead expect the following: whereas for U>D each 
perturbed state can be related to a unique unperturbed state (Cy.5mg); 
for U>D any perturbed state ®, consists largely of a mixture of unper- 
turbed states from the range (#,—D, #,+D) in roughly equal proportions, 
so that |c,,,|? is at most of the order of 1/D. A more definite result of 
this kind will now be obtained (§ 4). 


$4. A THEORY OF STRONG PERTURBATIONS 
(‘ Mepium INTERACTION ’) 


The term strong perturbations is perhaps paradoxical, for an effect 
that can be treated as a perturbation must be small compared with 
something: what is assumed here is that the levels are very closely 
spaced (e.g. D=100 ev), and that, whereas the perturbation is much 
too strong merely to displace the levels by an amount <D and instead. 
mixes up many adjacent levels, it is nevertheless weak enough for the 
contamination by distant levels to be slight. Thus, if the perturbation 
is gradually increased from zero to its full value, a perturbed level ®, 
(energy E,) retains hardly anything of the individual characteristics of 
®, except of course its serial number s which is preserved in consequence 
of the non-crossing rule of v. Neumann and Wigner. 

The result we need is not an exact algebraic theorem about the unitary 
matrix {c,,,} connecting the eigenstates of two Hermitian operators 
which differ slightly, but rather involves probability theory in a way 
which makes a physical argument more appropriate. As a convenient 
starting point we may take the analysis given by Weisskopf and Wigner 
for the problem of line-broadening by radiation damping. Just as with 
the radiation field in a large enclosure, H® has many closely spaced 
eigenvalues, and only certain average properties of the eigenfunctions 
are required. 

We analyse the system H-+H’ first in terms of the eigenfunctions 
@ of H, and then in terms of those of H-+-H’. 

Starting from the equation 


ih Ob/Ot=(H + A), See er (2) 
and seeking a solution of the form 
p= Lan (t) Pn \eXP (7 OL ammemae rs (Lo) 
one obtains the equations 
y= (1/18) pot exp [1(L,' —ELy t/t], ans Beene (a) 
Gy= (Ltt) XH om Gm EXP [oC 0) at MOL ae commas (1400) 


The initial conditions are a)(0)=1, with all the other a,,(0)=0. The 
usual text-book procedure is to put a(¢)=1 in (14 @) and to ignore (14 5), 
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which is good enough for small t. Weisskopf and Wigner (1930), however, 
showed how to extend the solution to larger values of ¢. Their procedure 


is to try putting 


ay=exp (—Fgi/n). (15) 
Then by integrating (14 @), : 
| ye 1 PLOSP (He iBy— toil (16) 
a,(t) =F x9 E,—H, +g o) ae 


This satisfies (14 b) too, provided that 4ig is equal to 
ey Pe 
—exp (‘Lo 1B mT xg )t/n aia [Hes pede 
En —Eot+ 3g 
which is assumed to be independent of m, or at any rate constant when 
averaged over any small range of energies. We therefore define g by 
g=2r Av| Hon’ PD. ~~. 

The distribution of the probability of occupation of the states when 

t>h/g follows from (16) 
|H m0. |? 


2 SE 
lanl > EP 
and one can check that 


2 
Elan) P= | ogres UE, “aestas (18) 


1 
| Hom’ P 
™m 


({ > «) 


as it should be. 

In the original application to the emission of light, g gives the natural 
breadth of the emitting state. The expression under the integral sign 
in (18) gives the probability distribution of the energy of the system. 

Now let + be analysed instead into eigenfunctions of H, i.e. the truly 
stationary states of the system. Since the initial condition of the 
system is 

Dy LieygPa, «s.r 

our solution (13) must be 

p=L),0, exp (—iH L/h). 3 hae ale crm 
The probability distribution of the energy of the system can thus be 
expressed alternatively in terms of the local average of | co, [?/D. Hence 
gD 1 
2a Lig 
This relation (21) supersedes (10), and in fact (21) yields ( 10) as a limiting 
case ; (21) is of course now consistent with (11), which (10) was not. 


AY |e, [Pas (21) 


$5. APPLICATION OF THE PERTURBATION THEORY TO 
ELASTIC SCATTERING 


(i) The General Nuclear Dispersion Formula 


For the sake of convenience a derivation of the nuclear dispersion 
formula will be given here with the sumplifying assumptions of § 3 (elastic 
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scattering of neutrons, J=0, etc). Detailed discussions can be found in 
the literature. 

The actual wave function Y at any chosen energy # can be decomposed 
according to (4) in the external region or channel, and it satisfies 


AX ronal, ete ere. ar os (22) 
everywhere. In the internal region it can be expanded in a series 
e200, Merde 1 cea. g kets» (23) 


which may be non-uniformly convergent on the boundary r=R, and in 
fact will be: the derivative w'(R) (cf. eqn. (4)) cannot be obtained by 
term-by-term differentiation, cf. (7). 

In order to bring in the derivative w’(&), one can use (6) and (22) to 
obtain 


i (WHO) *—O,*HV} dry? deo dr=(E,—Eyb, (24) 
When the left-hand side of (24) is simplified by means of Green’s theorem, 
the only non-vanishing part comes from the term 
h? 0" e? 
poe a *_ py ee 
Wi ee rO.*—1r®, aya rv} 
in the integrand ; integrating this by parts, using (4) and (7), the relation 
(24) becomes 


Oa cide ONES hao 


—a — % . — -—= f 2 
s—~ B,—E 2M »/(4n) | J fe eeu pn  *) 
where 
k 9 
= 4/ (477) | [xo dr; dw . . . . . . (26) 


and can be considered as the average of ©, over the entrance (r=) 
of that channel which corresponds to incidence and elastic scattering. 

Inserting the value of Y(r,, R) from (23) and (25) in (4) now gives the 
standard dispersion formula connecting the wave function with the 
reduced widths y,? : 


w(R) Via 
tats “ Ton? meee 2 £00.27) 
hh 
y= aM lds [ 5 . ° ° A 5 (28) 


From this, the Breit-Wigner formula for slow neutron scattering can be 
deduced, and the neutron width 1, identified as T.=2hy?2 which in this 
case has a simple physical interpretation : I’, is proportional to | d, |? 
which is a measure of the probability of finding the neutron in the channel 
entrance, and to the velocity #k/M with which it might escape. 


(ii) Application to the Medium Interaction Model 
It will greatly simplify the analysis if we adopt the assumption of a 
sharp boundary to the nucleus (the necessary correction. can be made 
afterwards), and if we further assume that the reaction radius is not 
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appreciably larger, so that r= represents both the surface of the nucleus 
and the sphere necessary to seal off the internal region from the reaction 
channels. 

The Hamiltonian can be divided into parts as in §3: 

H=HO+LHY +A’, Acne . (5) 
when H is the exact Hamiltonian of the target nucleus, and has eigen- 
functions x; 

Hy (r)=H, x;. 

E,‘) can be conveniently taken to be zero. H™ describes a neutron 
in a suitably chosen potential well 
HY —(h?/2M)V7?+ V(r), 
so that H’ can be treated as a perturbation. 
We now have two complete orthonormal sets of functions in the 
internal region, as in (6), (7) and 8: namely the eigenfunctions of the 
unperturbed problem (H“ +H) which can be written more explicitly 


as 
v,(r) (5) sin K,r 
STE 4 r; - —_- = y. Fi = rr EO A ‘< (29 
Lane im XM) (den) OLR) ee) aac) 
4)7/R (n==O31, Soe <a eee eee 


and the eigenfunctions of the avtual problem (H) which are 
®, ies Cr OR. 


Gn,s 


When @, is analysed in this of into a combination of the y,v,, the 
only terms which contribute anything to its ‘mean’ value ¢, at the 
channel entrance are the terms in y v, (target left in its ground state). 
Thus from (26) 


$= (2/R)¥/2 OB ree ae 
n 


The statistical result (21) then gives, for the average value of the reduced 
width, 


Oe 2MgD ., 1 
UR TER ,-)(K@—-K,)+aPgye (31) 
where K denotes the internal wave Bete [2M(W+£,)]*/2/h, and W is 
the well depth. 

The summation can be performed by means of the formula 

ie tan (€-+-ic) _>¥ 4é0 
Ste nao [(m-+$)?0?— &2+ 07? 46808” 

which can either be proved from Mittag—Leffler’s theorem or derived 
from the series cot z=21/(z—mz). The result is 


AV t= Av ral roel 


(32) 


ye DE oy tan (§+t0) _ DR & sinh 2c 
s Et+ico a €&-+o7 cos 2+ cosh 2c 


2o0~—s sin 2€ 
x {I= sores See ° . ° . (33) 
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Here € and o are positive numbers defined by means of 
Etjo=(Ki+iMgf)2R. . ... . . (34) 
In practice o?<&, so that KR, and the first two factors in (33) 


yield D/7K ; 2é is large enough for the last factor to be neglected. The 
mean level width for neutron emission is then found to be 


2Dk sinh 2o 
7K cosh 2c-+cos 2KR 


which agrees with the Weisskopf formula (3) in the limit o +09. 

It is now necessary however to allow for the effect of the simplifying 
assumption made above, that the nucleus had a sharp boundary. The 
gradual nature of the transition of the field V(r) from the well-depth 
—W to zero can be taken into account as in Scott (1954 b), and has two 
effects : it replaces the factor 1/K by a larger constant S, and it adds a 
constant to the phase angle KR, so that €2KR—«. In the experimental 
reports, the velocity-variation of I’ is often allowed for and the value of 
PM=J4/(B,/E) is quoted, calculated for a standard energy of #,=1 ev. 
With this convention, 

AGT | YASS) sinh 2c 
YD a cosh 2o-+cos 2€’ 
where £=KR—a«a=(0-6A1/3—0:8)m and k,=(2MB,)?/h=2-2 x Lovemimt: 


T=2hy= (35) 


(36) 


§6. NumERIcAL VALUES AND DIscUSSION 


The first factors in (36) correspond to the 2k,/7K of the sharp-boundary 
theory (Blatt and Weisskopf 1952, and eqn. (35)), which would be 
1-4 10-4 with a well depth of 20 mev as used in the literature, or 0-9 x EOss 
with the value 51 Mev which we believe to be more correct (Scott 1954 b). 
With the diffuse boundary however, using S=18-9 x 10-14 cm as in the 
latter paper, the actual value of the factors in (36) is rather larger, namely 
26S. 

The value of the interaction-strength o can be estimated by the method 
given previously (Scott 1954 b), but only roughly because of inadequate 
statistics: for nuclei with 50<N<82, tanh o=0-16-0-04, and for 
N>s2, tanh c=0-7+0-25. This suggests that I/D should have a 
maximum of 4-41 near A=165, a somewhat larger maximum at A=)6, 
and a broad minimum of 0:4-+0-1 at A=102. This is in fair agreement 
with preliminary results from Brookhaven (Harvey, Porter and Hughes 
1954). 

When the widths of enough levels have been measured it may be 
possible to work backwards and deduce o from J. But, in considering 
experimental values of I’,/D, the statistical uncertainties in both D and 
I, when deduced. from a sample of say 5 or 10 levels must be borne in 
mind, as well as the experimental error of an individual I, which often 
depends on an assumed value for the radiation width, so that the task 
of determining I',/D accurately will be difficult and laborious. 
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In certain cases, e.g. Zn, the predicted mean level widths are rather 
insensitive to the choice of a value for oc, and a better prediction is possible. 
For Zn, (36) gives [)/D=3-7x 10-4, with a probable error of ~ 10% 
arising mostly from S. This may be compared with (51:2) x 10-4 from 
four levels in °7Zn observed by Dahlberg and Bollinger (1954). 

It is not possible to compare the present theory with the cloudy 
crystal ball of Feshbach, Porter and Weisskopf, since details of the 
application of the latter to the resonance region are not yet available 
to us. Some of our formulae, e.g. (34) and (35) do indeed suggest an 
optical analogue with an extinction coefficient oc o/R or a complex 
potential W+4ig but it should be noted that o and g do not correspond 
to the neutron being lost from the entrance channel to other channels 
(inelastic scattering or radiative capture). 

It is noteworthy that the argument from the graph of z(#), which 
gave (3), also suggests that when two or three levels happen to be unusually 
close together (#,,,—H,<D), their widths should be below average, 
whereas the present analysis does not lead one to expect any such 
correlation. This point could probably be tested experimentally. It is 
also natural to suppose that the real and imaginary parts of cy, , in (31) 
will follow a Gaussian distribution, and therefore that J° will follow 
dP=(1/I’) exp (—I/TI’) dW, but this does not seem to be strictly deducible 
from the considerations given here. 


It is a pleasure to thank Professor D. J. Hughes for several valuable 
discussions. 


APPENDIX 
Estimation of o from the Thermal Scattering Data 


In this appendix the probability distribution of the coherent scattering 
length will be obtained. It is a consequence of the present theory 
which was quoted without proof in Scott (1954 b) (referred to as IT). 
According to IT, §3, this problem is equivalent to finding the distribution 
of the values of 8, the phase of the neutron wave function when extra- 
polated inwards to the centre. 

The relation between the notation of II and that of Blatt and Weisskopf 
(1952) is as follows : 


II=S8. (1954 b) |B. & W. (1952) 


Phase just outside —60 2(0) —4a 
Phase extrapolated to centre £(0) —40 
Difference of these KR-a=é KR 
Phase shift in traversing nuclear surface &—a7 0 


The result of the present paper can be interpreted in terms of the graphs 
of the Feshbach-Weisskopf functions 2(Z) and ¢(Z), and their systematic 
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fluctuations in slope. These fluctuations of average slope are periodic in 
C=z—fE=2z—KR-+«. The results (33) and (35) gave the average value 
(y?/S) of dH/dz at those points where z is a multiple of 7 (i.e. the resonances) 


eee oS ra a a a 
di ~= dz ~= cos 2¢-++cosh 20° 
Now if the energy of a thermal neutron, #0, falls randomly among 
the resonance levels, the corresponding value of ¢, namely ¢(0)=$+ 37, 
does not fall randomly on the scale of £: instead, the chance that it hes 
in a range df is proportional to d#/d¢, and the probability distribution 
of 8 is therefore 
sinh 2c dp 
a cosh 20—cos 26° 
This is the formula which was used in II to estimate o. 
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CXLIX. The Characteristic Temperature of Platinum from X-Ray 
Reflections 


By K. Atrexopouxtos and P. EutHymrov 
Physical Laboratory, University of Athens * 


[Revised draft received August 10, 1954] 


SUMMARY 


The characteristic temperature 6 of solids is generally derived from 
measurements of the specific heat. It can also be derived from the 
temperature factor of the diffraction of x-rays. The first method gives 
uncertain results in the case of the transition metals, because in these the 
electronic heat constitutes a considerable part of the total specific heat, and 
therefore corrections have to be made in order to obtain the specific heat 
due to the lattice vibrations alone. In the present paper the intensity of 
the Debye-Scherrer lines diffracted from platinum were determined at 
different temperatures. The gradual change of the ratio of intensities of 
two lines gives the desired value of 6. 

The intensity of the diffraction lines was determined from the density of 
blackening on a photographic film. The temperatures were determined 
from the expansion of the lattice constant. 

Exposures were made at 13 different temperatures between 298° and 
1133°K. From the ratio of the integrated intensities of Debye—Scherrer 
lines (111) and (311) a value of 0=152+-9° results. From the pair of lines 
(200) and (311) the more uncertain value of 192+29° resulted. The 
weighted mean of the two measurements after taking systematic errors in 
consideration gives 156+ 26°. 

These results can be compared to values of @ from specific heat measure- 
ments. 

In the 7? range specific heat measurements give 6=233°. At higher 
temperatures the thermal measurements have to be corrected for electronic 
heat. Accepting for the degeneration temperature 7')=1750°K and for 
the number of electrons g=0-18, values around §=235°K result. The 
6 values from x-ray diffraction are thus decidedly lower than the thermal 
values. The same disagreement has been observed with the metals Li 
and Na, 


$1. IvrrRopUCcTION 


ACCORDING to the Debye theory the thermal vibrations of solids have a 
spectrum that ends abruptly at a maximum frequency, v,. The ex- 
pression fy,,/k is called the Debye characteristic temperature where h and 


k are Planck’s constant and Boltzmann’s constant respectively. Once 


eerie ee OR 
* Communicated by the Authors. 
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8 is known the above theory permits the calculation of the specific heat C, 
of the solid (from now on called lattice heat) at all temperatures. Con- 
versely measurements of the specific heat gives the value of the parameter 
6. According to the above theory, the values of @ obtained from these 
measurements should be constant at all temperatures. For a large 
number of solids this has indeed been found to be the case but in certain 
interesting instances @ varies. Blackman (1937) has given convincing 
reasons for these deviations. Another point that complicates the 
correlation between specific heat measurements and _ characteristic 
temperature is the existence of an electronic specific heat. This is, for 
most solids, small in relation to the lattice heat (except at temperatures 
near the absolute zero) so that this complication is not important. 
The transition metals however, in which the density of the 3d-states at the 
highest occupied level is high, is an exception to the above rule and will 
have considerable electronic heat. It is clear that to calculate the 
characteristic temperature of the lattice vibrations in this last category of 
elements, we must first subtract the electronic heat from the measured 
specific heat. The residue is due exclusively to the lattice heat and can be 
used to calculate the characteristic temperature 6. This has been success- 
fully done by Clusius for Pd (1947), Ni, Fe, Co (1952 a) and in (1952 b). 

Another way of investigating the lattice vibrations is based on a com- 
pletely different principle : The reflection of x-rays on a lattice plane at an 
appropriate Bragg angle, assumes that all the atoms are at their equili- 
brium position, i.e. on the reflecting lattice plane. As temperature rises 
the heat waves will constitute a perturbation and the resulting reflection 
will be less intense. Debye and Waller calculated the change of the 
integrated intensity J of a diffraction line against temperature and 
correlated it with a characteristic temperature 0, according to the 
formula [=J, exp (—2M), where 

6h? sin? 3 /d(z) 1 
tian ( a ;) 

This calculation was made for an isotropic simple cubic lattice but 
experiments showed that it is also valid for a large number of other 
lattices. Thus James and Brindley (1928, 1929) and others found that it 
could be applied to NaCl, KCl and Al throughout a wide range of temper- 
atures. The values of 6y that resulted from their measurements practi- 
cally coincided with the values of 6 obtained from specific heat measure- 
ments. X-ray measurements have also been made by Andriesen (1935) 
on Ag and by Owen and Williams (1947) on Cu, Au, and Al in the range 
293-920°x. The accuracy obtained in the measurements by Owen and 
Williams was sufficient to detect the small gradual change in the value 
of 0,, due to the thermal expansion of the lattice. In all these cases, the 
characteristic temperature derived in this way was in good agreement with 
the characteristic temperature derived from specific heat measurements. 

On the other hand, other measurements are much less satisfactory. 
Arakatsu and Scherrer (1930) applied the above formula to the diffraction 
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lines of Liat room temperature and obtained 0;,=510°, although a value 
around 350 can be shown to fit their results as well. A more sensitive 
approach was undertaken by one of the authors (Alexopoulos 1932), who 
compared the intensity of lines at two different temperatures, 90° and 
293°K, obtaining 0,;,=355-+18°. Pankow (1936) repeated the same 
measurement and got 352--12°. He also compared the intensities in 
smaller temperature intervals and obtained 0,,=315+11° in the interval 
90-190°K, and 0,,=378-+16° in the interval 190-293°K. These results 
show that 0,, gets larger with increasing temperature and 0,,=352° should 
be considered as a mean for the two other values. These results should be 
compared to the values of @ derived from specific heat measurements. 
Simon and Swain (1935) show that the curve C,,=/(7’) cannot be described 
with a constant value of @. At very low temperatures (7’=15°K) @ is 
approximately 330° and it gets larger as temperature rises to reach 430° at 
190°k from where on it stays practically constant. In the range of the 
X-ray measurements it has a mean value of about 410° in the interval of 
90° to 190°K and 430° in the interval from 190° to 293°K. 

A second similar case is Na. Dawton (1937) using X-rays investigated 
the region between 117 and 291°kK and found @;,=140° at the lowest 
temperature and a gradual drop to @;,=120° at room temperature. These 
values are all smaller than 6=159° that results from specific heat measure- 
ments (Simon and Zeidler 1926) in the temperature range between 20 and 
200°K. 

The present paper reports on the investigation of Pt in the range between 
298° and 1133°K. For this transition metal the total specific heat consider- 
ably exceeds, at higher temperatures, the classical value of 5, 96 cal/deg. 
mol.,indicatingalarge contribution from the electrons. Hence, before using 
specific heat values for the determination of @, an appropriate correction 
for the electronic heat has to be made. 


§ 2. DeTERMINATION OF 0 FROM DiFFRACTION EXPERIMENTS AT DIFFERENT 
TEMPERATURES 


The power of a beam of x-rays diffracted under a Bragg angle 9 from a 
crystalline powder at a temperature 7’ is 
le*r8n?j dV 1-+-cos? 29 


P=[,——_—- ——_-_—_ 
9 167rm2ct sin 9. sin 29 


XA(A, 9). [ZF exp {2ri(ha,+hy,+lz,)}?.exp(—2M), . . (1) 


where J) is the intensity of the incident x-radiation ; e¢, m, ¢, A, aV (h, kb) 
have their usual significance, J is the length of the Debye-Scherrer line, 
j is the multiplicity factor for the crystal planes (h, k, 1) corresponding to 
the glancing angle %, » is the number of cells per unit volume, r is the 
radius of the Debye-Scherrer camera, A is the absorption factor depending 
on A and 9 and I, is the scattering factor of an atom k at rest whose 
reduced coordinates in the unit cell are x,, y,, 2 ke 
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If we compare the powers of two different lines at glancing angles 9, and 
% taken at the same temperature 7’, but with different wavelengths K , 
and K,, and calculate the logarithm of their ratio, L, we have 


= 12h? /sin? 9, sin?9,\ 12h2 

ens aol Werte: )-ae 
sin? 3, sin? S2\ $(8u/7) 5 
Ae? pe 2 6y/T ? Ole eg (2) 


where C is a constant practically independent of 7. Only the third term 
contains the temperature and is an apparently complicated function of 7’. 
Nevertheless if we limit the temperature range from 300°K to 1000°K and 
keep the values of the parameter 0,, between 100° and 250°, the function 
T/0y° .d(0y/T') plotted against temperature, gives a family of curves 
which are very approximately rectilinear. Their slopes depend on the 
respective value of the parameter 0). 

According to the above calculation, if the ratio of the intensity of two 
diffraction lines is measured at different temperatures and is plotted on a 
semilogarithmic scale against 7’ a straight line should result with a slope 
equal to the slope of the third term of eqn. (2). Thus by comparing the 
slope of the experimental line L=/f(7’) with that of the corresponding 
theoretical line we get the value of 6. 


§3. EXPERIMENTAL PART 


The power of the Debye—Scherrer lines diffracted from platinum was 
measured by the blackening on a photographic film. 

The scattering specimen was a vertical platinum wire of 200 » diameter 
(fig. 1) continuously rotated by means of a pulley. The rotation was 
necessary because at high temperatures a certain amount of recrystalli- 
zation occurred resulting in spottiness of the Debye-Scherrer lines. The 
wire was kept taut by a small weight which dipped into a container with 
mercury. This served at the same time as the second lead for the electric 
current that heated the wire. Water circulating in the brass body kept it 
cool. The film was pressed around the camera’s external surface by means 
of a thin brass sheet covered with black velvet. Elastic tape held every- 
thing in position. The camera was exposed to unfiltered iron radiation, 
working at a voltage of 40 kv anda current of 5 ma. The duration of each 
exposure was | hour. 

The temperature, 7’, of the wire was determined from the current a that 
flowed through the wire. In a first series of measurements (Euthymiou 
1952)the calibrating curve 7'=f(7) was determined experimentally bymeans 
of a thermoelectric method: a junction of the platinum wire with a thin 
constantan wire (40) was heated in a suitable oven to various temper- 
atures, and the resulting thermoelectric force was measured, Then the 
wire was heated under normal conditions in the scattering chamber and 
the thermoelectric power measured against heating current. In this way 
the relation 7'=f(i) was obtained. This method was found to introduce 


SER. 7, VOL. 45, NO. 37I.—DEC. 1954 48 


1336 K. Alexopoulos and P. Euthymiou on the Characteristic 


serious systematic errors, as the temperature at the junction of the heated 
platinum wire was found to be considerably lower than along the rest of its 
length, where the scattering volume was situated. This method was then 
abandoned in favour of a direct determination of the temperature from the 
expansion of the lattice constant. Owen and Yates’ measurements (1934) 


Fig. 1 
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of the lattice constants of platinum at different temperatures, permitted 
us to calculate the distance of two suitable Debye-Scherrer lines on the 
photograph and to plot distance against temperature. The temperature 
for each value of the heating current i was found by measuring the distance 
s on the corresponding film. Four exposures were made at each value of 7 
resulting in eight values of s. Because of errors in the measurements of s on 
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a comparator, there was a large scattering in the temperature values thus 
determined. Nevertheless, this is not very troublesome, as the method is, 
at least, exempt from systematic errors. The radius of the cylindrical . 
film does not enter into the calculation since exposures were made at room 
temperature, and the distances s referred to them. 

Exposures were taken at 13 different heating currents in the range from 
room temperature to about 1133°K. No higher temperature could be 
attained because the rate of recrystallization became excessive. At each 
value of the heating current four exposures were made, each giving two 
semi-circumferential films, so that a total of eight films could be used for 
each temperature. The density curve of each Debye—Scherrer line was 
measured with a microphotometer. The lines (111), (200), (220), (311), 
produced by the K, radiation and the line (331) produced by the Kg 


Fig. 2 
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radiation were measured. In order to find the relation between the 
indications of the galvanometer and the quantity of radiation (intensity sé 
time), a number of films were exposed to direct blackening with the aid of a 
rotating, graded sector. With the aid of such a calibrating curve the 
density curves for each Debye-Scherrer line were reduced to intensity 
curves and then their areas were measured with a planimeter. Each area 
represents the integrated intensity I of the corresponding Debye—Scherrer 
line, which is proportional to the power of the diffracted x-radiation. 
Figure 2 is a plot for lines (111) and (311), where the abscissa is temper- 
ature, and the ordinate is In (I441/Ls11) calculated from the mean values on 
eight films. The slope of the most probable straight line through these 
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points was found by the method of least squares. Each point was given a 
different weight depending on the divergence of the eight values from the 
mean value. By comparing this slope with the slopes of the calculated 
family of straight lines (§ 2) we found the value of 6y=152°. The probable 
error in the slope was found to produce a probable error of 4@y—=-+9°. 
This value of 40, results from fig. 2 if the values of the temperature are 
considered without error. In reality errors exist because the temperature 
was computed, as mentioned, from the ‘ most probable curve ’7=f(T') that 
passes through scattered points of the lattice expansion measurements. 
On both sides of this curve two curves were traced so as to contain all the 
scattered points. The distance of these from the ‘most probable curve ’ 
gives indication of the maximum possible error in temperature. At 
1130°K this error is about 120°. If the temperature of 1130°K is thus 
changed to 1250°K, and if the other intermediate temperatures are 
correspondingly corrected, new values of 6y result. By considering in 
this way the maximum possible systematic error in the temperature 
measurement, 6,, could either be larger or smaller by 8°. In fig. 3 we 


Fig. 3 


plot experimental results for the second pair of lines (200) and (311). The 
same calculations as above were made for this pair of lines. All the results 
are grouped in the table. We have also plotted the results for the other 
three pairs of lines (220)-(311), (222)-(311) and (331)-(311), but the 
points were so broadly scattered and the resulting error so large that 
results would be without significance, 

The two measurements of the table after being given appropriate weights 
give the mean value 4,,=156°-- 23°. If the systematic error is also taken 
in consideration, the total error reaches 26°. 
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$4. Discussion oF RESULTS 


In §3 we showed that the change of the intensity of x-rays diffracted 
from platinum leads to a characteristic temperature of 156°. This value 
will now be compared to values of 6 found in other ways. The most 
important methods are specific heat measurements in the 7? range and 
in the medium temperature range up to 70. 


Maximum 
Lines 6% | Random error systematic 
error 


(111)-(311) 152 +9 +8 
(200)-(311) 192 29) +10-6 


The specific heat of a solid is given by two terms, the first due to the 
lattice and the second to the free electrons. The first term is given by 
C,, lattice=6R/2 x D(T/@) where D is the Debye function. At very low 
temperatures (i.e. Z'<6) this term simplifies to C,,, lattice =446(7/6)°. 

The electronic heat cannot be easily given in an explicit form. At temper- 
atures that are low relatively to the degeneration temperature 7’) (approxi- 
mately up to 7',/3) the electronic heat is a linear function of T 

C= (3) RqT/T=y7, earl act 8 (3) 
where g is the number of free electrons per atom. At higher temperatures 
C, asymptotically takes the classical value C,=(3/2)R . q. 

The curve C,=f(T') in the whole range of temperatures is given in a large 
number of books (e.g. Mott and Jones 1936, p. 180). 

Specific heats at very low temperatures will thus follow a law 

C= 446(T/0)+yT. 
Measurements by Kok and Keesom (1936) in the liquid He-range give 
§@—233° and after a correction by Wohlfarth (1948) one gets T'9/g=7000. 
This latter value in conjunction with the magnetic measurements of 
Hoare and Matthews (1952) gives 7’)=1750° and g=0, 25. In order to 
determine 6 from the specific heat measurements at medium temperatures 
the electronic heat should be subtracted from the experimental C,. 
We notice that if 7’) is equal to 1750°, measurements from about 500°K 
upwards should not be corrected according to the formula (3) but by using 
the exact expression. Since this last point although understood, has not 
been carried out by others exactly, we will undertake here a new evaluation 
of the values of @ at different temperatures. The total specific heat of 
platinum has been measured in the range from 20 to 200°K by Simon and 
Zeidler (1926) and up to 1870°K recently by Jaeger and Rosenbohm (1939). 
Curve Cot of fig. 4 shows the best fit through their measurements. If the 
electronic heat is now subtracted, a curve for the lattice heat only will 


result. 
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By using 7)=1750°K and g=0-18 we obtain a curve for the lattice heat 
that approaches the classical value 5955 cal/deg. mol. fairly well for temper- 
atures up to 800°K. It should be noticed that if the value q=0-25 
proposed by Hoare and Matthews had been used, the agreement would 
have been less satisfactory. 


Fig. 4 
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In the range where the slope is sufficient, the corrected curve can be 
used to calculate 6 at different temperatures. The results are given in fig. 5 
On the same plot the values of @ are shown as obtained by Simon and 
Zeidler from their measurements without correcting for the electronic heat 
We notice that in this range 0 should be given higher values around 235° 
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than the value 225° generally accepted. This gives also a better agreement 
with the value 233° found in the 7? region. On the same figure we insert 
the results of this paper in the form of a shaded area. 

The position of this area shows that the characteristic temperature 0),, 
derived from the temperature factor of x-ray experiments, is decidedly 
lower than the value derived from specific heat measurements. In figs. 2 
and 3 lines with a slope corresponding to 0y=225° have been traced. 
It is obvious that in no case can this value be brought to agree with the 
X-ray results. 

A possible reason for such a discrepancy could be the large difference 
in the temperatures at which the characteristic temperatures were 
determined from specific heat data and from x-ray intensities. The 
resulting change in 6 due to expansion can be accounted for in the 7’'>0 
region by assuming an Hinstein model for the vibrations with a frequency 
equal to k6/h. By using the formula (Mott and Jones, p. 19) 


@ IN Vmax/d IN V=—aVo/x Cy, 


the change of 6 is calculated to dd=—«6(«V>/y,C,) dT. By inserting the 
experimental values (Mott and Jones, p. 318) a change of 9° (from 235 
to 226) results when the temperature is changed from 0 to 700°k—the 
mean temperature of the x-ray experiments. This change is not sufficient 
to explain the far smaller value 0,,;=156° found from the x-ray experiments. 

The possibility of the anharmonic terms of vibration being excessively 
large should also be considered. Since no direct way of calculating these 
terms for Pt exists, the amplitude of vibration of Pt was compared to the ~ 
amplitudes of Al, KCl, Cu, Ag, and Au for which the experimental results 
have been found to be well expressed by a Debye-Waller formula up to a 
quite high limiting temperature. The ratios of mean square displace- 
ments at the limiting temperature to the mean square displacements at 
the melting point have been computed for the cases mentioned and all 
fall in the region 0-55 to 0-67. The same ratio for Pt is 0-44 thus making it 
improbable that the anharmonic terms should be more important in Pt 
than in the other solids. Also the ratio of mean square displacements at 
the limiting temperature to the square of the interatomic distances are 
0-:15% for Pt against 0-23 to 0-48% for the other cases. 

Platinum with Li and Na thus constitutes a third case of disagree- 
ment between 0,,and 6. All three give a 4, value that is too low. 
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The Lattice Thermal Conductivity of Dilute Copper Alloys 
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[Received September 20, 1954] 


Tu lattice thermal conductivity, K,, of metals at low temperatures has 
a particular interest as it appears to be limited at such temperatures by 
‘free electron scattering of the lattice waves or phonons, and therefore 
yields information on, the interaction between phonons and free electrons. 
Apart from measurements on copper-nickel alloys (see, for example, 
Estermann and Zimmerman 1952) and recent work on silver-palladium 
alloys (Kemp et al. 1954), there appears to be little experimental knowledge 
of this interaction. It has,been suggested (e.g. Gurevich 1945, Frederikse 
1953, MacDonald 1954) that the phonon-electron interaction may play 
an important part in explaining the departures of the ‘observed thermo- 
electric forces from those theoretically predicted. 

As part of an investigation of the electrical resistance minimum and 
thermoelectric power, dilute copper-iron alloys having a relatively large 
residual electrical resistance have been produced in this laboratory ; 
these seemed to us ideal subjects for study of the low temperature thermal 
conductivity and subsequent deduction of the lattice component. 

The measurements of thermal and electrical conductivity were made 
between 2 and 150°K on rods of two different alloys in an apparatus 
which will be described more fully elsewhere. The method of measurement 
of the thermal conductivity was similar to that described by White 
(1953 a) and the electrical resistance was measured with a galvanometer 
amplifier (MacDonald 1947). 

The total thermal conductivity is given by the sum of the lattice 
component, K,, and the electronic component, K,, 


Koa Kei, Je toe rear at) 
where the electronic thermal resistivity 1s 
1/K = W = Wot W;,. . . . . . . ° (2) 


W. is the ideal electronic thermal resistivity due to scattering by lattice 
vibrations ; W, is the residual electronic thermal resistivity due to 
impurity scattering and should be related to the residual electrical 
resistance py at a temperature 7’ by the relation 

pol Wol=L, Ne aes) 
where L is the Lorentz number. 


eres ers s 5 2 a re ae Fae eae 
* National Research Laboratories Postdoctorate Fellow, at present on leave 
from Commonwealth Scientific and Industrial Research Organization (Australia). 


1344 Correspondence 


The two alloys investigated were homogenized and annealed speci- 
mens containing nominally 0:056% and 0-0043% iron which’ <5 
residual electrical resistances at helium temperature of 0-56 x 10~® and 
0-041 10-® ohm cm respectively (cf. pure copper at room temperature, 
p=1:72X10-%ohmem). The respective alloys exhibited a rT 
minimum at approximately 25° and 15°K, the resistance at 2°K being 
5 to 8% higher than at the minimum. 
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Experimental results for thermal conductivity of two dilute copper-iron alloys, 
and values deduced for electronic and lattice conductivities. 


At temperatures below 20°x, W,;<W, and hence K=K,+K,. This, 
together with the fact that we expect K,ocT? (e.g. Makinson 1938, 
Klemens 1954) and to be very much less than K o below 4°K, enables us 
to obtain K,/T or W,T' directly from the experimental values. This 
latter quantity, like po, should be constant if we neglect the possible 
effect of an anomaly associated with the resistance minimum. For 
the 0-056%, and the 0-0043°%, alloys we obtain values for Wl" of 21-7 
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and 1-52 cm deg*;watt respectively, which give values for L (eqn. (3)) 
within a few per cent of the theoretical value of 2:45 10-8, 

Values for W, used to calculate the electronic conductivity, K,, at 
temperatures above 20°K have been taken from previous work on 
pure copper (Berman and MacDonald 1952 and White 1953b). This 
work indicates that for 7<50°K, W,; can be well represented by 
71x 10-874 watt/em deg. The calculated curves for K, in the figure 
are then the function 

K,=1/W,=(Wo+T1 x 10-6724)-1 
where W)=21-7/7 for the 0-056% alloy and 1-52/7' for the 0-0043% 
alloy. Values for K, calculated from (1) are shown only for the 0-056% 
alloy in which the lattice contribution is comparable to the electronic 
contribution and hence the calculated results will not be seriously affected 
by small uncertainties in W, or W,. 

In the region below 15°K, where the lattice waves are scattered pre- 
dominantly by free electrons, we expect K, to vary with 7”. The results, 
shown in the figure, are well represented by 

ese 1021 (T=. 16°R), Peer bes. oa) 
This compares with the value of 1-4 10-3 T? calculated by Klemens from 
the ideal electronic conductivity K,, assuming that the free electrons 
interact equally with transverse and longitudinal phonons. 

It is of interest that if K,=1-8x10-°7?, the phonon mean free path 
below 20°K will be approximately given by 2 10-°7’-1cm. The relatively 
strong phonon-electron interaction which we infer from these results 
lends substance to the suggestion of Gurevich that the equilibrium of 
phonons in the lattice may be disturbed by a flow of electrons at low 
temperatures. This may not only affect the thermoelectric power, but 
also the ideal electronic thermal conductivity, there being a marked 
disparity between theory and experiment in both these quantities. 
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(Butterworths Scientific Publications, London.) [Pp. ix+191.] Price 30s. 


Tue study of the related fields covered by the title, after a long period of 
inactivity, has revived in the last decade, and the extensive advances made 
justify this detailed review. Dr. Hall’s monograph is written ina self-contained 
form, including chapters on elementary crystallography and on the * way of 
life’ which has grown up around the experimental investigation of metallic 
single crystals. 

The crystallographic treatment of twinning mostly covers classical material 
(however, some of the twinning elements of titanium listed in Table V and 
referred to on p. 75, seem to have been arrived at by an unsatisfactory method 
of inspired guesswork), but recent advances make up the bulk of the chapter 
on Formation of Twins under Stress. The fragmentary nature of the informa- 
tion here presented is no reflection on the author, but rather serves to emphasize 
that this topic is at the frontier of present knowledge. By contrast, in the 
chapter on annealing twins a much more unified picture is presented. The 
chapter on Theories of Twin Formation is prefaced by an introduction to 
dislocation theory which is, perhaps necessarily, too condensed to be of much 
help to the complete beginner. The dislocation theory of twinning is well 
presented, and the author reflects the present incomplete state of theoretical 
understanding by referring to the Russian elastic treatment of twinning and 
commenting that a synthesis between the two approaches remains to be made. 

The chapter on diffusionless transformations deals largely with crystal- 
lography ; Dr. Hall expounds this involved subject fully and clearly. Purely 
metallurgical aspects, as for instance the stabilization of austenite, are excluded. 
No unified theory of the crystallographic nature of these transformations 
was available to the author when he wrote, and it is a measure of the pace of 
current research that in the past year no fewer than three rival theories have 
been published. 

The book can be confidently recommended as a reliable source of information 
and as a stimulating introduction for research workers entering the field. 
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Theoretical Elasticity. By A. E. GREEN and W. Zprna. (Oxford: University 
Press.) Price 50s. 


InrTEREsT in the theory of elasticity as distinct from the applications of the 
theory has, during recent years, very greatly increased, largely for two reasons. 
The first is the development of complex variable methods in this country by 
Stevenson, Green and others and to the belated discovery of the work of 
Muskhelishvili in Russia. The second is the work of Green, Rivlin and others 
on finite deformations facilitated by the systematic use of the tensor notation. 
These two developments and elastic shell theory are the main topics discussed 
in this very useful and timely book which should be widely read. It is, as the 
title implies, a book on the theory of the subject and is written for those whose 
interests are primarily mathematical. 

_ It is clear and as easy to follow as the nature of the subject permits though, 
in the reviewer's opinion, it would be improved by the addition of a biblio- 
graphy at the end of each chapter. The book is well produced and is at a 


reasonable price and is a worthy successor to the recent volumes of applied 
mathematics from the same publishers. W.M.S. 
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The Theory of Cohesion. By M. A. Jaswon. (Pergamon Press.) [Pp. 245.] 
Price 37s. 6d. 


“ AN outline of the cohesive properties of electrons in atoms, molecules and 
crystals ’’, is the subtitle of this book. It is intended to be an introduction to 
the subject, preparatory to reading the standard treatises; but surely three 
initial chapters on wave mechanics are superfluous nowadays. The various 
approximations for the helium atom and hydrogen molecule are given at 
length, with a chapter on Molecular Orbitals. After a sketch of the theory of 
metals, the results are given of the cellular approximation applied to alkali, 
divalent, and noble metals, and certain alloys. The final chapter deals with 
covalent structures. But the text is very incomplete. For the student, 
there is no general account of types of binding forces (cohesive energy itself 
is only defined in a footnote), so that he will find it difficult to see the relevance 
of each succeeding topic. On the other hand, some modern work of great 
interest, such as the theory of plasma oscillations, is deliberately excluded. 
There is no reference at all to Slater’s MIT reports. This book may provide 
useful summaries of certain topics for research workers, but is too inelegant 
and incoherent to be recommended as an introduction. ve MaZe 
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